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Abstract 

The  optimal  employment  of  autonomous  search  and  destroy  vehicles  is  ad¬ 
dressed.  The  results  apply  to  air,  land,  or  water  vehicles  with  1,  k,  or  infinite 
warheads.  The  specific  scenarios  considered  involve  an  air  vehicle  searching  a  battle 
space  for  stationary  targets  in  the  presence  of  false  targets.  Encounters  are  modelled 
with  uniform,  Poisson,  and  normal  distributions.  Linear  and  circular  search  patterns 
are  examined.  All  relevant  parameters  are  extracted  from  intelligence  information, 
the  sensor  performance  specification,  and  the  air  vehicle  performance  specification. 
Analytic  system  effectiveness  measures  are  derived  using  applied  probability  theory. 
The  effectiveness  measures  derived  in  this  dissertation  handle  time-varying  parame¬ 
ters  which  characterize  the  battle  space  environment  and  the  performance  of  the 
system.  This  allows  the  formulation  and  solution  of  optimization  problems  that 
maximize  the  probability  of  target  attacks  while  at  the  same  time  constraining  the 
probability  of  false  target  attacks.  Optimal  schedules  for  controlling  sensor  threshold 
and  area  coverage  rate  during  a  mission  are  derived  and  compared  to  the  constant- 
parameter  results.  These  schedules  establish  a  system  operating  characteristic.  An 
increase  in  system  effectiveness  is  demonstrated  when  parameters  are  dynamically 
controlled  during  a  mission.  Plots  depicting  sensitivity  to  the  constraint  on  false 
target  attacks  and  sensitivity  to  the  number  of  warheads  are  generated  to  give  deci¬ 
sion  makers  the  complete  trade  space  for  either  designing  new  systems  or  operating 
existing  systems. 
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Establishment  of  a  System  Operating  Characteristic 

for 

Autonomous  Wide  Area  Search  Vehicles 


I.  Introduction 

1 . 1  Motivation 

Carl  Von  Clausewitz  stated,  “War  is  an  extension  of  politics  by  other  means”  [10]. 
At  the  heart  of  wartime  politics  is  blue-force  casualties  and  collateral  damage  to 
non-combatant  entities.  Recent  advances  in  autonomous  target  recognition  and  au¬ 
tonomous  navigation  have  given  political  leaders  the  potential  to  prosecute  targets 
in  dangerous  environments  without  having  to  worry  about  blue-force  casualties.  Col¬ 
lateral  damage,  however,  still  remains  a  concern.  One  can  imagine  the  news  headline 
when  an  autonomous  search  and  destroy  vehicle  takes  out  a  school  bus  filled  with 
innocent  children.  Manned  systems  have  relied  on  humans  as  additional  checks  to 
minimize  collateral  damage.  In  the  case  of  autonomous  systems,  the  human  checks 
are,  by  definition,  removed.  Prior  to  employing  any  such  system,  one  must  determine 
how  well  the  system  can  perform.  System  effectiveness  is  conditioned  on  the  system 
operating  characteristic,  upon  which  the  operating  point  is  defined  by  area  coverage 
rate  and  the  sensor  threshold  used  for  target  declaration.  System  effectiveness  also 
depends  on  the  battle  space  environment  in  terms  of  numbers  and  types  of  targets 
and  false  targets. 

Determining  or  predicting  system  effectiveness  can  be  done  using  mathematical 
analysis,  computer  simulation,  experimentation,  or  operational  flight  test.  This  re¬ 
search  establishes  a  system  operating  characteristic  for  autonomous  wide 
area  search  vehicles  that  can  be  used  to  analytically  calculate  system  ef- 
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fectiveness.  The  presence  of  false  targets  as  well  as  a  sensor’s  likelihood  of  misclas- 
sifying  false  targets  is  acknowledged  up  front.  False  targets  can  be  non-combatant 
entities,  decoys,  or  objects  of  little  strategic  value.  Inputs  to  the  mathematical  “black 
box”  in  Figure  1.1  come  from  the  order  of  battle  intelligence,  the  sensor  performance 
specification,  and  the  air  vehicle  performance  specification.  The  order  of  battle  intel¬ 
ligence  provides  battle  space  dimensions,  target  distributions,  and  false  target  distri¬ 
butions.  Sensor  performance  comes  from  Receiver  Operating  Characteristic  (ROC) 
curves  (determined  experimentally),  which  relate  true  positive  fraction  to  false  posi¬ 
tive  fraction.  Air  vehicle  performance  involves  airspeed  limits  and  endurance  values. 
Outputs  from  the  black  box  are  system  effectiveness,  a  sensor  threshold  schedule, 
and  an  area  coverage  rate  schedule.  Any  objects  reported  as  targets  (correctly  or 
incorrectly)  are  attacked.  The  number  of  false  target  attacks  during  a  mission  is 
denoted  m,  and  the  number  of  target  attacks  is  denoted  n.  System  effectiveness 
can  be  quantified  using  the  probability  of  at  least  m  false  target  attacks,  denoted 
P  [in  >  rh),  and  the  probability  of  at  least  n  target  attacks,  denoted  P  (n>  n ).  The 
values  for  rh  and  h  are  specified  by  the  designer  or  mission  planner.  The  contents  of 
the  black  box  will  be  revealed  in  Chapter  II. 

* 

Figure  1.1  Input-Output  Diagram. 

Clearly  one  wants  the  expected  number  of  target  attacks  to  be  high  and  the 
expected  number  of  false  target  attacks  to  be  low.  Unfortunately,  due  to  the  ROC, 
adjusting  the  sensor  threshold  to  increase  the  number  of  target  attacks  also  increases 
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the  number  of  false  target  attacks.  For  fixed  mission  durations,  lower  area  coverage 
rates  increase  sensor  performance,  but  decrease  total  area  searched.  These  tradeoffs 
must  be  worked  to  find  solutions.  Once  the  internal  definition  of  the  black  box 
in  Figure  1.1  is  known,  optimization  problems  can  be  formulated  and  solved.  The 
general  optimal  control  problem  will  be  of  the  form 

Max:  P  (n  >  h) 
subj  to:  P  ( m  >  m)  <  b 

where  the  upper  bound  b  is  set  by  the  designer  or  mission  planner.  The  control 
variables  are  area  coverage  rate  and  sensor  threshold,  which  are  further  defined  in 
the  next  section. 

1.2  Search  and  Destroy  Mission 

The  search  and  destroy  mission  involves  searching,  detecting,  classifying,  and 
attacking  [1],  This  can  be  done  a  number  of  ways.  For  example,  a  sensor  craft  can 
locate  and  classify  objects,  then  pass  the  information  to  a  separate  bomber  aircraft 
that  performs  the  attacks.  An  Unmanned  Combat  Aerial  Vehicle  (UCAV)  can  locate, 
classify,  and  destroy  objects  on  its  own.  An  expendable  autonomous  wide  area  search 
munition  locates,  classifies,  and  destroys  an  object  as  well  as  itself  in  the  process. 
All  of  the  above  concepts  involve  sensors  with  automatic  target  recognition  (ATR) 
algorithms,  and  warheads.  The  mathematical  foundation  developed  in  this  research 
handles  systems  with  1,  k ,  or  infinite  warheads.  Ironically,  the  infinite- warhead 
case  can  be  used  to  mathematically  evaluate  the  zero-warhead  sensor  craft  case. 
The  relation  comes  from  the  fact  that  a  sensor  craft  broadcasts  target  coordinates 
assuming  a  shooter  exists  to  destroy  the  targets.  A  sensor  craft’s  mission  does  not 
terminate  because  it  has  no  more  warheads.  It  continues  to  search  and  classify  until 
it  is  sent  home  or  shot  down.  Thus  a  sensor  craft  is  mathematically  equivalent  to  a 
vehicle  having  an  infinite  number  of  warheads. 
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1.2.1  Searching.  Searching  requires  a  pattern  and  area  coverage  rate.  Two 
search  patterns  are  addressed  in  this  research.  The  first  pattern  is  the  linear  pattern 
depicted  in  Figure  1.2(a).  Although  Figure  1.2(a)  shows  a  constant  swath  width, 
the  linear  pattern  can  can  be  used  for  any  linear-symmetric  battle  space.  Further, 
any  randomly-shaped  battle  space  can  be  broken  up  into  piecewise  linear-symmetric 
segments.  The  second  pattern  involves  using  concentric  annuli  emanating  from  the 
origin.  Figure  1.2(b)  shows  one  such  annulus  at  radius  r.  Searching  a  disc  of  radius 
R  using  concentric  annuli  emanating  from  the  origin  approximates  an  outward  spiral 
search  pattern  and  is  mathematically  tractable. 


(a)  Linear 


Figure  1.2  Search  Patterns 


Area  coverage  rate  can  be  varied  a  number  of  ways.  For  const  ant- altitude  flight, 
area  coverage  rate  is  the  product  of  velocity  and  swath  width.  Thus,  changing  either 
parameter  changes  the  area  coverage  rate.  Swath  width  is  constrained  by  gimbal 
limits.  For  fixed  degrees  of  azimuth  and  fixed  velocity,  changing  altitude  changes 
the  area  coverage  rate.  In  this  research,  the  flight  condition  in  terms  of  altitude, 
velocity,  and  swath  width  will  all  be  rolled  up  into  area  coverage  rate,  which  will  be 
one  of  the  two  control  variables.  Specific  application  of  whether  one  varies  altitude, 
velocity,  or  swath  width  is  up  to  the  designer  or  mission  planner.  An  example  of 
varying  velocity  to  hold  area  coverage  rate  constant  is  given  in  Section  1.3. 
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1.2.2  Detecting.  Detections  will  occur  based  on  assumed  target  and  false 
probability  distributions.  The  seven  battle  space  scenarios  conceived  by  Jacques  and 
Pachter  [26]  are  examined  in  this  research.  Scenarios  1-4  use  linear-symmetric  battle 
spaces,  and  Scenarios  5-7  use  circular  battle  spaces.  The  differences  in  the  scenarios 
come  from  the  target  and  false  target  probability  distributions  as  well  as  the  number 
of  targets.  For  all  scenarios,  both  targets  and  false  targets  are  assumed  fixed  at 
unknown  locations. 

Scenario  1  consists  of  a  single  target  uniformly  distributed  in  a  linear-symmetric 
battle  space  among  a  Poisson  field  of  false  targets. 

Scenario  2  consists  of  a  linear-symmetric  battle  space  with  a  Poisson  field  of 
targets  and  a  Poisson  field  of  false  targets. 

Scenario  3  consists  of  N  targets  uniformly  distributed  in  a  linear-symmetric 
battle  space  among  a  Poisson  field  of  false  targets.  Note:  Scenario  1  is  a  special  case 
of  Scenario  3  with  N  —  1. 

Scenario  4  consists  of  N  targets  and  M  false  targets  uniformly  distributed  in 
a  linear-symmetric  battle  space. 

Scenario  5  consists  of  N  targets  normally  distributed  in  a  circular  battle  space 
among  a  Poisson  field  of  false  targets.  Normally  distributed  refers  to  a  circular- 
normal  distribution  centered  at  the  origin  with  a  target  standard  deviation  ay. 

Scenario  6  consists  of  N  targets  and  M  false  targets  normally  distributed  in 
a  circular  battle  space.  Normally  distributed  refers  to  circular-normal  distributions 
centered  at  the  origin  with  a  target  standard  deviation  <jy  and  a  false  target  standard 
deviation  ax. 

Scenario  7  consists  of  a  single  target  normally  distributed  in  a  circular  battle 
space  among  a  Poisson  field  of  false  targets.  Note:  Scenario  7  is  a  special  case  of 
Scenario  5  with  N  —  1. 
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Other  combinations  of  target  and  false  target  distributions  exist.  Table  1.1 
shows  where  the  seven  scenarios  fall  in  a  combination  matrix  for  uniform,  Poisson, 
and  normal  distributions.  The  techniques  and  results  in  this  dissertation  could  easily 
be  applied  to  the  remaining  combinations. 


Tab’ 

le  1.1  Scenario  Matrix 

uniform 

targets 

Poisson 

targets 

normal 

targets 

uniform  false  targets 

Scenario  4 

— 

— 

Poisson  false  targets 

Scenarios  1,3 

Scenario  2 

Scenarios  5,7 

normal  false  targets 

— 

— 

Scenario  6 

1.2.3  Classifying.  In  practice,  detected  objects  are  classified  to  a  certain 
level  of  discrimination.  For  example,  one  may  classify  on  object  as  either  air  breath¬ 
ing  or  ballistic.  A  finer  level  of  discrimination  may  be  a  specific  type  of  air-breathing 
or  ballistic  object.  Regardless  of  the  level  of  discrimination,  there  is  a  point  where 
the  sensor  reports  a  detected  object  as  either  a  target  thereby  authorizing  an  attack, 
or  a  false  target  thereby  commanding  no  attack.  Sensor  performance  is  judged  by 
how  often  the  sensor  is  correct.  The  probability  of  a  target  report  is  the  probabil¬ 
ity  the  system  correctly  reports  a  target  when  a  target  is  encountered.  Selecting  a 
sensor  threshold  is  tantamount  to  selecting  the  probability  of  a  target  report.  The 
sensor  can  also  be  wrong.  By  using  a  ROC  curve  to  model  sensor  performance,  the 
probability  of  attacking  a  false  target  can  be  written  in  terms  of  the  probability  of 
a  target  report.  Therefore,  the  second  control  variable  will  be  the  probability  of  a 
target  report. 

1.2.4  Attacking.  Any  objects  reported  as  targets  (correctly  or  incorrectly) 
are  authorized  to  be  attacked.  An  autonomous  munition  would  fly  itself  into  the 
reported  target,  an  autonomous  UCAV  would  launch  one  of  its  warheads,  or  an 
autonomous  sensor  craft  would  broadcast  target  coordinates  to  shooters. 
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1.3  Hypothetical  Operational  Example 

To  help  illustrate  the  relevance  of  this  research,  a  hypothetical  operational 
example  is  given.  An  intelligence  source  states  that  a  strategic  target  is  known 
to  be  in  the  area  designated  Ab  in  Figure  1.3,  which  is  one  of  the  most  heavily 
defended  areas  in  the  region.  It  can  be  assumed  the  target  is  equally  likely  to 
appear  anywhere  in  Ab-  It  is  also  known  that  some  subareas  of  Ab  have  more  false 
targets  than  others.  The  Battle  Commander  tasks  his  analysts  to  provide  options 
involving  an  autonomous  search  and  destroy  vehicle.  He  would  like  the  probability  of 
attacking  the  target  to  be  greater  than  0.50.  However,  he  would  like  the  probability 
of  attacking  any  false  targets  to  be  less  than  0.30. 


Figure  1.3  Hypothetical  Operational  Example.  Map  downloaded  from  [52], 

Based  on  the  order  of  battle  intelligence  given,  the  analysts  use  a  uniform 
distribution  to  model  the  target  encounter.  Since  false  targets  are  often  encountered 
at  certain  rates,  they  use  a  Poisson  distribution  based  on  the  expected  false  target 
densities  (denoted  by  a)  shown  in  Figure  1.4  to  model  false  target  encounters. 
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Although  Ab  is  not  exactly  a  linear-symmetric  battle  space,  the  analysts  as¬ 
sume  the  heading  changes  and  additional  area  searched  are  small  enough  to  make 
a  linear-symmetric  assumption  valid.  Hence,  the  operational  example  resembles 
Scenario  1.  The  sensor  performance  specification  gives  the  ROC  curve  shown  in 
Figure  1.5  for  a  constant  area  coverage  rate  of  5  The  ROC  curve  was  generated 
from  experimental  data  using  various  signal-to-noise  ratios.  True  positive  fraction 
(TPF)  is  the  probability  of  the  sensor  reporting  a  target  when  a  target  is  encoun¬ 
tered.  False  positive  fraction  (FPF)  is  the  probability  of  the  sensor  reporting  a  target 
when  a  false  target  is  encountered.  The  analysts  use  the  curve  fit  in  Figure  1.5  to 
mathematically  relate  TPF  and  FPF. 


Figure  1.5  Receiver  Operating  Characteristic  Curve. 
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It  should  be  noted  that  TPF  has  the  same  definition  as  the  probability  of  a 
target  report.  Thus  scheduling  TPF  is  the  same  as  scheduling  the  probability  of  a 
target  report.  The  sensor  and  ATR  algorithm  were  designed  for  the  constant  area 
coverage  rate  of  5  therefore,  the  analysts  decide  to  fix  area  coverage  rate  and 
use  the  single  ROC  curve  provided.  Adjusting  sensor  threshold  over  time  moves 
the  operating  point  along  the  ROC  curve.  This  amounts  to  a  schedule  for  TPF 
and  a  corresponding  schedule  for  FPF.  Assuming  a  constant  area  coverage  rate  of  5 
the  analysts  must  schedule  airspeed  and  find  the  best  schedule  for  either  TPF. 
The  assumed  path  is  left  to  right;  whereby,  the  vehicle  starts  in  the  area  where 
a  =  4  CPCfa»2 ers  •  Using  differential  equations,  the  analysts  schedule  airspeed  to  keep 
area  coverage  rate  constant  at  5  UU.  No  optimization  is  required.  The  resulting 
schedule  is  shown  in  Figure  1.6.  As  one  would  expect,  higher  velocities  are  needed 
at  the  narrow  portions  of  Ag,  and  lower  velocities  are  needed  at  the  wider  portions. 
A  velocity  schedule  like  that  shown  in  Figure  1.6  is  plausible  for  an  unmanned  air 
vehicle. 


Figure  1.6  Velocity  Schedule. 
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Using  the  results  from  this  research,  they  calculate  a  number  of  schedules  for 
TPF  based  on  various  b  values.  Three  examples  for  an  autonomous  search  and 
destroy  vehicle  with  one  warhead  are  shown  in  Figure  1.7.  In  areas  where  false 
target  density  is  high,  TPF  is  low  and  relatively  flat.  In  areas  where  false  target 
density  is  low,  TPF  is  high  and  gradually  increases  as  the  time  remaining  decreases. 
The  increase  is  more  when  the  constraint  on  false  target  attacks  is  loose  and  less 
when  it  is  tight. 


Figure  1.7  Probability  of  a  Target  Report  (or  True  Positive  Fraction)  when  One 
Warhead  is  Available. 
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The  analysts  repeat  the  calculations  for  2,  3,  and  5  warheads.  The  results  for  5 
warheads  are  approximately  equal  to  the  results  for  an  infinite  number  of  warheads; 
therefore,  they  conclude  considering  more  than  5  warheads  makes  little  sense.  They 
generate  a  summary  plot  in  Figure  1.8  of  probability  of  attacking  the  target  versus 
probability  of  attacking  any  false  targets.  In  other  words,  they  generate  the  system 
operating  characteristic  broken  down  by  number  of  warheads  and  upper  bound  on 
false  target  attacks.  Each  point  on  the  graph  represents  an  optimal  schedule  for 
sensor  threshold.  Had  area  coverage  rate  been  allowed  to  vary,  each  point  would 
have  represented  the  combined  optimal  schedules  for  area  coverage  rate  and  sensor 
threshold.  Figure  1.8  shows  the  initial  desire  of  getting  the  probability  of  attacking 
the  target  to  be  greater  than  0.50  while  keeping  the  probability  of  attacking  any 
false  targets  less  than  0.30  is  not  achievable.  However,  the  Commander  now  has  a 
complete  depiction  of  the  trade  space.  He  or  she  can  look  at  one  figure  and  see  the 
effects  of  changing  constraint  level  or  changing  number  of  warheads. 


Figure  1.8  System  Operating  Characteristic. 
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1.4  Previous  Work 


The  overarching  problem  of  multi-vehicle  cooperative  control,  which  the  present 
problem  falls  under,  has  a  tremendous  amount  of  literature.  Programs  such  as 
DARPA’s  Taskable  Agent  Software  Kit  (TASK)  are  attempting  to  establish  common 
infrastructure,  mathematical  foundations,  and  metrics  for  building,  controlling,  and 
understanding  agent-based  systems.  Models  range  from  biological  [11]  to  physics- 
based  control  models  [19].  Concepts  like  reactive  behavior  [18],  game  theory  [29], 
and  decision  theory  [21],  [35]  are  being  examined.  In  addition  to  military  operations, 
applications  include  robot  navigation  [46],  highway  transportation  [6],  and  drug  in¬ 
terdiction  [9].  Each  research  area  has  its  own  wealth  of  literature.  The  present 
research  area  involves  analytic  expressions  for  effectiveness  of  autonomous  search 
and  destroy  vehicles.  The  literature  related  to  this  area  falls  into  the  following  five 
categories: 


•  Single  Target,  No  False  Targets 

•  Single  Target,  False  Targets 

•  Multiple  Targets,  No  False  Targets 

•  Multiple  Targets,  False  Targets 

•  Moving  Targets 

In  this  dissertation,  stationary  targets  and  stationary  false  targets  are  examined. 
Since  ATR  performance  is  the  focus,  false  targets  are  required.  Thus,  the  category 
of  multiple  targets,  multiple  false  targets  is  most  relevant.  For  continuity,  samples 
of  literature  found  in  each  category  as  well  as  classic  works  and  surveys  arc  listed. 

1.4.1  Classic  Works.  While  many  have  contributed  to  search  theory,  Koop- 
man,  Stone,  and  Washburn  seem  to  be  the  most  influential.  Koopman  laid  the 
foundation  in  the  1940s,  when  searching  for  U-boats  in  World  War  II  motivated  re¬ 
searchers  to  analyze  the  fundamental  search  problem.  His  work  is  summarized  in  [30] 
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and  contains  the  much  referenced  exponential  detection  law.  Stone  wrote  the  first 
unified  and  comprehensive  presentation  of  search  theory  results  in  1975  [40].  More 
recently,  he  wrote  a  text  book  on  Bayesian  Multiple  Target  Tracking  [42],  Washburn, 
currently  a  Professor  of  Operations  Research  at  the  Naval  Postgraduate  School,  has 
written  text  books  on  both  search  theory  [49]  and  game  theory  [51]. 

Two  comprehensive  review  articles  on  search  theory  exist.  Stone  provides  the 
historical  perspective  with  his  “What’s  Happened  in  Search  Theory  Since  the  1975 
Lanchester  Prize?”  article  [41],  Benkoski,  et  al.  provide  the  topic-by-topic  perspec¬ 
tive  with  their  “A  Survey  of  the  Search  Theory  Literature”  article  [4],  Together,  the 
articles  give  progress  in  search  theory  up  through  the  1980s,  where  the  stationary 
target  problem  was  thought  to  have  “reached  a  mature  state”  [45]. 

The  classic  works  of  Koopman,  Stone,  and  Washburn,  as  well  as  the  litera¬ 
ture  listed  in  the  surveys,  provide  an  analytic  foundation  for  basic  search  problems. 
However,  the  research  does  not  address  multiple-target,  multiple-false-target  prob¬ 
lems  nor  the  disposable-munition  and  multi-warhead-UCAV  problems.  In  addition, 
the  research  does  not  allow  sensor  thresholds  or  operating  conditions  to  vary  with 
time. 


1-4-2  Single  Target,  No  False  Targets.  The  single-target,  no  false-target 
problem  is  the  simplest  of  the  search  and  destroy  problems.  This  does  not  mean  the 
problem  is  trivial.  Concepts  such  as  target  distribution,  search  effort,  and  timing 
need  to  be  addressed. 

Haoi  and  Leondes  [22]  tried  to  eliminate  the  need  for  knowing  a  target’s  a 
priori  probability  density  function.  Their  goal  was  to  maximize  the  probability  of 
detecting  the  target  using  a  single-try  (non-redundant)  search  which  was  a  function 
of  search  effort  and  target  location.  They  used  Dobbie’s  [16]  extension  of  Koopman’s 
exponential  detection  law,  whereby  search  effort  was  a  function  of  the  searcher’s  loca¬ 
tion.  They  developed  a  search  plan  that  traded  off  performance  in  terms  of  detection 
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probability  for  insensitivity  to  the  target  a  priori  probability  density  function.  Their 
minimax  solution  guaranteed  a  non-zero  detection  probability,  but  not  much  more 
than  that.  This  work  is  another  example  of  how  adding  robustness  decreases  per¬ 
formance.  In  the  present  research,  search  plans  are  limited  to  linear  and  circular 
plans,  and  a  priori  probability  density  functions  are  assumed  correct.  That  does  not 
prohibit  sensitivity  studies  to  parameter  changes  in  the  probability  distributions. 

Iida,  et  al.  [23]  were  able  to  relax  the  assumption  of  local  effectiveness  of 
searching  effort,  whereby  only  objects  at  a  point  (not  objects  in  a  neighborhood)  are 
able  to  be  detected.  They  assumed  search  effort  was  dependent  on  location  and  not 
time.  Thus  image  and  data  processing  time  was  assumed  to  be  zero. 

1.4-3  Single  Target,  False  Targets.  When  two  types  of  objects  exist,  de¬ 
tection  alone  no  longer  suffices.  One  must  be  able  to  classify  an  object  as  either 
a  target  or  false  target.  Along  with  classification  comes  the  chance  to  misclassify. 
Attacking  a  false  target  or  not  attacking  a  target  can  have  grave  consequences.  For 
this  dissertation,  attacking  false  targets  is  the  focus.  Specifically,  a  maximum  is 
set  on  the  probability  of  at  least  m  false  target  attacks.  Detections  are  modelled 
based  on  target  and  false  probability  distributions.  The  detection  process  involves 
scanning  an  area  for  objects.  Once  an  object  is  detected,  a  classification  is  made. 
Correct  classification  and  misclassihcation  are  modelled  with  a  Receiver  Operating 
Characteristic  curve. 

The  approach  of  scanning  areas  for  objects,  then  classifying  all  detected  objects 
was  done  by  Stone  and  Stanshine  [43],  [44]  in  the  early  1970s.  They  used  a  broad 
search  density  function  to  calculate  search  effort.  They  used  Poisson  distributions 
to  model  false  target  encounters.  They  also  brought  in  the  idea  of  optimization. 
Specifically  they  examined  the  problem  of  minimizing  the  mean  time  to  find  the 
target.  However,  the  probability  of  classifying  a  false  target  as  a  target  was  assumed 
to  be  zero.  Either  an  object  was  correctly  classified  after  a  certain  amount  of  effort 
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or  the  searcher  moved  on.  In  addition,  their  search  plan  did  not  depend  on  the 
number  of  false  targets  found. 

Dobbie  [17]  looked  at  allowing  search  plans  to  depend  on  the  number  of  false 
targets  found.  He  allowed  for  duplicative  search,  where  the  same  spot  could  be 
searched  more  than  once.  His  results,  however,  applied  to  a  particular  case  where 
the  number  of  false  targets  was  bounded. 

1-4-4  Multiple  Targets,  No  False  Targets.  Multiple  targets  can  be  inter¬ 
preted  two  ways.  One  can  have  multiple  targets  of  the  same  type,  or  one  can  have 
multiple  types  of  targets.  In  this  dissertation,  the  focus  is  on  multiple  homogeneous 
targets.  Cozzolino  [12]  examined  multiple  targets  of  different  sizes.  He  used  Poisson 
distributions  for  target  encounters.  Since  he  used  a  sensor  craft,  with  no  limits  on 
the  number  of  target  declarations,  he  was  not  concerned  with  limiting  the  number  of 
misclassihcations.  He  simply  classified  an  encounter,  either  correctly  or  incorrectly, 
then  moved  on  with  the  search.  A  more  realistic  problem  involves  limiting  the  num¬ 
ber  of  declarations  (and  warheads  for  subsequent  attacks),  thus  conditioning  success 
on  the  number  of  previous  attacks.  Special  attention  must  be  given  to  false  target 
attack  rate  which,  if  too  high,  makes  autonomous  operation  less  attractive. 

1.4-5  Multiple  Targets,  False  Targets.  Koopman,  Stone,  Washburn,  and 
others  cited  thus  far  made  great  strides  in  developing  analytic  expressions,  but  they 
did  not  specifically  address  scenarios  involving  multiple  targets,  multiple  false  tar¬ 
gets,  and  air  vehicles  with  multiple  warheads.  Richardson  [37]  summarized  assump¬ 
tions  made  about  measures  of  effectiveness,  target  motion,  and  ways  in  which  search 
effort  is  characterized.  According  to  Richardson,  two  of  the  most  common  measures 
of  effectiveness  are  probability  of  detection  and  expected  time  to  detection.  While 
these  are  appropriate  for  certain  problems,  they  do  not  address  the  real  operational 
environment  where  false  targets  and  decoys  exist.  In  his  remarks  for  future  research, 
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Richardson  addressed  the  need  to  examine  scenarios  with  multiple  targets  and  false 
targets. 

Jacques,  Pachter,  Decker,  Jeffcoat,  and  Slater  have  all  recently  published  lit¬ 
erature  on  analytic  expressions  for  effectiveness  of  autonomous  search  and  destroy 
vehicles.  Jacques  and  Pachter  [26]  developed  analytic  expressions  for  autonomous 
wide  area  search  and  destroy  munitions  for  the  seven  scenarios  listed  in  Table  1.1. 
They  offered  a  general  approach  for  defining  task  benefits  for  cooperative  behavior. 
However,  their  expressions  assume  the  ATR  probabilities,  area  coverage  rate,  target 
parameters,  and  false  target  parameters  are  all  constant. 

In  2004,  Decker  [15]  extended  the  work  of  Jacques  and  Pachter  to  include 
multi-warhead  autonomous  search  and  destroy  vehicles.  He  proved  the  munition, 
UCAV,  and  sensor  craft  problems  could  all  be  evaluated  using  a  common  mathe¬ 
matical  framework.  However,  his  expressions  also  assume  the  ATR  probabilities, 
area  coverage  rate,  target  parameters,  and  false  target  parameters  are  all  constant. 

The  recent  expressions  derived  by  Jacques,  Pachter,  and  Decker  [26],  [15]  were 
the  starting  point  for  this  research.  Many  of  the  expressions  were  in  closed  form  and 
ideal  for  parameter  optimization  problems.  However,  the  battle  space  environment 
need  not  be  constant.  Moreover,  the  system  operating  point  in  terms  of  sensor 
threshold  and  area  coverage  rate  need  not  be  constant.  A  natural  extension  of 
their  work  involves  deriving  generalized  expressions  that  can  handle  time-varying 
parameters. 

Jeffcoat  [27]  examined  the  effect  of  cueing  on  the  probability  of  target  detection. 
He  derived  analytic  expressions  using  Markov  chains  and  established  an  upper  bound 
on  the  benefits  of  cueing  for  two  search  munitions  against  stationary  targets.  His 
results,  arrived  at  through  a  different  approach,  are  similar  to  those  of  Jacques  [25]. 
His  method  of  using  Markov  chains  was  incorporated  into  Decker’s  work  [15]. 
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Slater  [38]  looked  at  task  allocation  where  search  munitions  can  either  perform 
battle  damage  assessment  or  continue  searching  for  new  targets.  If  a  search  munition 
determines  a  target  is  still  alive,  it  mounts  a  second  attack.  Otherwise  it  continues  to 
search.  His  results  show  advantages  to  cooperative  behavior  are  heavily  dependent 
on  the  probabilities  assumed  for  the  initial  discovery  process,  the  effectiveness  of  the 
attack,  and  the  value  associated  with  striking  real  targets  versus  false  targets.  These 
results  are  in  line  with  those  produced  by  Jacques  [25].  Slater  derived  analytic 
expressions  for  two  search  munitions.  However,  his  expressions  also  assume  the 
problem  parameters  are  all  constant  throughout  the  mission. 

1-4-6  Moving  Targets.  Although  moving  targets  are  not  considered  in  this 
research,  they  exist  in  the  battle  space  and  will  eventually  need  to  be  addressed. 
Moving-target  problems  are  formidable.  Assuming  the  initial  target  distribution  is 
good  for  all  time  no  longer  works.  Data  association  and  tracking  must  be  addressed. 
Time  of  information  needs  to  be  recorded.  Models  for  target  motion  need  to  be 
derived.  The  most  common  way  to  handle  moving  targets  is  to  assume  the  motion 
is  Markov  [3],  [5],  [50].  This  makes  the  mathematics  more  tractable  and  allows  use 
of  common  techniques  such  as  Kalman  filters  [28],  [34],  Recently,  Dambreville  [14] 
departed  from  the  standard  Markov  assumption  and  considered  the  set  of  available 
target  trajectories.  He  reasoned  set  trajectories  were  more  realistic  than  simple 
probabilistic  models.  The  problem  resulted  in  games  between  the  target  and  the 
searcher.  The  games  were  solved  using  approximation  methods.  Some  proofs  were 
provided  for  special  cases  such  as  convex  evaluation  functions.  How  Dambrevillc’s 
results  scale  to  more  complex  problems  is  unclear. 

Mahler  [20],  [31],  [33]  presents  a  promising  comprehensive  approach  to  sensor 
management.  He  uses  finite  set  statistics  derived  from  point  process  theory  [13],  [47] 
and  reformulates  the  multi-sensor,  multi-target  problem  into  a  single-sensor,  single- 
target  problem.  Specifically,  all  of  the  sensors  are  bundled  into  one  giant  meta-sensor, 
and  all  of  the  targets  are  bundled  into  one  giant  meta-target.  All  individual  char- 
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acteristics  are  retained  in  the  process.  The  mechanics  to  get  solutions  involves  “set 
derivatives”  and  “set  integrals”.  Practical  solutions,  however,  involve  approximation 
techniques.  One  example  is  the  probability  hypothesis  density  filter  [32] ,  [48]  which  is 
a  1st  order  approximation  of  the  Bayes  multi-target  filter.  While  Mahler’s  approach 
seems  to  account  for  everything  (moving  targets,  false  targets,  heterogeneous  agents 
and  targets  popping  in  and  out,  etc.),  Stone  remarks  “it  is  unnecessary  and  may 
even  be  counterproductive  to  use  random  sets  to  model  the  probability  distribution 
of  all  the  ways  that  contacts  could  have  been  produced  but  were  not”  [42], 

1-4-7  System  Operating  Characteristic.  The  idea  of  establishing  a  system 
operating  characteristic  as  opposed  to  a  receiver  operating  characteristic  is  not  new. 
In  fact,  a  “receiver”  in  most  applications  is  really  a  system  of  subcomponents.  Pa¬ 
pers  by  authors  such  as  Bar-Shalom  exist  that  describe  a  transition  from  receiver 
operating  characteristic  to  system  operating  characteristic  [2] .  Regardless  of  the  ap¬ 
plication,  the  “system”  must  be  clearly  defined.  For  the  present  research,  the  system 
is  defined  as  the  area  coverage  rate  and  automatic  target  recognition  threshold.  The 
combination  of  the  sensor  and  the  vehicle,  whose  motion  together  with  the  instanta¬ 
neous  sensor  footprint,  determines  the  area  coverage  rate.  No  literature  was  found 
using  these  definitions.  Clearly  one  could  define  a  different  system  and  establish  a 
different  system  operating  characteristic. 

1.5  Research  Statement 

The  key  element  covered  in  this  research  not  covered  in  the  literature  is  time- 
varying  parameters  which  leads  to  function  optimization  as  opposed  to  parameter 
optimization.  The  result  is  a  generalized  mathematical  framework  that  handles  spa¬ 
tially  different  battle  space  environments  and  allows  system  operating  points  to  vary 
with  time.  The  mathematical  framework  applies  to  the  munition,  UCAV,  and  sensor 
craft  problems.  A  schedule  for  sensor  threshold  and  area  coverage  rate  to  maximum 
P  ( n  >  n)  for  a  given  maximum  allowable  P  (m  >  fri)  is  possible.  These  are  tactical 
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solutions  needed  by  the  hands-on  mission  planner  to  maximize  operational  effective¬ 
ness.  The  generalized  mathematical  framework  can  handle  the  following  three  types 
of  optimization  problems: 

•  Fixed  Area  Coverage  Rate,  Fixed  Threshold 

•  Fixed  Area  Coverage  Rate,  Variable  Threshold 

•  Variable  Area  Coverage  Rate,  Variable  Threshold 

Up  to  now,  only  the  first  type  was  possible. 

Expressions  needed  to  calculate  P  ( n  >  n)  and  P  (m  >  rh )  for  the  seven  sce¬ 
narios  listed  in  Table  1.1  are  derived.  To  gain  insight  into  the  problem,  analytic 
conclusions  using  the  munition  and  sensor  craft  results  for  Scenarios  1,  2,  and  7 
are  provided.  Using  these  tractable  instances,  the  problem  can  be  bounded  and  nu¬ 
merical  results  can  be  verified.  Numerical  results  from  dynamic  problems  are  also 
checked  against  analytic  results  from  fixed  problems. 

1.6  Applicability 

The  results  of  this  dissertation  can  be  used  in  the  design  phase  or  operational 
phase  of  a  system.  The  system  can  be  manned,  unmanned,  airborne-based,  land- 
based,  or  water-based.  The  sensor  can  be  RADAR,  SONAR,  IR,  LASER,  or  any 
other  type  of  sensor  whose  performance  can  be  modelled  with  a  ROC  curve.  In 
the  design  phase,  the  results  can  be  used  to  develop  system  requirements,  compare 
competing  designs,  or  verify  simulation  results.  Parameter  sensitivity  studies  can  be 
done  to  determine  the  cost  effectiveness  of  improving  sensors,  warheads,  or  vehicle 
performance.  The  mathematical  framework  enables  sound  systems  engineering.  In 
the  operational  phase,  tactical  decisions  in  the  form  of  setting  operating  points  and 
determining  the  number  of  warheads  are  possible  using  the  results  of  this  dissertation. 
The  hypothetical  operational  example  in  Section  1.3  showed  how  summary  plots  of 
P  (n  >  n)  versus  upper  bound  on  P  (rri  >  rh)  give  Commanders  a  complete  picture 
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of  the  trade  space.  He  or  she  can  look  at  one  figure  and  see  the  effects  of  changing 
constraint  level  or  changing  number  of  warheads. 

1.7  Outline  of  Document 

In  Chapter  II,  the  black  box  in  Figure  1.1  is  defined.  The  sequential  events 
method  used  to  calculate  system  effectiveness  is  explained.  Everything  needed  to 
calculate  P  (n  >  h )  and  P  (rn  >  rh )  for  the  seven  scenarios  listed  in  Table  1.1  is 
provided.  How  a  ROC  curve  relates  ATR  parameters  and  how  to  setup  an  optimal 
control  problem  are  shown.  Chapter  II,  together  with  AppendixA  and  Appendix  B, 
contain  the  new  analytic  contributions  to  the  field  where  the  fundamental  problem 
is  generalized  to  handle  time-varying  parameters. 

In  Chapter  III,  analytic  conclusions  using  the  munition  and  sensor  craft  results 
for  Scenarios  1,2,  and  7  are  provided.  The  analytic  conclusions  are  derived  using  well 
known  methods  from  optimal  control  theory.  The  analytic  conclusions  are  then  used 
to  verify  the  numerical  results  in  Chapter  IV,  where  various  types  of  optimization 
problems  listed  in  Section  1.5  are  examined.  Conclusions  and  recommendations  for 
future  research  are  given  in  Chapter  V. 
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II.  Problem  Setup  and  Definitions 

In  this  chapter,  the  black  box  in  Figure  1.1  is  defined.  The  sequential  events  method 
used  to  calculate  system  effectiveness  is  explained.  Everything  needed  to  calculate 
P  (■ n  >  n)  and  P  ( m  >  m)  for  the  seven  scenarios  listed  in  Table  1.1  is  provided.  How 
a  ROC  curve  relates  ATR  parameters  is  shown.  Finally,  how  to  setup  an  optimal 
control  problem  is  covered  along  with  a  sample  formulation. 

Thus  far,  notation  for  things  like  the  battle  space  and  the  probability  of  at 
least  n  target  attacks  have  been  defined.  To  fill  in  the  black  box  and  setup  optimal 
control  problems,  considerably  more  definitions  are  needed. 

2.1  Sequential  Events  Method 

The  purpose  of  the  sequential  events  method  is  to  produce  a  probability  density 
function  that  can  be  integrated  to  calculate  an  overall  effectiveness  or  expected  level 
of  collateral  damage.  To  model  the  search  and  destroy  mission,  the  following  area 
definitions  which  are  illustrated  in  Figure  2.1  are  needed: 

Ab  =  battle  space 

A  =  area  covered  through  time  t  or  radius  r 

As  =  area  covered  through  terminal  time  T  or  terminal  radius  R 

dA  =  area  of  sensor  footprint  at  time  t  or  radius  r 
Af  =  As-  (A  +  dA) 

Choosing  to  define  the  sensor  footprint  as  the  infinitesimal  area  dA  allows  for  a 
physical  interpretation  of  the  mathematical  operation  of  integration.  One  can  think 
of  dragging  the  sensor  footprint  over  an  area  as  an  ordered  search  effort.  For  an  event 
of  interest  (e.g.  attacking  a  target)  to  occur  in  the  sensor  footprint  at  time  t,  certain 
events  must  occur  prior  to  the  event  and  certain  events  must  occur  after  the  event. 
For  example,  if  all  warheads  were  used  prior  to  time  t,  there  is  no  way  an  attack 
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A  dA  Ay 


Figure  2.1  Area  Definitions. 

can  occur  at  time  t.  Therefore,  events  must  be  sequenced  based  on  conditions.  The 
various  events  of  interest  and  associated  numbering  scheme  is  defined  below. 

i  =  number  of  false  target  encounters  in  A 
j  =  number  of  target  encounters  in  A 

k  =  number  of  warheads  available  at  the  beginning  of  a  mission 
m  =  number  of  false  target  attacks 
M  =  number  of  false  targets  in  Ab 
n  =  number  of  target  attacks 
N  =  number  of  targets  in  Ab 
Ai  =  events  involving  false  target  enconnters 
A/"  =  events  involving  target  encounters 

T  =  events  where  targets  are  correctly  classified  and  attacked 
T  =  events  where  false  targets  are  misclassified  and  attacked 
T  =  events  where  targets  are  misclassified  and  not  attacked 
T  =  events  where  false  targets  are  correctly  classified  and  not  attacked 

For  notation,  the  number  of  events  occurring  in  a  given  area  is  subscripted 
along  with  the  applicable  area.  For  example,  n  —  1  target  attacks  in  A  is  denoted 
Tn_M.  To  calculate  the  probability  of  at  least  a  certain  number  of  attacks,  one 
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needs  to  know  the  probability  of  an  exact  number  of  attacks.  Decker  showed  one 
must  distinguish  between  the  case  where  all  warheads  are  used  and  the  case  where 
warheads  are  left  over  [15].  The  probability  of  exactly  m  false  target  attacks  and 
n  target  attacks  when  m  +  n  —  k  is  denoted  The  probability  of  exactly 

rri  false  target  attacks  and  n  target  attacks  when  m  +  n  <  k  is  denoted  Pm?nn<k- 
Decker  showed  both  cases  are  needed  to  calculate  the  overall  probability  of  an  exact 
number  of  target  and  false  target  attacks  in  As.  In  general  terms,  for  distributions 
involving  M  false  targets  and  N  targets,  he  showed 

P%T=k)(As)  =  I  [p  a  n  Tj_{n_1)tA  n  JVjv-i-j.A/  n  t14A)  x 

J  As 

P  {P m,A  D  P i—m,A  D  AA.M—i,Af)  T 

P  (%i,A  H  T j-n,A  n  Mn-3,A} )  x 

P  (Pm-1,A  H  H  H  F\,<Ia)\  (2.1) 

and 

P(rn+n<k)  (As)  =  j  [p(T„_unrHn_1MnT,_HA/nT14)x 

Jas 

P  m,A  H  J~ i—m,A  H  P M-i,Af)  + 

P  (%i,A  H  T j-n,A  n  T N-j,Af  )  X 

P  (Pm-1,A  H  Pi~(m-1),A  H  PM-l-i,Af  H  Pl,dA)\  ,  (2.2) 

where  i  is  the  number  of  false  target  encounters  in  A  and  j  is  the  number  of  target 
encounters  in  A.  The  subtle  difference  between  Eqs.  (2.1)  and  (2.2)  involves  what 
must  happen  after  time  t  or  radius  r  in  Af.  When  m  +  n  =  k,  once  the  final 
warhead  is  expended,  one  doesn’t  care  if  the  vehicle  comes  across  any  more  attack 
situations.  An  attack  situation  is  one  where  an  attack  would  occur  if  the  vehicle 
had  a  limitless  supply  of  warheads.  Therefore,  one  only  cares  about  the  number 
of  remaining  encounters  in  Af.  When  m  +  n  <  k,  there  can  be  no  attacks  in  Af. 


2-3 


Therefore,  all  subsequent  target  encounters  must  be  misclassified,  and  all  subsequent 
false  target  encounters  must  be  correctly  classified.  In  other  words,  ATs  and  A/”’s  in 
Eq.  (2.1)  become  T' s  and  T’s  respectively  in  Eq.  (2.2). 

Calculating  the  probability  of  at  least  a  certain  number  of  attacks  involves 
summing  all  the  possible  mutually  exclusive  probabilities.  Thus, 


min  (M,k—  1) 

P  (m  >  rh)  = 


min  (N,k—m—l) 


E 

71=0 


p(rn+n<k ) 
rm,n 


E 


dE-T’  (2.3) 


m= max  (rh,k—N) 


and 


min  (N,k— 1) 

P  (n  >  n)  =  ^ 


min  (M,k— n—  1) 


Ep(m+n<k) 
r m,n 


771=0 


E 


P, 


(m-\-n=k) 


k—n,m 


(2.4) 


71= max  ( h,k—M ) 


The  cases  when  either  m  =  0  or  n  =  0  must  be  considered  separately.  Either 
case  results  in  a  negative  subscript  in  Eqs.  (2.1)  and  (2.2).  If  a  subscript  is  negative, 
the  corresponding  probability  is  simply  zero.  When  all  warheads  are  used, 


p{m+n=k)  /  , 
rm.=0, n=k 


[p  (T'k-i^A  n  T j-(fc-i),A  n  MN-i-j,As  n  Ti^a)  x 

P  {Fi,A  n  MM-i,Af)\ 


(2.5) 


and 


p(m+n=k)  /  , 
Jrm=k,n=0  V^*-s 


[P  {J~k—\,A  n  Fi~[k-i),A  n  M.M-i-i,A}  n  T\ ,<m)  x 

P  (T j,A  n  AfN-j,Af)\  (2.6) 


When  warheads  are  left  over, 


-p(m+n<k) 

771=0, 71</c 


(A.)  =  P(Fo,a.) 


/  p{Tn  — 1  ,A  Cl  P' j— {n—  1),A  C  1_ jtAf  C  P\ ,c(a) 

Jas 


(2.7) 
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and 


P^n<Kn=a  (^s)  =  P  (%,As)  I  P  (Pm-1,A  H  ^-(m-  1),A  H  T M-l-i,Af  O  T\ ,<m)  (2.8) 


An  equivalent  way  to  calculate  P  (m  >  rh)  is  using 


P  (m  >  rh)  =  1  —  P  (m  <  rh) 


min  (M,m— 1) 


£ 


m= 0 


min  (N,k—m—l) 


Ep(m+n<k ) 
r m,n 


n= 0 


min  (M,m— 1) 

Ep(m+n=k) 

rm,k—m  ’  \/j"J ) 

m= max  (0,k—N) 


and  an  equivalent  way  to  calculate  P  (n  >  h)  is  using 


P  (n  >  h)  —  1  —  P  (n  <  h) 


min  (N,h—  1) 

i-  £ 


71—0 


min  (M,k—n—l) 

Ep(m+n<k ) 
r  m,n 


771—0 


min  (N,h— 1) 

E  d-diA  <2'10> 

n= max  (0  ,k—M) 


When  using  either  Eqs.  (2.9)  or  (2.10),  the  case  when  both  m  =  0  and  n  =  0  must 
also  be  considered  separately.  This  can  only  occur  when  warheads  are  left  over 
(assuming  k  >  1).  The  resulting  probability  is 


p{m+n<k)  /  . 
-rm=0,n=0 


P  P 


(2.11) 


Viewing  Eqs.  (2.1)-(2.11),  12  “elemental”  probabilities  emerge  that  need  to  be 
determined  regardless  of  the  probability  distributions  used.  Table  2.1  lists  the  12 
elemental  probabilities  and  the  corresponding  conditions. 

2.2  Elemental  Probabilities  for  Scenarios  1-7 

One  can  now  apply  Eqs.  (2.1)-(2.11)  for  Scenarios  1-7.  Details  for  calculat¬ 
ing  the  elemental  probabilities  for  each  scenario  are  given  in  Appendix  B.  Prior  to 
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Tabic  2.1  Twelve  Elemental  Probabilities 


Name 

Elemental  Probability 

Condition 

Pi 

p  (%i-i,a  n  P n  MN-\-3,Af  ^  aa) 

-p(rn-\-n=k) 

-Lm,n 

P2 

P  {P m ,A  ft  P i—m,A  ft  A ^ M —i,Af) 

-Pm,n 

P3 

P  {Pn,A  ft  P j—n,A  ft  JpN-j,Af  ) 

Y~)(m-\-n=k) 

-Pm,n 

Pa 

P  {P m— 1,A  ft  P i— (m— 1),A  ft  P^M—l—i,Af  ft  P 1,<£a) 

jj(m-\-n=k) 

-Pm,n 

P5 

P  {P i,A  ft  M.M-i,Af) 

T~>(m-\-n=k) 

m=0,n=k 

Pe 

P  {P j,A  H  JpN-j,Af  ) 

j~t(m-\-n=k) 
m=k,n= 0 

Pi 

P  (Pn-1,A  ft  P j—(n—l),A  H  P N-l-j,Af  H  Pl,dA) 

p(m+n<fc) 

-Lm,n 

PS 

P  {P in, A  ft  P i—m,A  ft  P M—i,Af) 

p(m+n<fc) 

-Lm,n 

P9 

P  (' Pn,A  ft  P j—n,A  ft  P N—j,Af) 

-p(m+n<k) 

-Pm,n 

Pi  0 

P  ( P m—l,A  ft  P i—(m—l),A  ft  P M—l—i,Af  ft  P 1  AA) 

-p(m+n>k) 

-Pm,n 

Pu 

P(F«,a.) 

p(m+n<fc) 

m=0,n<k 

P12 

p(m+n<fc) 
m<k,n= 0 

calculating  an  elemental  probability,  the  probability  of  an  exact  number  of  classi¬ 
fied  encounters  occurring  up  through  a  given  time  or  radius  for  a  given  probability 
distribution  is  needed.  These  derivations  are  given  in  Appendix  A. 

When  a  sensor  encounters  an  object,  it  compares  the  image  to  a  stored  tem¬ 
plate  or  pattern  and  either  declares  the  object  a  target  or  a  false  target.  In  practice, 
detected  objects  are  classified  to  a  certain  level  of  discrimination.  For  example,  one 
may  classify  on  object  as  either  air  breathing  or  ballistic.  A  hirer  level  of  discrimi¬ 
nation  may  be  a  specific  type  of  air-breathing  or  ballistic  object.  Regardless  of  the 
level  of  discrimination,  there  is  a  point  where  the  sensor  reports  a  detected  object  as 
either  a  target  thereby  authorizing  an  attack,  or  a  false  target  thereby  commanding 
no  attack.  Sensor  performance  is  judged  by  how  often  the  sensor  is  correct.  The 
probability  of  a  target  report  ( Ptr )  is  the  probability  the  sensor  correctly  reports  a 
target  when  a  target  is  encountered.  The  probability  of  a  false  target  report  ( Pftr ) 
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is  the  probability  the  sensor  correctly  reports  a  false  target  when  a  false  target  is 
encountered.  Together,  Ptr  and  Pftr  determine  the  entries  of  the  binary  “confu¬ 
sion  matrix”  shown  in  Table  2.2  which  can  be  used  to  determine  the  outcome  of  a 
random  draw  each  time  an  object  is  encountered  in  simulation. 


Table  2.2  Simple  Binary  Confusion  Matrix 


Declared  Object 

Encounl 

Target 

;ered  Object 
False  Target 

Target 

Ptr 

1  —  Pftr 

False  Target 

1  ~  Ptr 

Pftr 

The  expression  (1  —  Ptr)  represents  the  probability  the  sensor  reports  a  false 
target  when  a  target  is  encountered.  This  type  of  error  results  in  a  target  not  being 
attacked.  The  expression  (1  —  Pftr )  represents  the  probability  the  sensor  reports  a 
target  when  a  false  target  is  encountered.  This  type  of  error  results  in  a  false  target 
being  attacked.  For  this  binary  confusion  matrix,  true  positive  fraction  is  Ptr,  and 
false  positive  fraction  is  (1  —  Pftr)-  If  multiple  types  of  targets  are  involved,  the 
dimension  of  the  confusion  matrix  can  be  higher  [24],  This  research  only  considers 
the  2x2  confusion  matrix  in  Table  2.2. 

Allowing  Ptr  and  (1  —  Pftr)  to  vary  with  time  means  integrals  will  appear 
in  the  expressions  for  the  elemental  probabilities.  To  ease  the  appearance  of  the 
equations  and  to  help  formulate  optimal  control  problems,  states  are  defined  as  the 
integrals.  The  state  q  represents  the  probability  of  encountering  a  single  object. 
The  state  x  represents  either  the  unconditional  probability  of  attacking  a  single 
false  target  (for  Scenarios  4  and  6)  or  the  Poisson  parameter  for  false  target  attacks 
(for  Scenarios  1,  2,  3,  5,  and  7).  The  state  y  represents  either  the  unconditional 
probability  of  attacking  a  single  target  (for  Scenarios  1,  3,  4,  5,  6,  and  7)  or  the 
Poisson  parameter  for  target  attacks  (for  Scenario  2). 
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In  the  remainder  of  this  section,  a  two-part  table  is  provided  for  each  sce¬ 
nario  containing  the  12  elemental  probabilities  and  corresponding  state  definitions. 
Temporal  variables  are  used  for  Scenarios  1-4  assuming  the  mission  ends  at  time 
T.  One  can  readily  re-derive  expressions  assuming  spatial  variables;  however,  area 
coverage  rate  would  no  longer  appear  in  the  expressions.  Search  paths  for  Scenarios 
1-4  are  linear  as  depicted  in  Figure  2.2(a).  Spatial  variables  are  used  for  Scenarios 
5-7  assuming  the  mission  ends  at  radius  R.  One  can  readily  re-derive  expressions 
for  Scenarios  5-7  assuming  temporal  variables.  Area  coverage  rate  would  enter  into 
the  expressions  in  the  form  of  Q  —  2i rrr.  A  state  for  radius  r  would  need  to  be 
added.  Bounding  Q  would  result  in  state-dependent  inequality  constraints  that  are 
undefined  at  r  =  0.  For  this  reason,  spatial  variables  are  used  for  Scenarios  5-7. 
Search  paths  for  Scenarios  5-7  rely  on  concentric  annuli  as  depicted  in  Figure  2.2(b). 


(a)  Linear 


Figure  2.2  Search  Patterns 
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Table  2.3  Scenario  1  Elemental  Probabilities  and  State  Definitions 


2.2.1  Scenario  1.  For  this  scenario,  one  target  is  uniformly  distributed 
among  a  Poisson  field  of  false  targets  with  density  a  (t).  The  air  vehicle  flies  along 
a  straight  path  with  area  coverage  rate  Q  (t).  Table  2.3  lists  the  12  elemental  prob¬ 
abilities  and  the  corresponding  state  definitions. 
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Tabic  2.4  Scenario  2  Elemental  Probabilities  and  State  Definitions 


Pi 

-  e  w(t)  ^Li)!  Ptr  (t)  Q  (t)  (3  (t)  dt 

P'2 

_  pwr 

e  to! 

P'i 

_  p-i/W  b(0]n 
e  n! 

Pa 

-  e  ( }  [1  Pftr  (£)]  Q  if)  a 

( t )  dt 

P5 

—  g~x0') 

Pe 

—  e-y(t) 

Pr 

=  e  y(T)  PTR{t)Q{t)P{t)dt 

Ps 

_  p—x(T)  [x(t)]m 
c  to! 

P9 

_  e-y(T)  [2/(01n 
n\ 

Pi  0 

=  e  X(T)  1  (to— 1)!  t1  PFTR  (*)]  Q  (*)  ° 

;  (t)  dt 

Pi  1 

=  e~x{T) 

P\  2 

=  e-y(P 

X  ( t ) 

=  Jo  [!  -  Pftr  (t)]  Q  (r)  a  (r)  dr 

y(t) 

=  fo  pTR  (t)  Q  (t)  (5  (r)  dr 

2.2.2  Scenario  2.  For  this  scenario,  there  is  a  Poisson  field  of  false  targets 
with  density  a  (t)  and  a  Poisson  field  of  targets  with  density  j3{t).  The  air  vehicle  flies 
along  a  straight  path  with  area  coverage  rate  Q  (£).  Table  2.4  lists  the  12  elemental 
probabilities  and  the  corresponding  state  definitions. 
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Table  2.5  Scenario  3  Elemental  Probabilities  and  State  Definitions 


2.2.3  Scenario  3.  For  this  scenario,  N  targets  are  uniformly  distributed 
among  a  Poisson  field  of  false  targets  with  density  a  (t).  The  air  vehicle  flies  along 
a  straight  path  with  area  coverage  rate  Q  (t).  Table  2.5  lists  the  12  elemental  prob¬ 
abilities  and  the  corresponding  state  definitions.  As  stated  previously,  Scenario  1  is 
a  special  case  of  Scenario  3  with  N  —  1. 
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Tabic  2.6  Scenario  4  Elemental  Probabilities  and  State  Definitions 


Pi 

sr-^N—  1 
Z-^/j=n—  1 

{f-'K/J  bW]""1  b(t)  -  y  wr1"-11  [1  -  9  W]"-1-'}  x 

t  PTK  (t)  Q  (0  dt 

P-2 

=  E"„,{( 

:")  (;)  &  <*)]’”  [«<*)  -  *  wrm  [1  - « wr_i} 

P3 

=  Ef.„{( 

hObMl”  [?(*)  -y(i)rn 

P4 

Z—ji=m—l 

{(T1)(m‘-i)b(f)]m"1 

[q  (<)  -  X  (t)]'-*™-1'  X 

[1  ~  Pftr  (/)]  Q  (t)  dt 

P5 

=  [1  ~x(t)]M 

PG 

=  [1  -y(t)]N 

P? 

Z-^j=n—  1 

X 

7 

£ 

T 

1 

7 

e 

ho' 

, 

1 

e 

[l-9(i)-y(T)  +  !/(t)]A'-1-J' 

>  £ Ptr  (t)  Q  (t)  dt 

Ps 

= 

;")  (,!,)  [x  (*)]”  [9  «  -  1  (t)tm  [1  -  q  it)  -  x  (: T )  +  x  («)]"-} 

P9 

=  E7.„{( 

j)  (i)  b  «)]”  [9  (<)  -  V  MF”  [1  -q(t)-y  (T)  +  y  (t)]w-d 

Pio 

{r.-1)  Lb)  For1 

[9  (t)  -  x  (i)]‘-<m-1>  X 

[1  -  q  (t)  -  x  (T)  +  x  (t.)}M  1 

\  Xb  —  ^ftr  (0]  Q  (*)  dt 

Pn 

=  [i-xC. t))m 

P\  2 

=  [i  -y(T)}N 

Q(t) 

r)  dr 

x  (t) 

=  /o'  *  I1 

-  Pftr  (t)]  Q  (t)  dr 

y(t) 

=  Io^pTR{r)Q{r)dT 

2.2.4  Scenario  4 ■  For  this  scenario,  N  targets  and  M  false  targets  are 

uniformly  distributed  in  the  battle  space.  The  air  vehicle  flies  along  a  straight  path 
with  area  coverage  rate  Q  (t).  Table  2.6  lists  the  12  elemental  probabilities  and  the 
corresponding  state  definitions. 
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Tabic  2.7  Scenario  5  Elemental  Probabilities  and  State  Definitions 


Pi 

= 

EjTn-1  {  ( Y)  Li i)  {y  (0]n_1  [%  (r)  -  V  (r)]J-(n-1}  [1  -  qy 

2 

xNPtr  (r)  dr 

P2 

= 

p-x(r)[x(r)r 

C  ml 

P3 

= 

Ef.n  {(1)  Q  [y  Ml"  [<Z»  (r)  -  y  (r)rn  [1  -  q,  (r)]^} 

Pa 

= 

e  >|r|  [  (m-l)l  l1  PfT«  Ml  2,rTO  W  dr 

P5 

= 

e~x(r) 

P6 

= 

[!  -  y  (r)]N 

Pi 

— 

Ef-T,  {  (V )  („£,)  [y  (r)]”-1  \q,  (r)  -  y  (r)y  <-'>  X 
[1  -  qy  (r)  -y(R)+y  (r )]7V-1-J  |  NPTr  ( r )  ^  dr 

Ps 

= 

p-x(R)  [*(r)]m 

C  ml 

Pi 

= 

EJL „  {(!)  (i)  [y  Ml"  [q,  (r)  -  y  (r)rn  [1  -  q,  (r)  -  y  (R)  +  y  M]^'} 

Pi  0 

= 

e  i(K)[  Ol-i)!  [i  PFTR(r)]2rrra(r)dr 

Pll 

= 

e~x(R) 

P\  2 

= 

[i  -y(R)]N 

Qy(r) 

= 

fo  ^e~^dp 

x  (r) 

= 

fo  t1  “  Pftr  (p)}  Zirpa  (p)  dp 

2 

y(r) 

= 

fo  Ptr  {p)  ^2 e  ^  dp 

2.2.5  Scenario  5.  For  this  scenario,  N  targets  are  normally  distributed 
among  a  Poisson  field  of  false  targets  with  density  a  (p).  Normally- distributed  refers 
to  a  circular  normal  distribution  with  standard  deviation  ay.  The  search  area  As  is  a 
circular  disc  of  radius  R.  The  air  vehicle  searches  the  disc  using  concentric  annuli  of 
radius  r  and  thickness  dr.  The  search  begins  at  the  origin  of  the  disc  and  progresses 
outward.  Table  2.7  lists  the  12  elemental  probabilities  and  the  corresponding  state 
definitions. 
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Tabic  2.8  Scenario  6  Elemental  Probabilities  and  State  Definitions 


Pi 

=  E"=n-1  (V)  (nil)  h  (r)r1  [qy  (r)  -  y  (r)]J“^  [1  - 

_ 

xNPxR^^e  dr 

Cy 

P2 

=  E"„,  (")  (;)  [x  (>')]”  fe  (r)  -  X  (r)]J-”  [1  -  q,  (r)]"-‘ 

P3 

=  Y.U  0  (»)  Ml”  fe  M  -  v  MF”  [1  -  % 

Pa 

=  E"F,  (M“‘)  (dl)  [*  Mr 1  fe  M  -  x  (r)]-"-1'  [1 

2 

-  Qx  (r)J 

xM[l-PmW]Je  dr 

P5 

=  [1  —  x(r)]M 

Pe 

=  [l-y(rf 

P7 

=  Ehd,  f-1)  (i ,)  [y  (r)]”-1  [q,  (r)  -  y  Mp1”''1  x 

|1  -  qy  (r)  -y(R)  +  J/MF_1_j  Wf™  (r) 

u  y 

Ps 

=  E"„.  (")  (;)  [x  Ml”  fe  M  -  x  (r)r">  [1  -  q,  (r)  -  x 

(. R)+x(r)}M~i 

P9 

=  EL,  m  Q  [</  (r)]”  fe  (r)  -  y  (r)p“  [1  -  q,  (r)  -  y  (R)  +  y  (r)]K^ 

Pi  0 

=  E^m-1  (V)  L-l)  [*  (0]m_1  [&  (0  -  *  (r)]*_(m_1)  X 

2 

[1  -  Qx  ( r )  -  X  (i?)  +  X  (r)]M_1_l  M  [1  -  Pftr  (r)]  T*e 

r 

2ct2  dr 

P\\ 

=  [!-*(*)]" 

P12 

=  [i-y(*)f 

Qx(r) 

=  fo  £e~^dp  qy  (r)  =  ff  £e  dp 

2  V 

x  (r) 

=  fo  t1  -  Pftr  (p)]  fze~^dp 

2  X 

y(r) 

=  j;pTR(p)$r^dp 

u  y 

2.2.6  Scenario  6.  For  this  scenario,  N  targets  and  M  false  targets  are 
normally  distributed  (circular)  with  standard  deviations  ay  and  ax  respectively.  Ta¬ 
ble  2.8  lists  the  12  elemental  probabilities  and  the  corresponding  state  definitions. 
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2.2.1  Scenario  1.  This  scenario  is  the  specific  case  of  Scenario  5  when 
N  —  1.  Table  2.9  lists  the  12  elemental  probabilities  and  the  corresponding  state 
definitions. 
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2.3  System  Operating  Characteristic 


Everything  needed  to  calculate  P  (to  >  to)  and  P  (n>  h)  has  been  provided. 
How  the  ATR  parameters  are  related  is  now  examined.  The  parameters  Ptr  and 
P ft n  characterize  the  sensor’s  performance  and  are  not  independent.  They  are 
often  linked  by  a  Receiver  Operating  Characteristic  (ROC)  curve.  A  mathematical 
representation  of  the  ROC  curve  produces  a  graph  of  true  positive  fraction  (Ptr) 
versus  false  positive  fraction  (1  —  PFTr)  that  starts  at  (0,0),  then  monotonically 
increases  to  (1,1).  The  ROC  curve  model  adapted  from  Ref.  [36]  is  used  in  this 
reasearch: 


where  the 
algorithm, 
geometry, 


(1  —  Pftr )  — 


P 


TR 


(2.12) 


(1  —  c)  Prr  +  c 

parameter,  c  €  [l,oo),  will  depend  on  the  sensor  and  data  processing 
It  will  also  depend  on  the  vehicle  speed  (dwell  time),  and  engagement 
which  includes  flight  altitude  and  look  angle.  A  family  of  ROC  curves 


Figure  2.3  Family  of  ROC  Curves. 

parameterized  by  c  is  shown  in  Figure  2.3.  As  c  increases,  the  ROC  improves.  As 
c  — >  (X),  the  area  under  the  curve  approaches  unity  indicating  perfect  classification. 

For  a  given  sensor  and  algorithm  with  a  constant  area  coverage  rate,  the  op¬ 
erating  point  on  a  ROC  curve  is  determined  by  the  sensor’s  threshold  for  declaring 
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target.  Although  typically  it  must  be  determined  from  experimentation,  it  will  be 
assumed  the  relation  between  Ptr  and  (1  —  Pftr)  is  known,  and  can  be  approxi¬ 
mated  using  Eq.  (2.12).  Thus,  our  goal  is  finding  the  optimal  probability  of  a  target 
report  (Ptr)  -  either  constant,  or  as  a  function  of  elapsed  time  or  space  covered. 
Optimal  solutions  will  be  denoted  with  *. 

If  area  coverage  rate  is  allowed  to  vary,  the  operating  point  in  essence  moves 
from  ROC  curve  to  ROC  curve.  Better  performance  is  expected  at  slower  area 
coverage  rates,  since  longer  dwell  times  allow  for  more  sensing  and  data  processing. 
Conversely,  worse  performance  is  expected  at  higher  area  coverage  rates.  One  way  to 
model  this  effect  is  having  the  parameter  c  be  inversely  proportional  to  area  coverage 
rate  Q.  It  is  assumed  experimental  results  will  yield  a  scaled  nominal  area  coverage 
rate,  Qn,  such  that  when  Qn  is  divided  by  Q,  the  ROC  curve  is  approximated  using 
the  ROC  parameter 


in  Eq.  (2.12).  There  is  still  the  requirement  of  c  >  1,  thus  Q  <  Qn.  The  goal  is  to 
End  the  optimal  combination  of  Ptr  and  Q  -  either  constant,  or  as  a  function  of 
elapsed  time  or  space  covered.  The  combination  of  Ptr  and  Q  establishes  a  system 
operating  characteristic. 

2.4  Optimal  Control  Formulation 

An  operator  wants  the  expected  number  of  target  attacks  to  be  high  and  the 
expected  number  of  false  target  attacks  to  be  low.  One  is  interested  in  a  priori  pre¬ 
dictions  of  effectiveness  and  thus  use  P  (n  >  h)  and  P  (m  >  rh).  Unfortunately  due 
to  the  ROC,  adjusting  the  sensor  threshold  to  increase  the  number  of  target  attacks 
also  increases  the  number  of  false  target  attacks.  Thus,  the  operator’s  objectives 
are  competing,  and  a  trade-off  situation  arises.  To  ensure  P  (m  >  rh)  stays  low,  a 
constraint  is  imposed  on  it.  The  general  optimization  problem  statement  is  then 
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Max:  P  (n  >  n) 
subj  to:  P  (m  >  m)  <  b 

where  the  upper  bound  b ,  threshold  for  target  attacks  n,  threshold  for  false  target 
attacks  m,  and  mission  duration  T  are  set  by  the  designer  or  mission  planner.  The 
decision  variables  are  Ptr  and  Q. 

Everything  needed  to  calculate  P  {n  >  n)  and  P  ( m  >  fri)  has  been  provided 
as  well  as  a  way  to  relate  Ptr  and  (1  —  Pftr )•  The  problem  can  now  be  formulated 
as  an  optimal  control  problem  whose  solution  method  is  given  in  various  text  books. 
The  notation  used  by  Bryson  and  Ho  [8]  is  used  and  the  problem  is  posed  in  their 
standard  form  which  is 

Min:  J  [s  (t)  ,  u  (t)]  —  (j)  [s  (T)]  +  JQT  L  [s  ( t )  ,  u  (£)  ,  t]  dt 
Subj  to:  s  =  f  [s  ( t ) ,  u  (£) ,  t],  s  (0)  =0 

Since  their  standard  form  involves  minimizing  a  cost  functional  J,  the  objective 
function  (maximizing  P  (n  >  n))  is  multiplied  by  minus  one.  The  integrand  L  then 
becomes  the  negative  of  the  probability  density  function  for  P  (n  >  h).  The  state 
vector  s(t)  consists  of  the  states  defined  for  calculating  elemental  probabilities  (g  (f), 
x(t),  and  y  (f))  augmented  with  a  state  representing  P  (m  >  m ),  which  is  denoted 
z{t)  for  brevity.  Initial  conditions  on  the  states  are  zero.  The  terminal  conditions 
on  g,  x,  and  y  are  free,  whereas  z  (T)  <  b.  The  vector  f  [s  ( t ) ,  u  (t) ,  t]  specifies  the 
dynamics.  The  control  vector  i i(f)  consists  of  the  decision  variables  Ptr  and  Q.  The 
performance  objective  0  is  a  function  of  the  the  final  states. 

The  Pontryagin  Maximum  Principle  [8]  is  invoked,  which  involves  forming  the 
Hamiltonian,  H,  given  by 

H  [s  (: t ) ,  u  (t) ,  t]  —  L  [s  (£)  ,  u  (£),£]  +  AT  (t)  f  [s  (t) ,  u  (t) ,  t]  (2.14) 
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where  A  is  the  vector  of  costates.  Note:  T  is  mission  duration  and  T  is  the  transpose 
operator.  The  costate  differential  equations  are 


dH  dL  <9f 

ds  ds  ds 


(2.15) 


The  constraint  involving  z  amounts  to  a  “path  constraint” ,  which  results  in  A^  being 
constant  for  all  time  [8].  Since  q(T),  x(T),  and  y  (T)  are  free,  A q(T)  =  A X(T)  = 
A y  ( T )  =  0.  Finally,  the  optimal  control  is  derived  using  the  first-order  necessary 


condition 


(2.16) 


The  differential  equations  from  this  formulation  are  coupled,  so  closed-form 
solutions  do  not  exist.  Instead,  the  state  and  costate  equations  must  be  numerically 
integrated.  Two-point-boundary-value  problems  of  this  form  can  be  solved  using  a 
“shooting  method”  [7]  or  a  discrete-time,  gradient-based  solver  such  as  Matlab’s® 
fmincon  routine. 


2.5  Sample  Formulation 

To  make  sense  of  everything  covered  in  Sections  2. 1-2.4,  a  sample  formulation 
is  provided.  But  first,  here  is  a  recap  of  where  the  various  parameters  come  from. 

Designer  or  Mission  Planner:  m,  n,  b,  k,  T ,  R 

Order  of  Battle  Intelligence:  a,  f3,  N,  M,  ax ,  ay,  Ab 

Performance  Specifications:  c,  Qn ,  Qmim  Qmax 

Let  m  —  n  —  1,  and  assume  T  and  b  are  designated.  In  addition,  assume  both 
the  number  of  targets  N  and  the  number  of  false  targets  M  are  both  greater  than 
the  number  of  warheads  k.  At  this  point,  the  actual  values  for  T,  b,  N,  M,  and  k 
are  not  needed  and  can  remain  parameters.  The  optimization  problem  becomes 
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Using  Eqs.  2.10  and  2.9  gives 


Max:  P  (n  >  1) 
subj  to:  P  (m  >  1)  <  b 


and 


P  (n  >1)  =  1  —  P  (n  =  0) 

fc-i 

_  i  \  '  p(m+n<fc)  p(m+n=fc) 

^  ~  /  j  dn, 0  U fc,0 

m= 0 


(2.17) 


P  (m  >1)  =  1  —  P  (m  =  0) 

fe-i 

=  i-ipi: 


(m+n</c)  p(m+n=/c) 

n  M),fe 


(2.18) 


n= 0 


Substituting  elemental  probabilities  gives 


fc~l  rT 


P  (n  >  1)  —  1  —  P11P12  —  f  Pi 0 Pi 2  —  f  P}P( 

m=l  '/0  -70 


(2.19) 


and 


/c  1  /»T  /*i 

P(m>l)  =  l-PnPi2-  V  /  P7P11-  P1P5. 

n=l  70  70 

Bryson  and  Ho  [8]  notation  becomes 


(2.20) 


0[s  (T)]  =  PnPi2  (2.21) 

fc-i 

P  [s  (t) ,  u  (£)]  dt  =  P4P6  +  P 10P12  (2.22) 

m=l 

u(t)  =  [PT*(*)  Q(t)]T  (2.23) 

s  (*)  =  [?(*)  ®(*)  y{t)  z  (t)]T  (2.24) 
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where  q(t),  x  (t),  and  y  (t)  are  defined  once  a  scenario  is  selected,  and 


z  (t)  =  1  -  P\\P\ 


12 


(2.25) 


For  example,  if  Scenario  2  elemental  probabilities  are  used,  then 


x(t)  =  [1  -  PFTr  (t)]  Q  (r)  a  (r)  dr 

Jo 


(2.26) 


y(t)=  /  Ptr  (r)  Q  (r)  f3  (r)  dr 
Jo 

0[s(T)]  =  e-x{T)e-y{T) 


(2.27) 

(2.28) 


\r 

L  [s  (i) ,  u  (t)]  =  e~x{t)— — — —  [1  -  PFTR  (t)]  Q  (t)  a  (t)  e~y(t)+ 

[m  l): 

k— 1  r  1) 

E  <=“*(T)  T^rw  [1  -  Pftr  «)]  Q  (t)  a  (t)  (2.29) 


m=  1 


(m  —  1)! 


z(t)  =  1  -  e~x^e~y^  -  [*  <j  e~vV  [v}T)]\JpTr  (t)  Q  (t)  f3  (r)  e"*W  + 

do  [  (n-1)! 

1  r  /  \i(n— 1) 


V  e-»»|gWI,',  pTR  M  Q  M  P  M  |>  dr  (2.30) 


n= 1 


(n  —  1)! 


Everything  thus  far  was  accomplished  without  specifying  a  ROC  model.  If  constant 
area  coverage  rate  is  assumed  constant  and  Eq.  (2.12)  is  used  for  ROC  model,  then 


ct  r 


x  ( t )  = 


Ptr (t ) 


(1  -  c)  PTR  (r)  +  c 


wva  (r)  dr 


(2.31) 


y(t)  =  /  Ptr  (r)  wv/3  (r)  dr 
Jo 

</>[s  (T)]  =  e~x(T)e-y{T) 


(2.32) 

(2.33) 
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L  [s  (t) ,  u  (t)]  =  e 


—  e-x(t) 


(m-  1)!  [(!  ~  c)  pTR .(t)  +  c 


Ptr  (t) 


k- 1 

E 

771=  1 


-*(t) 


>  (0] 


(m— 1) 


-Pur  (£) 


(m  —  1)!  |_(l-c)PTR(f)+c 


wva[t)e  y^+ 

wva  (t)  e~y ^  (2.34) 


z(t)  =  1  -  e~x{t)e-yW  -  [  I  e-wwMl)L lpTR  (T)  ^  (r)  + 

io  [  (n-1)! 

fc— 1  r  /  \i(n— 1)  'j 

e~y(t)  ^wr_  1)!  ^  (r)  wvfi  (r)  e~x{t)  r  dr  (2-35) 

In  this  case,  the  only  control  variable  is  Ptr-  The  densities  a(t)  and  (5  (t)  would 
be  specified  from  the  order  of  battle  intelligence.  The  ROC  parameter  c  would  be 
taken  from  the  sensor  performance  specification  for  the  given  area  coverage  rate. 
This  completes  the  sample  formulation.  The  solution  can  be  found  using  a  shooting 
method  or  a  discrete-time,  gradient-based  solver.  The  system  effectiveness  would  be 
known  along  with  the  schedule  for  sensor  threshold.  While  jumping  straight  to  a 
numerical  solver  may  be  tempting,  insight  would  be  lost  into  the  how  the  solution 
should  behave.  Thus  in  Chapter  III,  the  optimal  control  is  analytically  calculated 
using  =  0  for  tractable  instances  of  Scenarios  1,  2,  and  7. 
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III.  Analysis 

In  this  chapter,  analytic  conclusions  are  provided  for  the  munition  and  sensor  craft 
problems  for  Scenarios  1,  2,  and  7.  Using  these  tractable  instances,  the  problem  can 
be  bounded  and  insight  can  be  gained  prior  to  calculating  numerical  solutions.  Ten 
theorems  are  presented  from  which  the  following  conclusions  are  drawn: 

1.  For  Scenarios  1  and  2  with  k  =  1  and  no  constraint  on  P  ( m  >  m),  if  T,  Q, 
a,  and  /3  are  constant,  the  unbounded  PfR  is  monotonically  increasing  with  time. 

2.  For  Scenarios  1,  2,  and  7  with  k  =  oo,  n  —  1,  and  no  constraint  on 
P  [in  >  m),  if  T  (or  R ),  Q,  ct,  and  (5  are  constant,  the  bounded  PIR  =  Pin,,,,,,  =  1- 

3.  For  Scenarios  1  and  2  with  k  =  oo  and  P  [yrt  >  m)  =  b,  if  T,  Q,  a,  and  (3 
are  constant,  then  P£R  is  constant  for  all  time. 

4.  For  Scenarios  1  and  2  with  k  =  1  or  k  =  oo,  if  T,  a,  and  [3  are  constant  and 
Q  ( t )  is  unbounded,  then  the  unbounded  PIR  is  given  by 


Ptr  (t)  ~ 


Qr 


2Qn  ~  Q*  (' t ) 


(3.1) 


5.  For  Scenario  7  with  k  =  oo  and  P  (m  >  m)  =  b,  if  R  and  a  are  constant, 
then  PIR  is  monotonically  decreasing  with  radius. 

The  remainder  of  this  chapter  contains  detailed  calculations  required  for  the  10 
theorems.  The  general  approach  involves  forming  a  Hamiltonian  for  each  problem, 
then  solving  ^  =  0  to  determine  the  form  of  the  optimal  control.  Next  the  derivative 
of  the  optimal  control  with  respect  to  time  or  radius  is  taken  to  determine  any  trends 
(e.g.  increasing,  decreasing,  or  constant).  The  reader  can  skip  to  Chapter  IV  for 
applications  that  demonstrate  the  conclusions  above. 
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3.1  Scenario  1:  Fixed  T,  Fixed  Q,  Variable  Ptr 

For  Scenario  1,  one  target  is  uniformly  distributed  among  a  Poisson  field  of 
false  targets  with  density  a,  which  is  assumed  constant.  The  air  vehicle  flies  along 
a  straight  path  with  area  coverage  rate  Q  ( t ).  For  now  Q  is  constant  and  specified. 
Therefore,  the  only  decision  variable  is  Ptr ■  Constant  Q  implies  the  ROC  parameter 
c  is  constant. 

3.1.1  Munition  Problem.  To  start,  in  addition  to  assuming  T,  Q,  and  a 
are  constant,  there  is  no  constraint  on  P  ( m  >  m). 

Theorem  1  (Scenario  1  with  k  =  1  and  no  constraint  on  P  (m  >  m ))  IfT,  Q, 

and  a  are  constant,  the  unbounded  PRR  is  monotonically  increasing. 

Proof:  Let  the  state  x  be 

x{t)=  [  aQ  [1  -  PFTr  (t)]  dr.  (3.2) 

Jo 

The  corresponding  L  is 

L  =  -^-PTR(t)e~xW.  (3.3) 

A-b 

The  derivative  of  the  state  is 


x  (t)  =  aQ[  1  -  Pftr  (i)] 


(3.4) 


After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 


H  =  - S-PTRe -  +  KaQ- - 

Ab  (1  -  c)  PTR  +  c 


(3.5) 


and  the  costate  differential  equation  becomes 


Ac  (t)  —  —-j-Ptr  (t)  e 
Ab 


(3.6) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  _  Q^x  XxaQc 

3j Ptr  A-b  r/1  °  1  ~l2 


[(1  —  c)  Ptr  +  cY 


Solving  -7^^  =  0,  the  optimal  control  is 


Ptr  (t)  — 


c  ±  a/ Xx  (t)  acABex(P> 


c  —  1 


(3.7) 


(3.8) 


Only  the  “minus”  root  is  of  interest,  since  the  “pins”  root  puts  Ptr  outside  of  [0, 1] 
for  all  time.  Taking  the  derivative  of  Eq.  (3.8)  with  respect  to  t  gives 


Ptr  (t) 


—  1  /  cxcAe 

2  (c  -  1)  y  \x  (■ t )  ex<t) 


\x  (t)  x  (■ t )  +  Xx  (■ t ) 


exit). 


(3.9) 


Using  Eqs.  (3.4)  and  (3.6)  for  x  (t)  and  Aa,  (t)  yields 


Ptr  (t)  -  Q\j  icA~XXx(t)  [Ptr  Wf  >  0-  (3.10) 

□ 

At  the  terminal  point  where  Xx  (T)  =  0,  P^R  (T)  =  ^  >  1.  Thus,  the 
constraint  Ptr  <  1  becomes  active  prior  to  the  end  and  remains  active  until  time 
T.  See  Appendix  D  for  further  analysis  on  when  the  constraint  becomes  active. 
This  represents  a  “go  for  broke”  tactic  in  the  end  game.  This  strategy  of  starting 
conservative  then  gradually  increasing  aggressiveness  until  the  end,  where  an  all  out 
effort  is  tried,  is  common  in  game  theory.  Without  a  constraint  on  false  target 
attacks,  the  system  has  nothing  to  lose  by  declaring  all  objects  as  targets  near  the 
end. 
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Next,  the  effects  of  adding  a  constraint  on  false  target  attacks  is 
With  a  constraint  on  P  ( m  >  rh ),  the  following  states  are  needed: 


x  (t)  = 

y(t)  = 

z{t)  = 


aQ  [1  -  PFTr  (r)]  dr 


>o 

[  -j-Ptr  (t)  dr 
I o  AB 


ft 


a 


Q[  1  -  Pftr  (t)]  [1  -y  (r)]  e  *(T)dr. 


The  corresponding  L  is  still 


L  =  -~j-PTr  (t)  e~*V. 

FiR 


The  derivatives  of  the  states  are 


x  (t)  —  aQ  [1  —  Pftr  (£)] 
y(t)  =  ~t~Ptr  (t) 

ar 

z  ( t )  =  aQ  [1  -  PFTr  (t)]  [1  -y  (t)]  e~x{t) 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

H  =  (x«- «-)  t/™ + ^(1-cml+c  A  - + m d  ■ -  y)  «-l 

and  the  costate  differential  equations  become 


Ax  (t) 

A y  (t) 
A *  (0 


A  z^Q 

\Z®Q 


Ptr  (t) _ 

(1  —  c)  Ptr  ( t )  +  c 

Ptr  (t) _ 

(1  —  c)  Ptr  (0  +  c 


[i  -y(t)} 

e~xd) 


e~xd)  _ 


-j-Ptr  (t)  e 


0 


examined. 

(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

~x(P  (3.19) 

(3.20) 

(3.21) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  _X^Q  [Ax  +  \z(l-y)e  x]aQc 

9PTr  ~  v  6  ’  Ab+  [(1  -c)PTi?  +  c]2 

Solving  =  0,  the  optimal  control  is 


Ptr  (J) 


c  ± 


{Aa:(t)+Az[l-3/(t)]e  “Wj-a&Afl 


e~<*) -\y(t) 


C  —  1 


(3.22) 


(3.23) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1] 
for  all  time.  While  taking  the  derivative  of  Eq.  (3.23)  with  respect  to  t  is  possible, 
there  are  too  many  terms  involved  to  discern  any  trends.  Nonetheless,  Eq.  (3.23) 
shows  the  form  of  the  optimal  control. 


3.1.2  Sensor  Craft  Problem.  The  sensor  craft  problem  does  not  condition 
P  (■ n  >  n)  on  the  number  of  previous  attacks.  For  Scenario  1, 

P(n  =  1)=  [  -t~PTr  (t)  dt.  (3.24) 

Jo  ab 

Theorem  2  (Scenario  1  with  k  =  oo  and  no  constraint  on  P  (to  >  m))  If  T , 

Q,  and  a  are  constant,  the  optimal  constrained  solution  is  PfR  =  PfRmax  —  1- 

Proof:  With  no  constraint  on  P  (rn  >  to),  the  optimal  control  is  clearly  PfR  = 
PfR  =  1  for  all  time. 

l  Tlmax 


□ 

This  is  not  a  practical  solution,  since  all  false  targets  would  be  declared  targets  and 
attacked.  With  an  infinite  amount  of  warheads  and  no  constraint  on  false  target 
attacks,  there  is  no  reason  not  to  attack  every  object  encountered. 
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Next,  the  effect  of  adding  a  constraint  on  false  target  attacks  is  examined. 
The  sensor  craft  problem  does  not  condition  P  ( m  >  rh)  on  the  number  of  previous 
attacks. 


Theorem  3  (Scenario  1  with  k  =  oo,  rh  =  1,  and  P  (rn  >  1 )  =  b)  If  T,  Q,  and 

a  are  constant ,  PfR(t)  =  (1_e)Xp^+aQT- 

Proof:  With  a  constraint  on  P  (rn  >  rh),  the  states  become 


x  (t)  —  /  aQ  [1  —  Pftr  (t)]  dr 


rh—  1 


w  =  i-EMTe-4. 


m= 0 


m\ 


(3.25) 

(3.26) 


The  corresponding  L  is 


L  = 


-j-Ptr  (t) 
A-B 


The  derivatives  of  the  states  are 


(3.27) 


x  (t)  =  aQ  [1  -  PFTr  (t)] 
z  (t)  =  aQ  [1  —  Pftr  (£)] 


>  (*)] 


(rh— 1) 


(rn  —  1)! 


D-x(t) 


(3.28) 

(3.29) 


After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 


u  Q  d  i  m  Ptr 

H  -  — t~Ptr  +  oiQ- - — - - — 

Ab  (1  —  c)  Ptr  +  c 


A x  +  A2 


x 


(rh—  1) 


(m 


1)! 


(3.30) 


and  the  costate  differential  equations  become 


A x  (t)  XzotQ 

A  z(t)  =  0 


Ptr  ( t )  (rh  -  1) 

(1  -  c)  PTr  (t)  +  c  x  ( t ) 


x  (<)] 


(rh— 1) 


(rh  —  1)! 


e-xW(3.31) 


(3.32) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


BH_  =  0 

dPTR  AB+  [(1  -c)PTR  +  cf  K'  ] 

Solving  jAP  =  0,  the  optimal  control  is 

c  ±  J |  Xx  ( t )  +  A*  }  e~x(A  |  acAB 

Ptr  it)  = - 1 - " - •  (3-34) 

Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Prr  outside  of  [0, 1] 
for  all  time.  Taking  the  derivative  of  Eq.  (3.34)  with  respect  to  t  gives 


Using  Eq.  (3.28)  for  x  (t)  with  Eq.  (2.12)  for  the  ROC  model  yields 

Hn(t)  =  0  (3.36) 


implying  P^R  is  constant  for  all  time.  When  target  and  false  target  densities  remain 
fixed  throughout  the  search  effort  and  the  amount  of  warheads  is  infinite,  the  problem 
becomes  an  infinite  horizon  problem,  where  optimal  solutions  are  constant.  Consider 
the  case  when  m  —  1.  For  constant  T,  Q,  a,  and  Ptr, 

P  (m  >  1)  =  1  -  e~aQT{1~PFTR)  <  b.  (3.37) 
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Applying  Eq.  (2.12)  for  the  ROC  model  and  assuming  the  constraint  is  active,  the 
optimal  control  is 


Ptr  {t) 


— cln  [1  —  b] 

(1  —  c)  In  [1  —  b]  +  aQT 


(3.38) 

□ 


3.2  Scenario  1:  Fixed  T,  Variable  Q,  Variable  Ptr 

The  assumption  that  Q  is  constant  is  now  relaxed.  Therefore,  the  ROC  para¬ 
meter  c  is  no  longer  constant.  Equation  (2.13),  where  c(t)  =  is  used.  Recall  Qn 
is  a  constant.  The  decision  variables  are  Q  and  Ptr- 

Theorem  4  (Scenario  1  with  k  =  1  or  k  =  oo)  If  T  and  a  are  constant  and  Q 
is  unbounded,  the  unbounded  PfR  (t)  =  • 

Proof:  Consider  both  the  munition  and  sensor  craft  problems. 

3.2.1  Munition  Problem.  For  k  —  1,  dynamic  Q,  and  dynamic  Ptr,  the 
states  becomes 

x(t)  =  [  aQ  (r)  [1  -  PFTr  (t)]  dr  (3.39) 

Jo 

y(t)  =  [  -t-Q  (t)  Ptr  (t)  dr  (3.40) 

Jo  ab 

z(t)  =  [  a<5(r)[l-  PFTr(t)][1-  y  (r)}e~x{-T)  dr.  (3.41) 

Jo 

The  corresponding  L  is 

L=-^-Q{t)PTR{t)e~<t\  (3.42) 
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The  derivatives  of  the  states  are 


x  (t)  =  aQ  (t)  [1  -  Pftr  (t)] 

y(t)  =  4 -Q(t)PrR(t) 

A-b 

z  (t)  =  aQ  (: t )  [1  -  PFTR  (t)]  [1  -y  (t)]  e~x{t) 


Equation  (2.12)  becomes 


[1  —  Pftr  (£)] 


_ Q  (t)  Ptr  it) _ 

Q  (t)  Ptr  (t)  +  Qn  [1  —  Ptr  (£)] 


(3.43) 

(3.44) 

(3.45) 


(3.46) 


Using  Eq.  (3.46),  the  Hamiltonian  becomes 

H =  ( V ~  e~)  +  aQPTjQPXl-PTR)  +  (1  -  e "I  (3'47) 


and  the  costate  differential  equations  become 


A 


y 


K 


,  „  Q2  (t)  Ptr  (t) 

a  rv - 

Q  (t)  Ptr  ( t )  +  Qn  [1  —  Ptr  (£)] 
A-B 


[1  -y(t)]e-x® 


A  7a 


Q2  (t)  Ptr  (t) 


Q  (t)  Ptr  ( t )  +  Qn  [1  —  Ptr  (t)] 


D-x(t) 


0 


(3.48) 

(3.49) 

(3.50) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  I  r\  -aq  1  [A*  +  K  (1  -  y)  e  x]  QnaQ 

QP  ^  1  VA:y  e  J  + 


TR 


Ab  [QPtR  +  Qn  (1  —  Ptr)Y 


(3.51) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Q  gives 


dH  n  f,x  1  ,  rx  ,  x  ,  _xl  Q2PTR  +  2QnQ(l-PTR) 

SQ-PTRVK-e  )^  +  «U  +  V(l-!,)e  1  {QPTR  +  Qn{1-PTR)f 


(3.52) 
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Ignoring  the  trivial  solutions  Ptr  —  0  and  Q  =  0,  the  necessary  conditions  become 


1  [Ag  +  Ag  (1  -y)e  x }  QnaQ  1  _  Q 
^4_b  [ QPtr  +  Qn  (1  —  Ptr)}2  J 


(3.53) 


and 


1 

Ar 


+  oc  [Ax  +  Az  (1  —  y)  e  x 


Q2 Ptr  +  2 QnQ  (1  —  Ptr)  ) 
[QPtr  +  Qn  (1  —  -Ptr)]2  J 


0  (3.54) 


Subtracting  Eq.  (3.53)  from  Eq.  (3.54)  gives 


a  [Ax  +  Az  (1  -y)e  x 


Q2Ptr  +  2QnQ  (1  —  Ptr) 
[QPtr  +  Qn  (1  —  Ptr)}2 


[Ag  +  Az  (1  -y)e  x}  QnaQ 
[QPtr  +  Qn  (1  _  Ptr)}2 

=  0  (3.55) 


Grouping  like  terms  yields 


aQ  \x  +  Az  (1  —  y)  e  x 


QPtr  ~  ZQuPtr.  +  Qn 
QPtr  +  Qn(  1  —  Ptr)}2 


0 


resulting  in  the  necessary  condition 


(3.56) 


[Ag  +  Az  (1  -  y )  e~x]  [QPtr  -  2 QnPTR  +  Qn\  =  0.  (3.57) 


Thus,  either 


Ptr  (t) 


Qn 

2Qn  -  Q*  (t) 


(3.58) 


or 

\x(t)  =  -\z[l-y(t)}e-xW.  (3.59) 

In  the  case  where  PRR  (t)  =  insight  can  be  gained  looking  to  the  relations 

between  the  parameters.  From  the  requirement  that  c  >  1,  Q*  (t)  <  Qn.  So  as  Q* 
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increases,  PfR  increases.  Further,  for  Q*  e  (0,  Qn),  PfR  £  (|,  l).  In  other  words, 


Ptr  (t)  >  0-5 


(3.60) 


For  Qn  >>  Q*  (which  is  the  case  when  c  =  100),  PfR  ~  0.5  for  all  time.  In  addition, 
the  relation  between  PfR  and  Q*  does  not  depend  on  a  or  the  value  of  b. 

The  case  where  Xx  (t)  =  —  A2  [1  —  y  (£)]  e~x ®  never  occurred  in  practice.  In  ad¬ 
dition,  the  case  where  PfR  (t)  =  ,?Q  sh°ws  up  again  and  again  in  this  Chapter 

and  in  practice.  Therefore,  it  will  be  the  focus. 


3.2.2  Sensor  Craft  Problem.  For  k  —  oo,  dynamic  Q,  and  dynamic  Ptr, 
the  states  become 


x  (t)  =  /  aQ  (r)  [1  -  PFTr  (r)]  dr 


rh—1 


(t)  =  i-E 


m= 0 


WTC-X(t) 

rri\ 


(3.61) 

(3.62) 


The  corresponding  L  is 


L  —  ~-{-Q  (t)  Ptr  (t)  ■ 

A-  B 


The  derivatives  of  the  states  are 


(3.63) 


x  (■ t )  =  aQ  (t)  [1  -  Pftr  (t)] 

Ht)  =  aQ  (t)  [i  -  m 


(3.64) 

(3.65) 


After  applying  Eq.  (3.46),  the  Hamiltonian  becomes 


t t_  j_np  _ Q2  Ptr _ 

AbQ  TR  QPTr  +  Qu(1-Ptr) 


Xx  +  Xz 


x  (t)] 


(m— 1) 


(m  —  1)! 


(3.66) 
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and  the  costate  differential  equations  become 


-\r  (f)  XzOi 


Q2  (t)  Fruit) 


Q  (t)  PtR  (' t )  +  Qn  [1  -  PtR  CO]  L  X  it) 


(to  —  1) 


X  it)] 


(m—  1) ' 


( to  —  1)! 


~^~Q  it)  PTR  it.) 
A-b 


Kit)  =  o 


z~x{t) 

(3.67) 

(3.68) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  ) 

9Ptr  |  Ar 


K  +  K 


\x(t) 

(m— 1)! 


Qn^Q 


[QPtR  +  Qn  (1  —  Ptr)Y 


(3.69) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Q  gives 


dH  n  1 

dQ-PTK\-A-B+a 


K  +  A2 


>  it)] 


(m— 1) ' 


Q2Ptr  +  2 QnQ  (1  —  Ptr) 


(to  —  1)! 

Ignoring  the  trivial  solutions  Ptr  =  0  and  Q  =  0,  the  necessary  conditions  become 


[QPtr  +  Qn  (1  —  Ptr)]2 

(3.70) 


— ; - l~ 


AI  +  A,ftE@1 


i(m-l) 


(to— 1)! 


QnOtQ 


Ab  [QPtr  +  Qn  (1  —  Ptr)Y 


=  0 


(3.71) 


and 


■  — - b  01 

Ab 


Xx  +  A2 


x  (t)] 


(m— 1) ' 


(to  —  1)! 


Q2Ptr  +  2QnQ  (1  ~  Ptr) 
[QPtr  +  Qn  (1  —  Ptr)]2 


=  0  (3.72) 


Subtracting  Eq.  (3.71)  from  Eq.  (3.72)  gives 


a 


K  +  A2 


>  it)] 


(m— 1) 


(to  —  1)! 


Q2Ptr  +  2 QnQ  (1  —  Ptr) 
[QPtr  +  Qn  (1  —  Ptr)]2 


Xx  +  A2 


(m— 1)! 


[QPtr  +  Qn  (1  —  Ptr)Y 


=  0  (3.73) 
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Grouping  like  terms  yields 


aQ 


\x  +  Xz 


>  (t)] 


(m— 1) ' 


(to  —  1)! 


QPtr  —  2 QuPtr  +  Qn 
QPtR  +  Qn  (1  ~  Ptr)X 


=  0 


(3.74) 


resulting  in  the  necessary  condition 


Xx  +  A* 


x  (t)] 


(m—  1) ' 


(m  —  1)! 


[QPtr  ~  ZQnPTR  +  Qn]  —  0.  (3.75) 


Thus,  either 


or 


Ptr  (t)  ~ 


Xx  (t)  Xz 


Qr 


2Qn  -  Q*  (' t ) 


%  (*)} 


[fh—  1) ' 


D-x(t) 


(■ m  —  1)! 

Equation  (3.76)  is  the  same  necessary  condition  as  the  single  warhead  case. 


(3.76) 


(3.77) 


□ 

Thus,  it  appears  the  necessary  condition  in  Eq.  (3.76)  holds  regardless  of  the  number 
of  warheads,  false  target  density,  or  maximum  allowable  P  (to  >  to).  This,  of  course, 
assumes  unbounded  design  variables.  If  Q  saturates  at  an  upper  or  lower  bound, 
the  fixed-Q  necessary  conditions  are  used.  Numerical  examples  in  Chapter  IV  will 
address  the  bounded-Q  case. 


3-13 


3.3  Scenario  2:  Fixed  T,  Fixed  Q,  Variable  Ptr 

For  Scenario  2,  there  is  a  Poisson  field  of  targets  with  density  /?  and  a  Poisson 
field  of  false  targets  with  density  a.  Both  f3  and  a  are  assumed  constant.  The  air 
vehicle  flies  along  a  straight  path  with  area  coverage  rate  Q  (t).  For  now  Q  is  assumed 
constant  and  specified.  Therefore,  the  only  decision  variable  is  Ptr ■  Constant  Q 
implies  the  ROC  parameter  c  is  constant. 

3.3.1  Munition  Problem.  To  start,  in  addition  to  assuming  T,  Q,  (3,  and 
a  are  constant,  there  is  no  constraint  on  P  (m  >  m). 

Theorem  5  (Scenario  2  with  k  —  1  and  no  constraint  on  P  (■ m  >  m))  IfT,  Q, 

(3,  and  a  are  constant,  the  unbounded  PRR  is  monotonically  increasing. 

Proof:  The  only  state  needed  is 

x(t)  =  [  {[1-  Pftr(t)}(x  +  PTr(t)  P}Qdr  (3.78) 

Jo 

which  is  the  combined  Poisson  parameter  for  both  target  attacks  and  false  target 
attacks.  The  corresponding  L  is 

L  =  -PTR(t)  PQe~x®.  (3.79) 

The  derivative  of  the  state  is 

x  (t)  =  {[1  —  Pftr  (£)]  ot  +  Ptr  if)  P}  Q ■  (3.80) 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

H  =  -PTRPQe~x  +  [  P™a  +  Ptr(3\  Q  (3.81) 

L(1  -c)PTR  +  C  J 
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and  the  costate  differential  equation  becomes 


K  (t)  =  -PTR(t)PQe-xV. 


(3.82) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


d  H 


dP, 


—  — (3Qe  x  +  A XQ 


TR 


ac 


[(1  -  c)  PTr  +  c]‘ 


+  P 


(3.83) 


Solving  =  0,  the  optimal  control  is 


c  ± 


^tr  (t)  — 


ac\x(t) 


p\e-x(t)_Xx(t) 


C  —  1 


(3.84) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1] 
for  all  time.  Taking  the  derivative  of  Eq.  (3.84)  with  respect  to  t  gives 


Ptr  (t)  ~ 

Using  Eqs.  (3.80)  and  (3.82)  for  x  (t)  and  (t)  yields 
Ptr  (t)  —  Q 


-1 

/  ac 

e-xit) -\x(t)] 

'  e-x(t) 

Az  (t)  x  (■ t )  +  Xx  ( t ) 

\ 

2  (c  -  1)  V 

f3Xx  (■ t ) 

[e-AO  -  A,  (f)]2  J 

(3.85) 


a/3e  2a:d) 


4cAx  ( t )  [e  x(*)  -  A*  (t)] 


[Ptr  W]  >  0- 


(3.86) 


□ 

At  the  terminal  point  where  Ax  (T)  =  0,  P^R  ( T )  =  >  1.  Thus,  the  constraint 

Ptr  A  1  becomes  active  prior  to  the  end  and  remains  active  until  time  T.  This 
represents  the  same  “go  for  broke”  tactic  seen  in  Scenario  1. 
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Next,  the  effect  of  adding  a  constraint  on  false  target  attacks  is  examined. 
With  a  constraint  on  P  ( m  >  m),  the  following  states  are  needed: 

x  (t)  =  [  { [1  -  PFTr  (t)]  a  +  Ptr  (j)  (3}  Qdr 

Jo 

z{t)  =  [  aQ  [1  —  Pftr  {t )]  e~x^dr. 

Jo 

The  corresponding  L  is  still 

L  =  —PTR  (t)  (3Qe~x^ .  (3.89) 

The  derivatives  of  the  states  are 


(3.87) 

(3.88) 


x  (t)  —  {[1  —  Pftr  +  Ptr  (t)  fl}  Q  (3.90) 

z  (■ t )  =  aQ  [1  —  Pftr  (£)]  e  (3.91) 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

H  =  (\x  —  e~x)  Ptr(JQ  +  Ptr*Q  [\x  +  Xze~x]  (3.92) 

(1  -  c)  PTR  +  c 

and  the  costate  differential  equations  become 

A*  (t)  =  A zaQ- - e~*V  -  PTR  (t)  (JQe~^  (3.93) 

(1  -  c)  PTr  (t)  +  c 

A  z{t)  =  0  (3.94) 


Taking  the  partial  derivative  of  PI  with  respect  to  the  decision  variable  Ptr  gives 


dH 

<JPtr 


(Ax  -  e~x)  (JQ  + 


[A*  +  A2e  x]  aQc 
[(1  —  c)  Ptr  +  c]2 


(3.95) 
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Solving 


0,  the  optimal  control  is 


Ptr  (0 


c  ± 


ac\\x(t)+\ze  *(*) 


me  xW-\x(t) 


C~  1 


(3.96) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1] 
for  all  time.  While  taking  the  derivative  of  Eq.  (3.96)  with  respect  to  t  is  possible, 
there  are  too  many  terms  involved  to  discern  any  trends.  Nonetheless,  Eq.  (3.96) 
shows  the  form  of  the  optimal  control. 


3.3.2  Sensor  Craft  Problem.  The  sensor  craft  problem  does  not  condition 
P  (■ n  >  n)  on  the  number  of  previous  attacks.  For  Scenario  2, 


P  (■ n  >  h) 


h— 1 

£ 

n= 0 


Jo  Ptr  ( t )  (3Qdt 


foT  PmWPQdt 


n\ 


(3.97) 


Theorem  6  (Scenario  2  with  k  =  oo,  n  =  1,  and  no  constraint  on  P  {m  >  m)) 

If  T ,  Q,  and  f3  are  constant,  the  optimal  constrained  solution  is  PfR  =  PRR  =  1. 

Proof:  In  the  case  of  n  —  1, 


P(n>  1)  =  1  -  e-fo  pTRd)PQdt^ 


(3.98) 


and  with  no  constraint  on  P  (rn  >  rh),  the  optimal  control  is  clearly  PfR  =  PfRq 
1  for  all  time. 


□ 

This  is  not  a  practical  solution,  since  all  false  targets  would  be  declared  targets  and 
attacked.  With  an  infinite  amount  of  warheads  and  no  constraint  on  false  target 
attacks,  there  is  no  reason  not  to  attack  every  object  encountered. 
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Next,  the  effect  of  adding  a  constraint  on  false  target  attacks  is  examined. 
The  sensor  craft  problem  does  not  condition  P  ( m  >  m)  on  the  number  of  previous 
attacks. 


Theorem  7  (Scenario  2  with  k  =  oo,  m  =  1,  and  P  (rn  >  1 )  =  b)  If  T,  Q,  and 

a  are  constant ,  PfR(t)  =  (1_e)Xp^+aQT- 

Proof:  With  a  constraint  on  P  (rn  >  m),  the  states  become 


x(t) 

y(t) 

z(t) 


aQ  [1  -  PFTr  (t)]  dr 


(3QPtr  (r)  dr 


rh—  1 


m= 0 


MTe-*w 

ml 


In  this  formulation,  <f>  is  used  and  L  —  0.  The  corresponding  <j)  is 

f  h—l 

*=-  i-E- 


n= 0 


n\ 


The  derivatives  of  the  states  are 


(3.99) 

(3.100) 

(3.101) 


(3.102) 


X  (t) 

—  aQ  [1  —  Pftr  (i)] 

(3.103) 

y(t) 

=  / SQPtr  (' t ) 

(3.104) 

z(t) 

■  . 

(3.105) 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 


H  =  aQ 


P 


TR 


(1  -  c)  PTR  +  C 


A-r  +  A2 


x 


(m-l) 


L(m-  !)!J 


e~x  }  +  A vPQPTr  (3.106) 
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and  the  costate  differential  equations  become 


K  (t)  = 

A  y(t)  = 
K  (t)  = 


-A zaQPTR{t)  f(m-l)  1  J  [; 


(1  -  c)  PTr  it)  +  c  |_  x  (t) 


( rh  —  1)! 


(3.107) 

(3.108) 

(3.109) 


Xv  ~ 


r  \—  _il  MZT  -W(T) 

~^[y(T)  \  n\ 


(3.110) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  aQe{A,  +  v[ggg], 

9Ptr  [(1  -c)PTR  +  c]2 

Solving  jvtt—  =  0,  the  optimal  control  is 


+  A  y/3Q 


(3.111) 


Ptr  (t)  ~ 


(t)  +  Az  ■ 


(m— 1)!  J  p\y 


C  —  1 


(3.112) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1] 
for  all  time.  Taking  the  derivative  of  Eq.  (3.112)  with  respect  to  t  gives 


Ptr  (t)  ~ 


Using  Eq.  (3.103)  for  x(t)  with  Eq.  (2.12)  for  the  ROC  model  and  Eq.  (3.107)  for 
A*  (t)  yields 

Ptr{P)  ~  0  (3.114) 
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implying  PRR  is  constant  for  all  time.  When  target  and  false  target  densities  remain 
fixed  throughout  the  search  effort  and  the  amount  of  warheads  is  infinite,  the  problem 
becomes  an  infinite  horizon  problem,  where  optimal  solutions  are  constant.  Consider 
the  case  when  m  —  1.  For  constant  T,  Q,  a ,  and  Ptr, 

P(m>  1)  =  1  -  e-aQT{1~PFTR)  <  b.  (3.115) 


Applying  Eq.  (2.12)  for  the  ROC  model  and  assuming  the  constraint  is  active,  the 
optimal  control  is 


Ptr  (t) 


— cln  [1  —  b\ 

(1  —  c)  In  [1  —  b]  +  olQT 


(3.116) 

□ 


Both  Scenario  1  and  Scenario  2  have  the  same  Poisson  held  of  false  targets  with 
density  a.  With  infinite  warheads,  the  constraint  on  false  target  attacks  will  produce 
the  same  optimal  solution  regardless  of  the  target  distribution.  Thus,  Eq.  (3.116) 
matches  Eq.  (3.38). 


3-4  Scenario  2:  Fixed  T,  Variable  Q,  Variable  Ptr 

The  assumption  that  Q  is  constant  is  now  relaxed.  Therefore,  the  ROC  para¬ 
meter  c  is  no  longer  constant.  Equation  (2.13),  where  c(t)  =  is  used.  Recall  Qn 
is  a  constant.  The  decision  variables  are  Q  and  Ptr- 

Theorem  8  (Scenario  2  with  k  =  1  or  k  =  oo )  IfT,  (5,  and  a  are  constant  and 
Q  is  unbounded,  the  unbounded  PRR  (t)  =  QQ*(t)  ■ 

Proof:  Consider  the  munition  and  sensor  craft  problems. 
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3-4-1  Munition  Problem.  For  k  =  1,  dynamic  Q,  and  dynamic  Ptr,  the 


states  becomes 

x  (t)  =  [  { [1  -  PFtr  (t)]  a  +  PTr  (t)  (3}  Q  (r)  dr 

Jo 

Z  (t)  =  [  aQ  (r)  [1  -  Pftr  (t)]  e~x{r)dT. 

Jo 

The  corresponding  L  is 

L  =  -PTr  ( t )  /?Q  (*)  e-*^. 

The  derivatives  of  the  states  are 


(3.117) 

(3.118) 


(3.119) 


x  (t)  —  {[1  ~  Pftr  (£)]  ol  +  Ptr  (t)  ft}  Q  (t)  (3.120) 

z(t)  =  aQ(t)[l-PFTR(t)}e-xW  (3.121) 


Equation  (2.12)  becomes 


[1  —  Pftr  (£)] 


_ Q  (t)  Ptr  (t) _ 

Q  (t)  Ptr  (t)  +  Qn  [1  —  Ptr  (£)] 


(3.122) 


Using  Eq.  (3.122),  the  Hamiltonian  becomes 


H=(  Xx- 


e  x)  QPtrP  +  ol 


Q2P 


TR 


QPtr  +  Qn(  1  —  Ptr) 


Xx  +  Xze~x]  (3.123) 


and  the  costate  differential  equations  become 


Q  (t)  Ptr  (t)  +  Qn  [1  —  Ptr  (£)] 
-Q(t)PTR(t)(3e-xW 

Xz  =  0 


(3.124) 

(3.125) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


0H  =  Q  /  (A,  -  e~x)  (5+  +  A"6  Qn°Q 


dP- 


TR 


[QPtR  +  Qn  (1  —  Ptr)Y 


(3.126) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Q  gives 


^  =  Pth{  (A*  -  e-)  13  +  a  [A,  +  Ke~’}  +  2Q„Q  (1  -  Ptr)  j.  {  } 

9Q  \V  ’  1  '  \QPTR  +  Qn(\-PTR)f  ‘  ' 


Ignoring  the  trivial  solutions  Ptr  =  0  and  Q  =  0,  the  necessary  conditions  become 


(Aa 


ca  o  ,  [Ax  +  Xze  ]  QnaQ  .  _ 
’  [ QPTR  +  Qn  (1  -  PTR)f 


(3.128) 


and 


(\  „—x\  Q  ,  „  i  \  „—x~\  Q  pTR  +  ZQnQ  (1  -  Ptr)  y  n  /Q1om 

(\x-e  )P  +  a[ Xx  +  Xze  J  rnp  n h>  bT-  r  =  9  (3-129) 

[W*TR  +  Vn  ft  —  RtR) \ 


Subtracting  Eq.  (3.128)  from  Eq.  (3.129)  gives 


a 


[A,  +  A,e-’]  Q2pTH  +  2Q..Q(l-Pra)  _  [Ax  +  A te~*]Q„aQ  =  Q 
1  {QPtr  +  QA3-  Ptr)}  \QPtr  +  Qn  (1  -  Ptr)} 


Grouping  like  terms,  yields 


olQ  Xx  +  A^e 


QPtr  ~  ZQuPtr  +  Qn 
QPtr  +  Qn  (1  —  Ptr)]2 


0 


resulting  in  the  necessary  condition 


[Ax  +  A^e  x  [QPtr  —  ZQuPtr  +  Qn]  —  0. 


(3.131) 


(3.132) 


Thus,  either 


PjTR  (t) 


Qn 

2Qn  -  Q*  (t) 


(3.133) 
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or 

\x(t)  =  -Ke~x{t)  (3.134) 

Ptr  (t)  =  2 q  QQ*(t)  the  same  necessary  condition  as  the  single  warhead  case  for 
Scenario  1.  Both  Scenario  1  and  Scenario  2  have  the  same  Poisson  held  of  false 
targets  with  density  a.  Thus,  it  makes  intuitive  sense  the  same  necessary  condition 
should  hold. 


3.4.2  Sensor  Craft  Problem.  For  k  =  00,  dynamic  Q,  and  dynamic  Ptr, 
the  states  become 


x  (t) 

y(t) 

z{t) 


aQ  (r)  [1  -  PFTr  (t)]  dr 


(3Q  (r)  PTr  (t)  dr 


rh—1 


_  [x  {t)]m  (t) 

^  ml 

m=0 

In  this  formulation,  f  is  used  and  L  —  0.  The  corresponding  <f  is 

(  h—  1 

*  =  -  '-E 


?i=0 


MliC  -,(T> 

n\ 


The  derivatives  of  the  states  are 


After  applying  Eq.  (3.122),  the  Hamiltonian  becomes 

n(m  —  1) 


H  =  a 


Q2P 


TR 


QPtr  +  Qn  (1  —  Ptr) 


Xx  +  A2 


xK 


L(m  —  1)! 


(3.135) 

(3.136) 

(3.137) 


(3.138) 


X  (t) 

=  aQ  (t)  [1  —  Pftr  (£)] 

(3.139) 

y(t) 

=  j3Q  (t)  PTr  (' t ) 

(3.140) 

z(t) 

=  a0(t)[l-P„s(t)]|h|^|e-*W 

(3.141) 

e~x  }  +  XySQPTR  (3.142) 
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and  the  costate  differential  equations  become 


Ax  (t)  = 


—X zaQ2  (t)  Ptr  ( t ) 

Q  (t)  Ptr  (t)  +  Qn  [1  —  Ptr  (£)] 

r(m-i)  i  /iqqtbfL 

x(t)  |  (m  —  1)!  [ 


A  y(t)  = 

Xz  (t)  = 


(3.143) 

(3.144) 

(3.145) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  _  (  [A-  +  A*  ([l£)|1})  e  *]  Q-aQ 

dP^R  Q  ]  [QPTr  +  Qn  (1  -  Ptr)}2  +  ^ 


(3.146) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Q  gives 


dH  n  \  x  />(*)]  t 

dQ~PTR  1“  A-  +  A-(  (m-1)!  I' 


_x  Q2Ptr  +  2 QnQ  (1  —  Ptr) 


+  A  y/3 


[QPtr  +  Qn  (1  -  Ptr)]2  V  J 

(3.147) 


Ignoring  the  trivial  solutions  Ptr  =  0  and  Q  =  0,  the  necessary  conditions  become 


1 V  ' - L - +  \  n  k  =  o 

[QPtr  +  Qn(  1  —  Ptr)]  ( 


(3.148) 


a  Ax  +  A2 


(m  —  1)! 


_x  Q2PTR  +  2QnQ(l-PTR) 
[QPtr  +  Qn(  1  —  Ptr)]2 


+  Xy(3  >  =  0  (3.149) 
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Subtracting  Eq.  (3.148)  from  Eq.  (3.149)  gives 


a 


A x  +  Xz 


(t)) 


(m— 1) 


(m  —  1)! 


Q2PTR  +  2QnQ(l-PTR) 
[QPtR  +  Qn(  1  —  PTR )]2 


As  +  Az 


(m— 1)! 


[Q-Pt.r  +  Qn  (1  —  -Pur). 


Qn&Q 

S - =  0  (3.150) 


Grouping  like  terms  gives 


aQ 


Xx  +  A2 


>  (*)] 


(to— 1) 


(m  —  1)! 


QPtr  ~  2QnPTR  +  Qn 
QPtr  +  Qn(  1  —  Ptr)]" 


=  0 


(3.151) 


resulting  in  the  necessary  condition 


Aa;  +  A^ 


>  (*)] 


(m— 1) ' 


(m  —  1)! 


[ QPtr  ~  2 QuPtr  +  Qn]  —  0. 


(3.152) 


Thus,  either 


Ptr  (t) 


Qn 

2Qn  -  Q*  (t) 


(3.153) 


or 

a-  « = (g;)  e"“,)  <3’i54) 

Ptr  (0  =  2Q  _Q*(t)  is  the  same  necessary  condition  as  the  single  warhead  case. 

Thus,  it  appears  the  necessary  condition  in  Eq.  (3.153)  holds  regardless  of  the 
number  of  warheads,  false  target  density,  or  maximum  allowable  P  (m  >  m).  This, 
of  course,  assumes  unbounded  design  variables.  If  Q  saturates  at  an  upper  or  lower 
bound,  the  fixed-Q  necessary  conditions  are  used. 
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3.5  Scenario  1:  Fixed  R,  Fixed  Q,  Variable  Ptr 

For  Scenario  7,  one  target  is  normally  distributed  among  a  Poisson  field  of 
false  targets  with  density  a,  which  is  assumed  constant.  Normally- distributed  refers 
to  a  circular  normal  distribution  with  standard  deviation  ay.  The  search  area  As 
is  a  circular  disc  of  radius  R,  which  is  prescribed.  The  air  vehicle  searches  the  disc 
using  concentric  annuli  of  radius  r  and  thickness  dr.  Each  annulus  is  covered  using 
a  constant  Q ;  therefore,  the  ROC  parameter  c  is  constant.  As  r  increases,  the  time 
it  takes  to  cover  the  corresponding  annulus  also  increases.  The  search  begins  at  the 
origin  of  the  disc  and  progresses  outward.  The  only  decision  variable  is  Ptr- 

3.5.1  Munition  Problem.  In  additional  to  assuming  i?,  Q ,  and  a  are 
constant,  the  constraint  on  P  (m  >  m)  is  initially  removed.  In  doing  so,  only  the 
following  state  is  needed: 

rr 

x  (r)  =  /  2irap  [1  —  Pftr  (p)]  dp.  (3.155) 

Jo 

The  corresponding  L  is 

r2 

L  =  -PTR  (r)  JJ-e~^e~x{r).  (3.156) 

( 7 1 
V 

The  derivative  of  the  state  is 

x  (r)  =  27TQ!r  [1  —  Pftr  (?’)]  •  (3.157) 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

r2  p 

H  =  -PTR-2e~^e-x  +  Xx2nar- - ™  (3.158) 

07  (1  -  c)  PTR  +  c 
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and  the  costate  differential  equation  becomes 


A*  (r)  =  -PTr  (r)  -^e  l(r). 


err 


(3.159) 


Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


d H  r  _ x 

= - -e  2ay  e  + 


\x2i rare 


OR 


TR  tJv 


[(1  -  c)  PTr  +  cY 


(3.160) 


Solving  —  0,  the  optimal  control  is 


Ptr  (r) 


c  ±  \j  \x  (r)  2iracGfle  L 

c  —  1 


x(r)+- 


(3.161) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1]. 
Taking  the  derivative  of  Eq.  (3.161)  with  respect  to  r  gives 


Ptr  (r) 


2(c-l) 


\ 


2tc  ac<Ty 


Xx  (r)  e 


x(r)+- 


x(r)+ 


X 


Ax  (r)  x  (r)  +  Ax  (r)  —  +  Ax  (r) 


Gi 


(3.162) 


Using  Eqs.  (3.157)  and  (3.159)  for  x(r )  and  Ax  (r)  yields 


Ptr  (r) 


\ 


irar * 


2c\x  (r) 


0-2  g 

UyC 


x(r)+- 


[ Ptr  (r)] 


1) 


Ax  (r)  e 


(3.163) 

Unlike  Scenarios  1  and  2,  there  is  no  discernable  trend.  Nonetheless,  Eq.  (3.161) 
shows  the  form  of  the  optimal  control. 
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Next,  the  effect  of  adding  a  constraint  on  false  target  attacks  is  examined. 
With  a  constraint  on  P  ( m  >  m),  the  states  are 


x  (r) 

y(r) 

z  ( r ) 


2nap  [1  -  PFTr  (p)]  dp 


Ptr  (p)  -xjC  dp 
ay 

2nap  [1  -  PFTr  (p)]  [1  -y  (p)]  e~x(p)dp. 


(3.164) 

(3.165) 

(3.166) 


The  corresponding  L  is  still 


L  =  -PTR(r)^-e  ^e~x{r). 


at 


(3.167) 


The  derivatives  of  the  states  are 


x  (r)  =  2nar  [1  -  P ft r  (r)] 

r  r 2 

y  (r)  =  PTR  (r)  —  e  *°v 
ay 

z  (r)  =  27 rar  [1  -  Pftr  (r)]  [1  -  y  (r)]  e~x(r) 


(3.168) 

(3.169) 

(3.170) 


After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

r 2  p 

H  =  (Ay  -  e~x)  PTR^e  +  27 rar- - — —  [A*  +  A2  (1  -  y)  e~x]  (3.171) 

CTy  0  —  c)  FTR  +  C 

and  the  costate  differential  equations  become 


A ^  (r) 


K  (r) 

K  (' r ) 


\~2nar 


Ptr  (r ) 


(1  -  c)  PTr  (r)  +  c 


[1  —  y(r)]e-’M 


—Ptr  (r)  Ae 


a ; 


y 


\~2i rar 


Ptr  (r ) 


(1  -  c)  PTR  (r)  +  c 


-,-x{r) 


(3.172) 

(3.173) 

(3.174) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


9H  f  r  [\x  +  Xz(l  -y)e  x)2tt arc 

— —  -  (Ay  -  e  )  —e  y  H  —  —  —2 

oPtr  ay  [(1  -  c)  PTr  +  cj 

Solving  —  0,  the  optimal  control  is 


(3.175) 


c  ± 


{Aa;(r)+A2[l-j/(r)]e  *(r)}-2?r accrue  2cts/ 


Ptr  ( r ) 


e  *(r)-Aj,(r) 


C-  1 


(3.176) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1]. 
While  taking  the  derivative  of  Eq.  (3.176)  with  respect  to  r  is  possible,  there  are  too 
many  terms  involved  to  discern  any  trends.  Nonetheless,  Eq.  (3.176)  shows  the  form 
of  the  optimal  control. 


3.5.2  Sensor  Craft  Problem.  The  sensor  craft  problem  does  not  condition 
P  (n  >  n)  on  the  number  of  previous  attacks.  For  Scenario  7, 


P  (n  =  1)  = 


rR  r  _  r2 

Ptr  (r)  —e  ^ dr . 


at. 


(3.177) 


Theorem  9  (Scenario  7  with  k  =  00  and  no  constraint  on  P  (m  >  m))  If  R  and 

a  are  constant,  the  optimal  constrained  solution  is  PfR  =  Pfrt,,,,,,  =  1- 

Proof:  With  no  constraint  on  P  (rri  >  m),  the  optimal  control  is  clearly  PfR  = 


P 


TRr, 


=  1. 


□ 

This  is  not  a  practical  solution,  since  all  false  targets  would  be  declared  targets  and 
attacked.  With  an  infinite  amount  of  warheads  and  no  constraint  on  false  target 
attacks,  there  is  no  reason  not  to  attack  every  object  encountered. 
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Next,  the  effect  of  adding  a  constraint  on  false  target  attacks  is  examined. 
The  sensor  craft  problem  does  not  condition  P  ( m  >  rh)  on  the  number  of  previous 
attacks. 


Theorem  10  (Scenario  7  with  k  —  oo  and  P  ( m  >  1)  =  b)  If  R  and  a  are  con¬ 
stant,  the  unbounded  PfR  is  monotonically  decreasing. 

Proof:  With  a  constraint  on  P  (m  >  rh ),  the  states  become 


nr 

X  (r)  =  /  2nap  [1  -  PFTr  (p)]  dp 

J  o 

(3.178) 

Ur)  =  1  AC1  MfVw 

w  ^  m\ 

m= 0 

(3.179) 

The  corresponding  L  is 

L  = 

The  derivatives  of  the  states  are 


—Ptr  (r) 


(3.180) 


x  (r)  =  2irar  [1  -  PFTr  (r)] 
z  (r)  =  27 rar  [1  -  PFTr  (r)] 


x  (r) 


(m— 1) 


D-x(r) 


(■ m  —  1)! 

After  applying  Eq.  (2.12)  for  the  ROC  model,  the  Hamiltonian  becomes 

Ptr 


(3.181) 

(3.182) 


H  =  Ptr  —r e  ^  +  27 rar- 

(?y  (1  -  c)  PTr  +  c 

and  the  costate  differential  equations  become 

Xz27 rar PTr  (r) 


Xx  +  A2 


x 


(m  —  1) 


(rh  —  1)! 


(3.183) 


(r)  = 


Az  (r)  =  0 


'  (m  —  1)  1 

J 

f[x(r)]<"*-iq 

x  (r) 

1 

[  (hi  —  1)1  j 

3-x(r) 


(3.184) 

(3.185) 
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Taking  the  partial  derivative  of  H  with  respect  to  the  decision  variable  Ptr  gives 


dH  r  - 

= - -e  2^  + 


2  27 rare  <  Ax  +  A 


(m— 1)! 


9Ptr 


[(1  —  c)  Ptr  +  c]' 


(3.186) 


Solving  —  0?  the  optimal  control  is 


^tr  (r)  — 


c  ±  \  ‘liuxco‘1  \  \x  (r)  +  A 


[x(r 


i(m-l) 


(m— 1)! 


e~x 0)  l  e^y 


(3.187) 


Only  the  “minus”  root  is  of  interest,  since  the  “plus”  root  puts  Ptr  outside  of  [0, 1]. 
Taking  the  derivative  of  Eq.  (3.187)  with  respect  to  r  gives 


Ptr  (r) 


-1 


2(c-l) 


2nac<jy 


x 


(r)  +  A 

r2 

A*  (r)  e2^  +  A2i  (r) 


[*(*■)] 


(m-l) 


^  I  (m-l)! 

(rh  —  1) 
a?  (r) 


f>—x(r)  V.  g2^J 

(m— 1) 


x  (r)] 


Ax  (r)  +  A~ 


( m  —  1)! 

X(r)]{7h~iy 
( m  —  1)! 


e-x{r)e^y 


k 

2 

e~x(r) 

-P—Ij  r 
e  v  —  ; 

1 

y  ) 

-2}  (3.188) 


Using  Eq.  (3.181)  for  x  (r)  with  Eq.  (2.12)  for  the  ROC  model  and  Eq.  (3.184)  for 
K  (r)  yields 


c  —  1 


<  0. 


(3.189) 

□ 


As  radius  increases,  the  amount  of  expected  false  target  encounters  in  an  annulus 
increases,  and  the  probability  of  encountering  the  target  in  an  annulus  decreases. 
With  infinite  warheads,  the  focus  becomes  meeting  the  constraint  on  false  target 
attacks.  Since  the  amount  of  expected  false  target  encounters  in  an  annulus  increases 
with  radius,  it  makes  intuitive  sense  to  tighten  the  operating  point  on  the  ROC  curve. 
That  is,  decrease  Ptr  as  radius  increases. 
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IV.  Applications 

In  this  chapter,  sample  applications  of  the  techniques  and  equations  provided  in  this 
dissertation  are  given.  First,  the  simplest  form  of  the  problem  where  all  variables 
are  fixed  is  examined.  Second,  the  effects  of  varying  Ptr  along  a  single  ROC  curve 
over  time  is  examined.  Finally,  varying  both  Ptr  and  Q  over  time  is  examined. 
Matlab’s®  fmincon  routine,  which  is  a  discrete-time,  gradient-based  solver  is  used. 
Analytic  conclusions  from  Chapter  III  are  used  to  verify  the  numerical  results. 

Effectiveness  measures  are  the  probability  of  at  least  one  target  attack,  P  (n  >  1), 
and  the  probability  of  at  least  one  false  target  attack,  P  (m  >  1).  The  following  con¬ 
strained  optimization  problem  is  solved: 

Max:  P  {n  >  1) 

Subj  to:  P  (m  >  1)  <  b 

Having  n  —  m  —  1  allows  one  to  compare  of  the  munition,  UCAV,  and  sensor  craft 
problems.  Thus,  the  sensitivity  of  effectiveness  to  the  number  of  warheads  can  be 
evaluated.  By  varying  b,  the  sensitivity  of  effectiveness  to  the  upper  bound  on  false 
target  attacks  can  be  evaluated.  Finally,  by  examining  a  variety  of  battle  space 
scenarios,  the  effect  of  the  target-to-false-target  ratio  being  fixed  or  a  function  of 
sensor  footprint  location  can  be  evaluated. 

4-1  Numerical  Solution  Method 

The  numerical  solution  method  comes  from  Bryson’s  gradient  algorithm  [7].  A 
synopsis  of  the  algorithm  is  now  given.  First  the  continuous  state  equations  are  put 
into  discrete  form  using  Euler’s  method.  The  problem  is  then  put  into  Mayer  form, 
where  the  state  vector  is  augmented  by  one  state  representing  the  cumulative  sum 
of  L  to  step  i.  The  problem  statement  then  becomes: 


4-1 


Min:  0  [s  (N)] 

subj  to:  s  (i  +  1)  =  f  [s  (i) ,  u  (i) ,  i],  i  —  0, . . . ,  N  —  1 
s  (0)  =  s0 
0  [s  (N)]  =  0 

where  N  is  the  number  of  time  or  distance  steps,  as  opposed  to  N,  which  is  the 
number  of  targets,  0  represents  terminal  constraints  and  is  some  function  of  the 
terminal  states.  In  this  case,  0  [s  (N)]  =  z  (N)  —  b.  Mayer  form  is  equivalent  to 
Bolza  form;  however,  Mayer  form  yields  simpler  expressions  and  is  easier  to  code 
for  numerical  solutions.  Matlab’s®  fmincon  routine  is  used,  so  one  does  not  have  to 
guess  step  size.  The  values  for  f,  fs,  fu,  0S,  and  0S  at  any  time  or  distance  step  must 
be  provided  in  a  subroutine.  The  algorithm  then  proceeds  as  follows: 

•  Guess  u  (i)  for  i  —  0, . . . ,  N  —  1 

•  Forward  propagate  the  state  equations  from  s  (0)  using  u  (z)  and  store  s  (i) 

•  Evaluate  0  and  0  and  set  A^  (N)  =  0^,  A  A  (N)  =  0j 

•  Backward  propagate  and  store  the  response  sequences  ( i )  and  (■ i ) 


A*  (i) 

=  faT  (i)  \*(i  +  1) , 

i  =  N  —  1, . . . ,  0 

A*  (i) 

=  faT  (i)  A*  (i  +  1) 

=  fuT(i)AUi  +  l) 

=  ft  (•)  A*  (i  +  1) 

The  process  is  repeated  until  the  terminal  constraint  error  and  the  gradient  sequence 
are  negligibly  small.  The  problem  becomes  a  parameter  optimization  involving  nc  x 
N  decision  variables,  where  nc  is  the  number  of  control  variables.  A  method  for 
nondimensionalization  was  proposed  in  a  handout  by  Pachter  and  is  included  in 
Appendix  C. 
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4-2  Fixed  Q,  Fixed  Ptr 

If  all  parameters  are  fixed,  the  results  from  Decker  [15]  or  Jacques  and  Pachter  [26] 
can  be  used.  To  illustrate  how  the  generalized  equations  derived  in  this  dissertation 
collapse  to  the  results  in  [15]  and  [26],  consider  a  Scenario  3  problem  and  the  calcu¬ 
lation  of  the  first  elemental  probability  P\.  From  Table  2.5, 


N- 1 

*=  £ 

j=n— 1 


N~l\(  3 

j  )  \n  -  1 


’  ~j~ptr  (t)  Q  (t)  dr 
Uo  ab 


n—  1 


X 


1  j-(n- 1) 


IJ  0  ab 


—Q  (r)  dr  - 


-j- ptr  (r)  Q  (r)  dr 


/  o  vi b 


x 


ft  1 

N-l-j 

1  -  /  -j-Q{r)dT 
L  Jo  ab  J 

l 

J 

~7  PTR  (t)  Q  (t)  dt  (4.1) 

Vi  B 


Assuming  fixed  parameters  gives 


N- 1 

*  =  E 

j=n- 1 


N~l\(  j 

j  )  V  n  -  1 


P  A 

Ptra~ 


i  n— 1  r 


B 


(1  -  PTr )  -t— 
A-b 


1  j-(n- 1) 


X 


A  ' 


N-l-j 


N 

-r—PrRQdt  (4.2) 

Vi  B 


Decker  showed  encounters  could  be  decoupled  from  classifications  resulting  in 

N-l-j 


N- 1 

*  =  £ 
j=n~  1 


'  A  ' 

j 

[j  "4  1 

IV  3  ) 

,ab_ 

As_ 

X 


J 


,  I  [PtrT  1  [(1  —  PTr)Y  (n  ^  \  -R—PTRQdt  (4.3) 
n-lj  I  Ab 


He  then  applied  a  truncated  binomial  conversion  giving 

Pi  = 


(N-l\ 

A  ' 

n—  1 

„  N  A  ' 

U-v 

PtRAb_ 

(1  -  PTr )  -7— 
ViB. 

N—n 


N 

-r—PrRQdt 
Vi  B 


(4.4) 
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Letting  N  —  1  and  k  —  1  gives 


P 1  —  -j —PrRQdt 

A-  B 


(4.5) 


which  matches  the  results  derived  by  Jacques  and  Pachter  [26]. 

A  fixed-parameter  problem  has  advantages.  The  operator  simply  sets  Ptr 
at  the  beginning  of  the  mission.  No  sophisticated  scheduling  software  or  hardware 
is  needed.  Finding  the  optimal  Ptr  can  be  done  looking  at  a  plot  or  table.  For 
example,  in  Scenario  1  with  k  —  1,  the  effectiveness  measures  become 

—  e_(1~ 

(1  —  PF1 

and 


Pftr)ccQT 


a  A 


B 


(4.6) 


P  (n  >  1)  =  P; 


TR~ 
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P  (m  >  1) 


1  —  Ptr 


^  _  g— (1  ~Pftr)ciQT 

(1  —  Pftr )  OlAr 


PtrQT\ 
Ab  J 


g— (1  -Pftr)(*QT 


(4.7) 


Equation  (2.12)  can  be  used  for  the  ROC  model  with  c  =  100.  Assume  the  munition 
covers  Ab  by  time  T.  Thus,  QT  =  Ab-  One  of  three  mutually  exclusive  events 
are  possible.  The  munition  either  attacks  a  target,  attacks  a  false  target,  or  attacks 
nothing  and  self  destructs.  Figure  4.1  shows  the  outcome  probabilities  versus  Ptr 
when  c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  hr. 

If  the  constraint  on  a  false  target  attack  is  removed,  Figure  4.1  shows  the  best 
unconstrained  solution  is  PFr  =  0.723  with  a  corresponding  P(n>  1)*  =  0.535 
and  P  (rri  >  1)*  =  0.318.  If  P  (m  >  1)  is  bounded  by  b  =  0.2,  the  best  constrained 
solution  is  PFr  =  0.563  with  a  corresponding  P  (n  >  1)*  =  0.483.  In  Figure  4.1, 
the  outcome  probability  functions  are  smooth  and  well-behaved.  The  function  for 
P  (n  >  1)  has  one  peak.  The  function  for  P  (rri  >  1)  is  monotonically  increasing, 
so  any  constrained  solution  will  be  unique.  Finally,  the  function  for  the  probability 
of  no  attack  always  starts  at  (0,1),  then  monotonically  decreases.  Since  the  three 
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Figure  4.1  Outcome  probabilities  versus  probability  of  a  target  report.  Scenario  1 
with  constant  parameters,  c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  [hr], 

outcome  probabilities  are  mutually  exclusive  in  a  munition  problem,  they  sum  to 
one  for  each  PTR  considered. 
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The  problem  can  be  solved  for  a  number  of  b  values  and  a  plot  of  PRR  versus 
b  can  be  generated.  If  b  >  0.318,  the  unconstrained  solution  PrRu  =  0.723  is  used. 
Figure  4.2  illustrates  the  sensitivity  of  the  solution  to  changes  in  b.  Figure  4.2  is  the 
plot  an  operator  needs  to  set  Ptr  at  the  beginning  of  a  mission,  assuming  b  has  been 
determined  by  the  commander. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  4.2  Optimal  probability  of  a  target  report  versus  maximum  allowable  prob¬ 
ability  of  a  false  target  attack,  b.  Scenario  1  with  constant  parameters, 
c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  [hr]. 

Similar  results  are  expected  for  Scenarios  2-4,  where  target-to-false-target  ra¬ 
tios  are  the  same  regardless  of  sensor  footprint  location.  For  Scenarios  5-7,  target-to- 
false-target  ratios  depend  on  the  sensor  footprint  location.  Therefore,  the  function  for 
P  {m  >  1)  is  not  always  monotonically  increasing.  Hence  solutions  are  not  unique. 
For  example,  Figure  4.3  shows  the  function  for  P  (m  >  1)  initially  increases,  then 
decreases,  then  increases  again  for  a  Scenario  7  problem  with  k  —  1.  For  values  of 
P  (m  >  1)  above  the  local  minimum  at  P  (m  >  1)  =  0.385,  the  unconstrained  solu¬ 
tion  is  used.  On  a  plot  of  PRR  versus  b,  a  discontinuous  jump  appears  as  shown  in 
Figure  4.4. 
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Figure  4.3  Outcome  probabilities  versus  probability  of  a  target  report.  Scenario 
7  with  constant  parameters,  c  =  100,  a  =  5  [1/km2],  a2  =  0.25  [km2], 
and  R  =  5  [km]. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  4.4  Optimal  probability  of  a  target  report  versus  maximum  allowable  prob¬ 
ability  of  a  false  target  attack.  Scenario  7  with  constant  parameters, 
c  =  100,  a  —  5  [1/km2],  a2  =  0.25  [km2],  and  R  —  5  [km]. 
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4-3  Fixed  Q,  Variable  Ptr 

Fixed-parameter  problems  are  readily  solvable  and  simple  to  employ  opera¬ 
tionally.  Improvement  when  parameters  are  allowed  to  vary  is  now  examined.  To 
begin,  the  case  where  Ptr  is  allowed  to  vary  while  Q  remains  hxed  is  considered. 
Having  a  hxed  area  coverage  rate  means  a  single  ROC  curve  is  used.  Moving  along 
a  single  ROC  curve  is  done  by  adjusting  the  sensor  threshold,  which  is  tantamount 
to  adjusting  Ptr- 

Consider  Scenario  1  with  k  =  1  and  Eq.  (2.12)  for  the  ROC  model  with  c  =  100. 
Assume  the  munition  covers  Ar  by  time  T.  Thus,  QT  —  Ar.  The  optimal  control 
problem  can  be  solved  for  a  number  of  b  values.  Figure  4.5  shows  optimal  Ptr 
schedules  for  three  constraint  values.  In  Chapter  III,  the  optimal  Ptr  was  shown  to 
increase  over  time  when  the  constraint  on  false  target  attacks  is  removed.  Figure  4.5 
shows  the  optimal  Ptr  increases  over  time  even  when  a  constraint  on  false  target 
attacks  exists.  The  increases  diminish  as  constraint  values  are  tightened  (lowered). 


Figure  4.5  Optimal  probability  of  a  target  report  versus  time.  Scenario  1  with 
k  =  1,  c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  [hr]. 
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The  problem  can  be  solved  for  a  number  of  constraint  values  and  the  objective 
function  values  can  be  compared  to  those  obtained  using  fixed  parameters.  In  the 
system  operating  characteristic  shown  in  Figure  4.6,  there  is  negligible  improvement 
for  constraint  values  below  0.15  and  marginal  improvement  (up  to  3.4  percent)  for 
constraint  values  above  0.15.  Depending  on  the  application,  this  improvement  may 
or  may  not  be  worth  the  added  complexity. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  4.6  System  Operating  Characteristic.  Dynamic  versus  constant  Ptr-  Sce¬ 
nario  1  with  k  —  1,  c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  [hr]. 

Up  to  now,  only  the  munition  and  sensor  craft  problems  have  been  considered. 
It  was  shown  in  Chapter  III  that  the  optimal  Ptr  increases  over  time  for  the  Scenario 
1  munition  problem,  when  the  constraint  on  false  target  attacks  is  removed.  It  was 
also  shown  in  Chapter  III  that  the  optimal  Ptr  is  constant  when  k  =  oo  and,  for 
Scenario  1,  is  given  by 


Ptr  {i) 


— chi  [1  —  b] 

(1  —  c)  In  [1  —  b]  +  olQT 


(4.8) 
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Now  consider  the  effect  of  multiple  warheads.  Figure  4.7  shows  the  system  operating 
characteristic  for  k  =  l,2,6,oo  and  rh  —  h  —  1.  The  results  for  k  =  oo  come 
analytically  from  Eq.  (4.8).  The  k  =  1  results  do  not  exactly  conform  to  other 
k  values  which  are  bound  by  k  =  oo.  The  munition  problem  has  three  mutually 
exclusive  outcomes;  whereas,  the  outcomes  for  the  k  >  1  problems  are  not  mutually 
exclusive.  Nonetheless,  Figure  4.7  can  be  used  by  commanders  to  make  decisions 
about  constraint  values  and  number  of  warheads.  For  b  <  0.20,  there  is  not  much 
to  be  gained  by  adding  warheads.  In  this  single-target  scenario,  there  are  more  false 
targets  present,  which  means  the  probability  of  the  first  attack  being  a  false  target 
attack  is  relatively  high.  Thus  it  does  not  matter  how  many  warheads  there  are.  The 
constraint  is  met  by  the  Erst  attack.  For  b  <  0.20,  one  can  simply  use  the  constant 
closed-form  Ptr  given  by  Eq.  (4.8).  In  Figure  4.7,  it  is  clear  the  benefits  of  more  than 
6  warheads  is  negligible  throughout  most  of  the  range  of  b.  These  types  of  insights 
are  possible  using  the  techniques  and  equations  provided  in  this  dissertation. 


Maximum  allowable  probability  of  attacking  at  least  1  false  target 

Figure  4.7  System  Operating  Characteristic.  Dynamic  Ptr ■  Scenario  1  with 
c  =  100,  aQ  =  50  [1/hr],  and  T  =  0.5  [hr]. 
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Consider  Scenario  7  with  k  —  1  and  Eq.  (2.12)  for  the  ROC  model  with  c  =  100. 
The  target-to-false-target  ratio  depends  on  the  sensor  footprint  location.  Figure  4.8 
shows  optimal  Ptr  schedules  for  three  constraint  values.  No  conclusions  were  made 
for  the  Scenario  7  munition  problem  in  Chapter  III.  Figure  4.8  shows  as  b  tightens 
(decreases),  the  sensor  is  “turned  off’  earlier  in  the  search. 


Figure  4.8  Optimal  probability  of  a  target  report  versus  radius.  Scenario  7  with 
k  =  1,  c  =  100,  a  =  5  [1/km2],  a2  =  0.25  [km2],  and  R  =  5  [km]. 

The  objective  function  values  of  the  dynamic  solutions  can  be  compared  to 
those  of  the  constant  solutions.  Figure  4.9  shows  the  poor  performance  of  the 
constant- Ptr  solution  when  b  <  0.385.  Unlike  the  Scenario  1  example  where  dynamic- 
Ptr  solutions  were  only  a  few  percent  better  than  constant- Ptr  solutions,  there  are 
substantial  improvements  when  b  <  0.385.  Since  the  dynamic  problem  allows  the 
sensor  to  be  “turned  off’  at  some  point,  solutions  are  similar  to  those  at  smaller  R 
values.  How  far  out  in  radius  to  search  is  clear.  The  dynamic  problem  “recognizes” 
the  benefit  of  staying  close  to  the  origin  when  a  target  is  normally  distributed.  The 
fixed-parameter  problem  requires  an  arbitrary  pick  for  R.  Thus,  performance  is  poor 
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when  R  is  chosen  too  large.  To  improve  performance,  one  would  need  to  make  R  a 
parameter  and  run  a  sensitivity  study. 


Maximum  Allowable  Probability  of  a  False  Target  Attack 

Figure  4.9  System  Operating  Characteristic.  Dynamic  versus  constant  Ptr .  Sce¬ 
nario  7  with  k  =  1,  c  =  100,  a  =  5  [1/km2],  a2  =  0.25  [km2],  and  R  =  5 
[km] . 


The  dynamic- Ptr  solutions  have  been  compared  to  fixed- Ptr  solutions.  The 
amount  of  improvement  varies  based  on  problem  parameters  and  whether  the  target- 
to-false-target  ratio  is  constant  throughout  the  search.  Next  the  requirement  of  fixed 
area  coverage  rate  is  relaxed. 
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4-4  Variable  Q,  Variable  Ptr 

Allowing  Q  to  vary  while  keeping  T  fixed  introduces  another  trade  space.  For 
a  given  b ,  larger  Q  means  more  area  is  covered  in  time  T  resulting  in  more  target 
encounters.  However,  there  is  less  chance  of  correctly  classifying  the  targets.  Con¬ 
versely,  for  a  given  b,  smaller  Q  means  less  area  is  covered  in  time  T  resulting  in  less 
target  encounters.  However,  there  is  more  chance  of  correctly  classifying  the  targets. 
In  Section  2.3,  a  way  was  presented  to  model  sensor  performance  by  making  the 
ROC  parameter  c  a  function  of  area  coverage  rate  Q.  Specifically, 

=  W)  =  |  (4.9) 

where  Qn  is  a  scaled  nominal  area  coverage  rate. 

Allowing  both  Q  and  Ptr  to  vary  means  the  trajectory  on  the  Ptr  versus 
(1  —  Pftr )  plane  is  bound  by  a  region  instead  of  a  single  ROC  curve.  The  region  is 
formed  by  a  Qmin  ROC  curve  and  a  Qmax  ROC  curve.  Qmin  might  be  tied  to  an  air 
vehicle’s  stall  speed,  which  is  available  in  the  air  vehicle  performance  specification. 
Qmax  might  be  tied  to  an  air  vehicle’s  maximum  speed,  girnbal  limits,  or  search  area. 
Assuming  the  real  Qmax  would  result  in  As  >  Ar,  an  artificial  limit  needs  to  be 
invoked.  Otherwise,  probabilities  could  be  greater  than  one.  To  ensure  As  <  Ar, 
one  can  define  Qmax  as 

Qmax  =  (4.10) 

In  Chapter  III  it  was  shown  that  when  Q  is  unbounded, 


^tr  (Q 


Qn 

2Qn  ~  Q*  ( t ) 


(4.11) 


regardless  of  a,  b ,  or  number  of  warheads.  Further,  for  Q*  e  (0,Qn),  P^R  G  (l,  l). 
In  other  words, 

Ptr  (t)  >  0-5  (4.12) 
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For  Qn  »  Q*,  PRR  «  0.5  for  all  time.  If  Q  is  at  its  lower  or  upper  bound  for 
a  particular  region,  then  it  can  be  considered  fixed  in  that  region.  Therefore,  one 
expects  to  see  fixed-Q  trends  in  regions  where  Q  is  at  a  bound. 

To  illustrate  what  is  going  on  in  the  dynarnic-Q,  dynamic- Ptr  problem,  one 
can  plot  the  trajectory  of  the  operating  point  on  the  ROC  plot  for  Scenario  1  with 
k  =  1,  a  =  2,  T  =  0.5,  Qn  =  1000,  Qmtn  =  10,  and  Qmax  =  20.  Figure  4.10  shows 
four  trajectories  for  different  values  of  b  on  one  plot.  Figure  4.11  shows  the  four 
trajectories  on  4  different  plots.  When  the  constraint  on  false  target  attacks  is  tight 
(low),  the  air  vehicle  must  fly  at  Qmm  to  get  the  best  sensor  performance,  and  PRR 
increases  with  time.  When  b  =  0.09,  the  air  vehicle  starts  at  Qmin ,  but  eventually 
gets  off  the  bound.  As  predicted,  this  occurs  at  a  Ptr  value  slightly  above  0.50. 
When  b  =  0.15,  the  air  vehicle  flies  the  entire  mission  with  Q  unbounded  and  the 
corresponding  PRR  ps  0.5.  When  b  =  0.25,  the  air  vehicle  starts  with  an  unbounded 
Q ,  but  then  hits  the  upper  bound. 


Figure  4.10  ROC  plot.  Scenario  1  with  dynamic-Q,  dynamic- Ptr,  k  —  1,  a  —  2 
[1/km2],  Qn  =  1000,  Qmin  =  10,  Qmax  =  20,  and  T  =  0.5  [hr]. 
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Figure  4.11  ROC  plot  broken  out  by  upper  bound.  Scenario  1  with  dynamic- 
Q,  dynamic- PTr,  k  =  1,  a  =  2  [1/km2],  Qn  =  1000,  Qmin  =  10, 
Qmax  =  20,  and  T  =  0.5  [hr]. 

Time  is  implicit  to  the  trajectories  in  Figures  4.10  and  4.11.  If  a  trajectory  is 
along  the  Qmin  or  Qmax  ROC  curve,  time  increases  as  the  trajectory  moves  up  and 
to  the  right.  If  a  trajectory  is  in  a  region  where  Q  is  unbounded,  time  increases  as 
the  trajectory  moves  to  the  right. 
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Optimal  schedules  of  the  design  variables  Q  and  Ptr  are  shown  in  Figures  4.12(a) 
and  4.12(b).  When  Q*  is  at  either  of  its  bounds,  PRR  behaves  similarly  to  the  fixed-Q 
results  shown  previously.  That  is,  PRR  is  monotonically  increasing  with  flatter  curves 
at  tighter  constraint  levels.  When  Q*  is  unbounded,  PRR  ~  0.5  as  predicted  and  Q* 
is  increasing.  This  represents  the  strategy  seen  before  of  staring  conservatively,  then 
increasing  aggressiveness  with  time. 


Time  Time 

(a)  Q  schedule  (b)  Ptr  Schedule 

Figure  4.12  Optimal  Schedules  for  the  Design  Variables.  Scenario  1  with  dynamic- 
Q,  dynamic- PTR,  k  =  1,  a  =  2  [1/km2],  Qn  =  1000,  Qrntn  =  10, 
Qmax  =  20,  and  T  =  0.5  [hr]. 
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A  better  indication  of  what  is  going  on  comes  from  looking  at  the  schedule  for 
(1  —  Pftr )•  Figure  4.13  shows  (1  —  Pftr )  is  always  smooth  and  increasing.  These 
are  the  optimal  functions  required  for  meeting  the  constraint  on  false  target  attacks. 
Q  and  Ptr  must  combine  at  each  time  such  that  (1  —  Pftr)  falls  on  one  of  the 
functions  in  Figure  4.13. 


Time 


Figure  4.13  Optimal  (1  —  Pftr)  schedule.  Scenario  1  with  dynamic- Q,  dynamic- 
PTR,  k  =  1,  a  =  2  [1/km2],  Qn  =  1000,  Qmm  =  10,  Qmax  =  20,  and 
T  =  0.5  [hr]. 
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The  tradeoff  situation  where  lower  Q  values  provide  better  sensor  performance 
at  the  expense  of  less  area  searched  was  mentioned  earlier.  All  this  depends  on  b. 
Figure  4.14  shows  the  transition  from  the  minimum  search  area  to  the  maximum 


search  area. 


Maximum  allowable  probability  of  attacking  at  least  1  false  target 


Figure  4.14  Area  searched  versus  constraint  on  false  target  attacks.  Scenario  1 
with  dynamic-Q,  dynamic- TV  r,  k  —  1,  a  —  2  [1  /km2] ,  Qn  =  1000, 
Qmin  =  10,  Qmax  =  20,  and  T  =  0.5  [hr]. 


4-18 


As  with  previous  results,  it  is  beneficial  to  compare  the  dynamic-parameter 
results  to  fixed-parameter  results.  This  can  be  done  two  ways.  First,  one  can  require 
a  fixed-parameter  solution  at  each  b  value,  whereby  each  parameter  is  optimized  at 
each  b  value.  Second,  one  can  fix  one  of  the  parameters  and  allow  only  the  other 
parameter  to  be  optimized  at  each  b  value.  Figure  4.15  compares  dynamic  and 
fixed  solutions  when  each  parameter  is  optimized  at  each  b  value.  Only  marginal 
improvement  is  seen  for  the  dynamic-parameter  problem. 


Maximum  allowable  probability  of  attacking  at  least  1  false  target 

Figure  4.15  System  Operating  Characteristic.  Dynamic  versus  constant  solutions. 

Scenario  1  with  k  —  1,  a  —  2  [1/km2],  Qn  =  1000,  Qmin  —  10, 
Qmax  =  20,  and  T  =  0.5  [hr].  Fixed-parameter  problem  had  both  Q 
and  Ptr  optimized  at  each  b 
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In  Figure  4.16,  operating  points  for  different  b  values  are  plotted  on  a  ROC 
plot  when  both  Q  and  Ptr  are  fixed  but  optimized  at  each  b  value.  The  operating 
points  were  determined  using  Eq.  (4.6)  and  Eq.  (4.7)  with 


(1  —  Pftr ) 


QPtr 

QPtr  +  Qn  (1  —  Ptr ) 


(4.13) 


A  grid  search  was  used  at  each  b  value  to  find  the  optimal  fixed  operating  point, 
whereby  Q  and  Ptr  were  both  varied  from  their  minimum  to  maximum  values. 
The  upper  bound  b  was  incremented  by  0.01  to  0.40.  The  trajectory  is  similar  to 
the  dynamic- Q,  dynamic- trajectory  in  Figure  4.11.  For  low  values  of  b,  the 
air  vehicle  flies  at  Qmin ■  Fc»r  high  values  of  b ,  the  air  vehicle  flies  at  Qmax ■  A 
“bridge”  occurs  somewhere  above  Ptr  =  0.5.  Though  not  dynamic,  this  behavior  of 
Ptr  ~  0.5  when  Q*  is  unbounded  matches  the  dynamic  analytical  predictions.  This 
kind  of  insight  can  help  designers  and  mission  planners. 


Figure  4.16  ROC  plot.  Scenario  1  with  k  —  1,  a  —  2  [1/km2],  Qn  =  1000, 
Qmin  =  10,  Qmax  =  20,  and  T  =  0.5  [hr].  Fixed-parameter  problem 
had  both  Q  and  Ptr  optimized  at  each  b 
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Next  consider  the  case  when  one  of  the  two  design  variables  is  fixed  for  all 
b  values.  Figure  4.17  shows  when  Q  —  15  for  all  b  values,  the  improvements  for 
the  dynamic-parameter  problem  are  more  pronounced.  Hence,  there  are  benefits 
to  allowing  both  Q  and  Ptr  to  be  optimized  at  each  b  value,  even  when  they  are 
required  to  be  fixed. 


Maximum  allowable  probability  of  attacking  at  least  1  false  target 

Figure  4.17  System  Operating  Characteristic.  Dynamic  versus  constant  solutions. 

Scenario  1  with  k  —  1,  a  —  2  [1/km2],  Qn  =  1000,  Qmin  =  10, 
Qmax  =  20,  and  T  =  0.5  [hr].  Fixed-parameter  problem  had  Q  =  15 
[km2 /hr]  and  Ptr  optimized  at  each  b 
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The  effects  of  additional  warheads  remains  to  be  examined.  All  the  example 
problems  thus  far  have  been  either  Scenario  1  or  Scenario  7.  Therefore  an  example 
using  Scenario  2  is  in  order,  since  analysis  of  Scenario  2  was  also  done  in  Chapter  III. 
Consider  Scenario  2  with  a  =  1  [1/km2],  (3  =  0.1  [1/km2],  Qn  =  1000  [km2/hr], 
Qmin  =  10  [km2 /hr],  Qmax  =  20  [km2/hr]  and  T  =  0.5  [hr].  Analysis  showed  if 
k  —  oo,  m  —  1,  and  P  ( m  >  1)  =  b  with  T,  Q,  and  a  constant,  then 


p *  (f)  =  ~Qn  In  (1  ~  b) 

TR  u  (Q  ~  Qn)  In  (1  -  b)  +  aQ2T 


(4.14) 


assuming  c  =  % . 
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Equation  (4.14)  can  be  used  to  gain  a  rough  idea  of  when  it  is  better  to  fly  at 
Qmin  and  when  it  is  better  to  fly  at  Qmax-  Figure  4.18  shows  a  crossover  in  the  trade 
space  around  b  =  0.10.  At  tight  constraints  on  false  target  attacks,  the  system  must 
be  on  the  best  ROC  curve  and  uses  Qmin.  As  the  constraint  loosens,  the  system  can 
transition  to  a  worse  ROC  curve  with  the  idea  of  encountering  more  targets  by  using 
a  larger  Q.  While  Figure  4.18  applies  to  k  —  oo  and  fixed  Q,  the  dynamic  problem 
should  approach  the  k  —  oo,  Q  —  Qmax  at  high  values  of  b. 


Figure  4.18  System  Operating  Characteristic.  Scenario  2  with  Q  fixed  at  either 
Qmin  or  Qmax:  Qn  =  1000,  k  =  OO,  OL  =  1  [1/klU2],  /3  =  0.1  [l/km2], 
and  T  =  0.5  [hr]. 
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Figure  4.19  shows  the  trade  space  for  the  dynamic- Q,  dynamic- Ptr  problem 
for  k  =  1,2, 4, 6.  As  seen  before,  increasing  warheads  for  a  tight  constraint  on  false 
target  attacks  provides  no  value.  As  the  constraint  loosens,  increasing  warheads 
eventually  adds  value.  For  this  problem,  there  is  minimal  benefit  to  adding  more 
than  4  warheads.  For  the  most  part,  the  4- warhead  solution  follows  the  asymptote 
generated  by  the  k  =  oo,  Q  =  Qmax  solution  from  Eq.  (4.14). 


Maximum  allowable  probability  of  attacking  at  least  1  false  target 

Figure  4.19  System  Operating  Characteristic.  Scenario  2  with  dynamie-Q, 
dynamic  -Ptr,  a  =  1  [1/km2],  (3  =  0.1  [1/km2],  and  T  =  0.5  [hr]. 
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One  can  look  to  trajectories  on  a  ROC  plot  to  gain  further  insight.  Figure  4.20 
shows  4  trajectories  for  4  different  b  values  when  k  —  1.  When  the  constraint  on 
false  target  attacks  is  tight  (low),  the  air  vehicle  must  fly  at  Qmin  to  get  the  best 
sensor  performance,  and  PRR  increases  with  time.  When  b  =  0.05,  the  air  vehicle 
starts  at  Qmim  but  quickly  gets  off  the  bound.  As  predicted  in  Chapter  III,  this 
occurs  at  a  Ptr  value  slightly  above  0.50.  When  b  =  0.10,  the  air  vehicle  flies  the 
entire  mission  with  Q  unbounded  and  the  corresponding  PRR  «  0.5.  When  b  =  0.15, 
the  air  vehicle  starts  with  an  unbounded  Q,  but  then  hits  the  upper  bound.  It  is 
interesting  to  note  that  the  crossover  from  Qmin  to  Qmax  occurred  around  the  same  b 
value  as  the  k  =  oo,  fixed-Q  results  in  Figure  4.18.  While  Ptr  varies  in  Figure  4.20, 
the  variation  is  small. 


Figure  4.20  ROC  plot.  Scenario  2  with  dynarnic-Q,  dynamic- Ptr,  a  =  1  [1/km2], 
/ 3  =  0.1  [1/km2],  T  =  0.5  [hr],  and  k  —  1. 
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Figure  4.21  shows  4  trajectories  for  4  different  b  values  when  k  =  2.  The 


variations  in  Ptr  are  even  smaller.  In  addition,  at  each  b,  the  magnitudes  of  Ptr 
and  Q  are  smaller  than  the  k  —  1  results.  This  makes  intuitive  sense  as  increasing 
warheads  requires  a  more  conservative  operating  point  to  meet  the  constraint  on  false 
target  attacks.  In  fact,  for  the  given  problem,  adding  warheads  actually  degraded 
performance  for  tight  constraints,  albeit  by  very  small  amounts  not  readily  visible 
in  Figure  4.19. 


Figure  4.21  ROC  plot.  Scenario  2  with  dynamic-Q,  dynamic- Ptr,  ot—1  [1/km2], 
/3  =  0.1  [1/km2],  T  =  0.5  [hr],  and  k  =  2. 
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V.  Conclusions  and  Recommendations 


5. 1  Summary 

The  purpose  of  this  research  was  to  analytically  quantify  operational  effective¬ 
ness  of  airborne  systems  performing  search  and  destroy  missions.  The  primary  con¬ 
tribution  to  the  field  is  the  ability  to  handle  time- varying  parameters,  which  broadens 
the  possible  battle  space  scenarios  and  enables  improved  performance  through  func¬ 
tion  optimization.  Prior  to  this  research,  one  could  only  analytically  solve  problems 
where  area  coverage  rate,  sensor  threshold,  target  density,  and  false  target  density 
where  fixed.  Now  one  can  analytically  solve  problems  where  any  or  all  of  the  above 
parameters  vary  with  time. 

Two  types  of  search  patterns  were  considered.  The  first  type  was  the  linear 
search  pattern,  which  can  be  used  to  cover  any  linear-symmetric  battle  space.  The 
second  type  involved  concentric  annuli  emanating  from  the  origin,  which  can  be  used 
to  cover  any  circular  battle  space.  Uniform,  Poisson,  and  normal  probability  distrib¬ 
utions  were  used  to  model  target  and  false  target  encounters.  Expressions  needed  to 
calculate  the  probability  of  at  least  n  target  attacks,  P  (■ n  >  h),  and  the  probability 
of  at  least  m  false  target  attacks,  P  (m  >  rh ),  were  derived  for  the  following  seven 
battle  space  scenarios: 

Scenario  1  consisted  of  a  single  target  uniformly  distributed  in  a  linear- 
symmetric  battle  space  among  a  Poisson  field  of  false  targets. 

Scenario  2  consisted  of  a  linear-symmetric  battle  space  with  a  Poisson  field 
of  targets  and  a  Poisson  field  of  false  targets. 

Scenario  3  consisted  of  N  targets  uniformly  distributed  in  a  linear-symmetric 
battle  space  among  a  Poisson  field  of  false  targets. 

Scenario  4  consisted  of  N  targets  and  M  false  targets  uniformly  distributed 
in  a  linear-symmetric  battle  space. 
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Scenario  5  consisted  of  N  targets  normally  distributed  in  a  circular  battle 
space  among  a  Poisson  field  of  false  targets.  Normally  distributed  referred  to  a 
circular-normal  distribution  centered  at  the  origin  with  a  target  standard  deviation 
<JZ. 

Scenario  6  consisted  of  N  targets  and  M  false  targets  normally  distributed  in 
a  circular  battle  space.  Normally  distributed  referred  to  circular-normal  distributions 
centered  at  the  origin  with  a  target  standard  deviation  az  and  a  false  target  standard 
deviation  ay. 

Scenario  7  consisted  of  a  single  target  normally  distributed  in  a  circular  battle 
space  among  a  Poisson  held  of  false  targets. 

Sensor  performance,  in  terms  of  a  system’s  ability  to  autonomously  distinguish 
between  targets  and  false  targets,  came  from  a  Receiver  Operating  Characteristic 
(ROC)  curve.  A  representative  ROC  curve  model  was  used;  whereby,  a  parameter  c 
dictated  the  locus  of  operating  points.  The  parameter  c  depended  inversely  on  the 
area  coverage  rate  Q.  Based  on  a  user-supplied  constant  Qn ,  sensor  performance  was 
completely  determined  by  Q  and  the  probability  of  a  target  report  Ptr-  Together, 
Q  and  Ptr  establish  the  system  operating  characteristic. 

Optimal  control  problems  were  formulated  to  maximize  P  (n>  n )  subject  to 
P  {m  >  m)  <  b,  where  b  was  set  by  the  user.  The  decision  variables  were  Q  and 
Ptr ■  Due  to  coupling  of  the  state  and  costate  equations,  the  two  point  boundary 
value  problems  were  solved  numerically  using  Bryson’s  gradient  algorithm  [7]  and 
Matlab’s®  fmincon  routine. 

Prior  to  numerically  solving  sample  problems,  analytic  solutions  were  provided 
for  the  munition  and  sensor  craft  problems  for  Scenarios  1,  2,  and  7.  Using  these 
tractable  instances,  the  following  conclusions  were  made  and  used  to  verify  numerical 
results: 


5-2 


1.  For  Scenarios  1  and  2  with  k  —  1  and  no  constraint  on  P  ( m  >  rh ),  if  T,  Q, 
a,  and  f3  are  constant,  the  unbounded  PRR  is  monotonically  increasing  with  time. 

2.  For  Scenarios  1,  2,  and  7  with  k  =  oo,  n  —  1,  and  no  constraint  on 
P  (to  >  to),  if  T  (or  R),  Q,  a,  and  (5  are  constant,  the  bounded  PRR  =  PRR max  =  1. 

3.  For  Scenarios  1  and  2  with  k  =  oo  and  P  (to  >  to)  =  b ,  if  T,  Q,  a,  and  /3 
are  constant,  then  PRR  is  constant  for  all  time. 

4.  For  Scenarios  1  and  2  with  k  =  1  or  A’  =  oo,  if  T,  a,  and  (3  are  constant  and 
Q  ( t )  is  unbounded,  then  the  unbounded  is  given  by 


Ptr,  (t)  ~ 


Qr 


2Qn  -  Q*  (' t ) 


(5.1) 


5.  For  Scenario  7  with  k  =  oo  and  P  (to  >  to)  =  6,  if  R  and  a  are  constant, 
then  PRR  is  monotonically  decreasing  with  radius. 

Variants  of  the  problem  were  solved  where  some  or  all  of  the  parameters  were 
assumed  fixed.  Some  problem-specific  trends  were  noticed.  Dynamic  solutions  were 
compared  to  constant  solutions.  The  effects  of  additional  warheads  was  also  exam¬ 
ined.  The  ratio  of  targets  to  false  targets  impacted  results  as  well  as  whether  the 
ratio  depended  on  sensor  footprint  location. 


5.2  Applicability 

The  results  of  this  dissertation  apply  to  search  vehicles  with  0,  1,  or  k  warheads. 
Although  the  real  motivation  for  this  research  came  from  autonomous  air  operations, 
the  results  can  be  applied  to  manned  systems,  land-based  systems,  and  water-based 
systems.  They  can  be  applied  in  the  design  phase  or  operational  phase  of  a  system. 
In  the  design  phase,  the  results  can  be  used  to  develop  system  requirements,  compare 
competing  designs,  or  verify  simulation  results.  Parameter  sensitivity  studies  can  be 
done  to  determine  the  cost  effectiveness  of  improving  sensors,  warheads,  or  vehicle 
performance.  The  mathematical  framework  enables  sound  systems  engineering.  In 
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the  operational  phase,  tactical  decisions  in  the  form  of  setting  operating  points  and 
determining  the  number  of  warheads  are  possible  using  the  results  of  this  disser¬ 
tation.  The  system  operating  characteristic  plots  of  P  (n>  n )  versus  P  (m  >  m) 
that  include  warhead  effects  give  mission  planners  a  complete  picture  of  the  trade 
space.  They  can  look  at  one  figure  and  see  the  effects  of  changing  constraint  level 
or  changing  number  of  warheads. 

5.3  Recommendations  for  Further  Research 

1.  Run  more  sample  problems.  Investigate  different  target  and  false  target 
densities.  Investigate  different  battle  spaces. 

2.  Change  from  an  airborne-based  mission  to  a  land-based  or  sea-based  mission. 
Investigate  any  problem-specific  trends. 

3.  Take  an  existing  computer  simulation  and  perform  Monte  Carlo  runs.  Com¬ 
pare  the  results  with  the  analytic  results  possible  from  this  research. 

4.  Develop  a  Matlab®  toolbox  with  a  graphical  user  interface  that  allows  the 
user  to  select  the  scenario,  the  parameters  allowed  to  vary,  the  number  of  time  steps, 
solution  tolerance,  etc. 

5.  Examine  more  sophisticated  ROC  curve  models.  Include  both  altitude  and 
velocity  effects. 

6.  Develop  expressions  for  combinations  of  uniform,  Poisson,  and  normal  dis¬ 
tributions  not  covered  by  Scenarios  1-7. 

7.  Derive  expressions  to  handle  multiple  types  of  targets  or  false  targets. 

8.  Examine  a  bivariate  normal  distribution  other  than  circular  normal. 

9.  Develop  scenarios  for  a  linear-symmetric  battle  space  where  the  ratio  of 
targets  to  false  targets  depends  on  sensor  footprint  location. 
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10.  Derive  expressions  for  Scenarios  5-7  assuming  the  search  starts  at  the  outer 
radius  and  progresses  inward. 

11.  Examine  multiple  warhead  types. 
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Appendix  A.  Probability  of  Exactly  j  Classified  Encounters 

Occurring  up  through  a  Specified  Time  or  Radius 

The  probability  of  exactly  j  classified  encounters  occurring  up  through  a  specified 
time  t  or  radius  r  is  of  interest.  A  classification  occurs  every  time  an  object  is 
encountered.  Poisson,  uniform,  and  circular  normal  distributions  are  considered. 
The  following  symbology  is  used: 

A  =  area  covered  up  through  time  t  or  radius  r 

£  =  events  involving  encounters 

C  =  events  involving  classified  encounters 

T  =  events  involving  target  attacks 

T  =  events  involving  false  target  attacks 

X  =  random  variable 

Pc  =  general  classification  probability 

Pc  can  be  Ptr,  Pftr,  (1  —  Ptr),  or  (1  —  Pftr)  depending  on  the  application.  If  only 
encounters  are  needed,  Pc  —  1.  For  each  distribution,  solutions  are  derived  assuming 
only  one  classification  probability  is  used  for  all  j  encounters.  Special  cases  include 
the  probability  of  no  target  attacks  and  the  probability  of  no  false  target  attacks.  For 
the  distributions  having  a  finite  number  of  objects  (uniform  and  circular  normal), 
solutions  are  also  derived  assuming  l  of  the  j  encounters  are  classified  with  one 
probability  and  the  other  j  —  l  encounters  are  classified  with  a  second  probability. 

A.l  Poisson  Distribution 

Let  X  be  a  random  variable  on  the  sample  space  S=  {0,1,2,...}.  If  the 
probability  function  of  X  is 

P({X  =  y})  =  e-^,  77  =  0,1,2, ...  (A.l) 

rj\ 


A-l 


where  v  >  0,  then  X  is  said  to  obey  the  Poisson  probability  law  with  parameter 
v  [39].  Equation  (A.l)  gives  the  probability  of  encountering  exactly  i]  objects  while 
searching  a  Poisson  field  of  objects.  The  Poisson  field  of  objects  is  characterized 
by  a  density  distribution  a  (or  (3 )  so  that  when  an  area  A  is  searched,  the  Poisson 
probability  law  parameter  is  v  =  aA.  This  is  tantamount  to  assuming: 

1)  The  probability  that  exactly  one  object  will  be  encountered  in  the  incremen¬ 

tal  area  A  A  is  approximately  aAA,  in  the  sense  that  it  is  equal  to  a  A  A  +  R3  (AA), 
and  -»•  0  as  AA  — ►  0. 

2)  The  probability  that  exactly  zero  objects  will  be  encountered  in  the  in¬ 

cremental  area  A  A  is  approximately  1  —  aAA,  in  the  sense  that  it  is  equal  to 
1  —  aAA  +  R-2  (AA),  and  — >  0  as  A  A  — >  0. 

3)  The  probability  that  exactly  two  or  more  objects  will  be  encountered  in  the 

incremental  area  A  A  is  equal  to  R3  (AA),  such  that  the  quotient  — >  0  as 

AA  ->  0. 

4)  If  an  area  is  subdivided  into  n  subareas  and  for  i  —  1, . . . ,  n,  Ei  denotes  the 
event  that  at  least  one  or  more  encounters  occurred  in  the  ith  subarea,  then,  for  any 
integer  n,  Ex, . . . ,  En  are  independent  events. 

Poisson  parameters  are  often  broken  down  into  the  product  of  an  encounter 
rate  /x  and  time  t.  If  an  air  vehicle  has  an  area  coverage  rate  Q,  then  //  =  aQ.  How 
to  handle  an  encounter  rate  that  varies  with  time  is  now  addressed. 

Let  /i  (r)  be  the  rate  of  encounters  at  time  r.  Divide  the  time  interval  [0,  t] 
into  n  short  time  periods  of  length  T\, ,  rn  such  that 

n 

J2Ti  =  t  (A-2) 

i—  1 

Let  the  mean  rate  of  occurrence  and  mean  probability  of  classification  in  the  ith  in¬ 
terval  be  ^  and  PCi  respectively.  From  the  fourth  assumption  above,  the  probability 
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that  exactly  jt  encounters  occur  in  the  interval  Tj,  i  —  1, . . . ,  n  is 


~lHTi 


i=l 


En 

_  ^  i= 1 


n 

i— 1 


\3i 


Ji'- 


(A.3) 


The  probability  of  no  encounters  occurring  (and  hence  no  classification)  during  [0,  t] 


is 


En 

_  i= 1  M iTi 


P  (So, a)  ~  e 

The  probability  of  one  classified  encounter  in  the  ith  interval  is 

(ViTi)1 


(A.4) 


P  (Ci.aaJ  —  e 


—  p  E?=  1  Miri  . 


1!  Pci' 


(A.5) 


The  probability  of  exactly  one  classihed  encounter  occurring  during  [0,  t]  requires 
considering  all  the  permutations  where  one  classihed  encounter  could  occur.  That  is 


P  (Si, a)  —  P  (Cl,AAi)  +  P  (Sl,AA2)  +  '  '  '  +  P  (Cl, A An) 


(A.6) 


So, 


P(Ci,, 


=  C~EL  i  HiTj(Wl)1  p  +e-E?=iWri(/f272 )lp  + 


1! 


1! 


- he 


En 

_  ^  i= 1  ^iTi 


-  E"=i  MiU  iVnTn)1 


1! 


Pr. 


1  'y  ^  iMTiPc, 


(A.  7) 


i=  1 

The  probability  of  exactly  two  classihed  encounters  occurring  during  [0,  t]  requires 
considering  all  the  permutations  where  two  classihed  encounters  could  occur.  That 
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IS 


P  {C2 ,a)  —  P  (Cl,AAi  n  CitAA2)  +  P  (Ci,AAi  n  Ci,AA3)  +  •  •  •  +  P  (Ci.AAl  H  Ci,AAn)  + 
P  (^1,AA2  n  Cl,AA3)  +  P  (Cl,AA2  ^  Cl, am)  +  '  '  '  +  P  (Cl,AA2  n  Ci;AA„)  + 

+  P  (Ci,AA„_i  n  Ci,AA„)  + 

P  (C2,AAi)  +  P  ((?2,A A2)  +  '  '  '  +  P  (C2,A. An)  (A. 8) 


So, 


P(C2A)  =  e-^i« 

En 

i= 1  ^iTi 


(A'm)1  p 
1!  C1 

(A*2  T2) 1  D 
1!  C2 


(^2^2)  p  (AbPn)  p 

1!  C2  l!  Cl 

(A^7^)  p  (/hi  Pi)  p 

1!  C3  l!  Cl 


+ 


+ 


3-  E?=i  Min  (Abi-lPt-l)  p  j^nTn)  p 

^  |  Cn  —  1  ^  |  Cn  P 

(MlTl)  p2  I  .  #  .  I  (lJ/nTn)  p2 

2!  ^Cl  '  '  2!  ^Cn 


g  X)i=  1  Mi7"! 

1 


=  -e 


En 

i=l  M*n 


lM*Tl  [2/liTiPCl  (/^27"2-Pc2  +  •  •  •  +  UnTnPcn)  + 


=  e 


2/^272 -fc2  (/i3r3-Pc3  +  '  '  '  +  A4 nJnPcn )  + 

‘  ‘  H~  2/Xn_iTrj_iPCn_1/inTrjPCr[-|- 

(/XiTiPCl)2  +  (/i2T2PC2)2  4 - h  (/VTnPcJ2] 

- £?=i Min  (Si=i  A kTiPg) 


(A.9) 


In  general,  the  probability  of  exactly  j  classified  encounters  occurring  during  [0,A]  is 


P(qA)  =  g-EIU Min  (Sk^n) 


(A.10) 
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Taking  the  limit  as  r*  — >  0,  i  =  1, . . . ,  n,  n  — ■>  oo  such  that  Yfi= 1  r*  =  t  yields 


P(Cj,a) 


(A.  11) 


Equation  (A.  11)  can  be  used  to  determine  the  probability  of  either  no  false  target 
attacks  or  no  target  attacks  occurring  up  through  time  t.  In  the  case  of  no  false 
target  attacks,  the  possible  mutually  exclusive  events  are 


no  encounters 

1  encounter  correctly  classified 

2  encounters  correctly  classified 


oo  encounters  correctly  classified 


Thus, 


P  (£o ,A  n  CM  n  .  .  .  n  Coo, a)  —  P  (Co, a)  +  P  (Cl, a)  +  •  •  •  +  P  (Coo, a) 


—  g-  So  m(t)cZt  e-  f*  n(T)dr 


fo  Pc  (t)  V  (t)  dr 

1! 


rt  v/o  Pc  (r)  A4  (r)  dr 

_ \-e~Io  M T)d ' 

oo! 


—  g-  Jo  M(T)dr 


fo  pc  (t)  h  (r)  dr 

1  +  ^ 


1! 


/o  pc  (t)  V  (r)  dr 


OO! 


_  g-  fo  f^(T)dTefo  Pc(r)ti(T)dT 


+ 


—  g“ /o  [1-Pc(T)MT)dT 


(A.12) 
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To  calculate  the  probability  of  no  false  target  attacks  for  a  straight-line  search, 
Pc  (r)  =  Pftr  (t)  and  ju  (r)  =  a  (t)  Q  (r)  in  Eq.  (A.  12)  giving 


P  (P0,A(t))  =  e~  ti[l-PFTR{T)]a{r)Q{T)dT  (A.  13) 

where  a  (r)  is  the  false  target  density.  The  Poisson  parameter  for  false  target  attacks 
(which  is  time-varying)  becomes 

V FT<yt)=  [  [1  -  PFTr  (t)]  a  (t)  Q  (r)  dr  (A.  14) 

Jo 

To  calculate  the  probability  of  no  target  attacks  for  a  straight-line  search,  Pc  (r)  = 
1  -  Ptr  (t)  and  /i(t)=/3  (r)  Q  (r)  in  Eq.  (A.  12)  giving 

P  (T0.A{t))  =  e~  JolPTR(rmr)Q(r)dr  (A.15) 


where  [3  (r)  is  the  target  density.  The  Poisson  parameter  for  target  attacks  (which 
is  time-varying)  becomes 


vr{t)=  [  [PTr  (t)}  (3  (t)  Q  (t)  dr  (A.  16) 

Jo 

Now  that  Poisson  parameters  for  false  target  and  target  attacks  exist,  more  general 
expressions  can  be  given.  The  probability  of  exactly  m  false  target  attacks  in  A  (t) 

is 

p  (A.  17) 

and  the  probability  of  exactly  n  target  attacks  in  A  ( t )  is 

P  (Tn,A(t})  =  (A. 18) 

n\ 

For  a  straight-line  search, 

A  (t)  =  f  Q  (t)  dr.  (A.  19) 

Jo 
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In  the  case  of  searching  a  circular  disc  of  radius  R  using  concentric  annuli  emanating 
from  the  origin,  it  can  be  shown  that 


uFt  (r)  =  2tt  [1  -  PFTr  (p)]  a  (p)  pdp 
Jo 


(A. 20) 


and 


isT  (r)  =  /  2t tPtr  (p)  (3  (p)  pdp 


(A.21) 


The  probability  of  exactly  m  false  target  attacks  in  A  (r)  is 


P  ( p~m,A(r ))  ^ 


_  „-vft(t)VFT  (r)' 


ml 


(A. 22) 


and  the  probability  of  exactly  n  target  attacks  in  A  (r)  is 


P  ( Pn,A(r ))  ^ 


_  „-”T(r)  Ut  (r)? 


n! 


(A. 23) 


For  a  circular  search  using  concentric  annuli  emanating  from  the  origin, 


A  (r)  =  nr2 


(A. 24) 


A. 2  Uniform  Distribution 

There  are  N  uniformly-distributed  objects  in  the  battle  space  Ap-  Divide  the 
time  period  of  length  t  into  n  short  time  periods  of  length  T\ . . ,  rn  such  that 

n 

YJu  =  t  (A-25) 

i—  1 

Let  the  mean  probability  of  classification  and  area  coverage  rate  in  the  ith  interval 
be  PCi  and  Qi  respectively.  The  incremental  area  is  given  by 

A  Ai  =  QiTi  (A. 26) 
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and  the  search  area  is  given 

n 

AS  =  Y<  AA-  (A. 27) 

i 

No  encounters  during  [0,  t]  means  the  first  object  was  not  encountered,  the  second 
object  was  not  encountered,  ...,  the  Nth  object  was  not  encountered.  These  events 
are  independent,  thus  the  probability  of  no  encounters  (and  hence  no  classification) 
involves  multiplying  the  individual  probabilities. 


P{So,a) 


4-  [  Q  (r)  dr 

ab  Jo 


N 


(A. 28) 


Let  j  —  1.  In  the  ith  interval,  either  the  first  object  or  the  second  object  or  ...  or 
the  Nth  object  is  encountered.  These  are  mutually  exclusive  events.  Therefore,  to 
calculate  the  probability  of  one  classified  encounter  in  the  ith  interval,  the  individual 
probabilities  are  summed.  That  is, 


P(CiAAi) 


PcAAi  PcAAi 

— -  +  •••  +  — - 

Ab  Ab 

' - v - ' 

N 


N- 


Pr-AAi 


A 


B 


(A. 29) 


The  probability  of  exactly  one  classihed  encounter  occurring  during  [0,t]  requires 
considering  all  the  permutations  where  one  classified  encounter  could  occur.  That  is 


P  (Cl ,a)  —  P  (Ci;AAi)  +  P  {C1AA2)  +  •  •  •  +  P  (Ci;A An)  (A. 30) 


So 


P{Ci,A) 


PgAAi  nPC2AA2 
Ab  A 


+  ---  +  N 


PCnAAn 

Ab 


N 

Ab 


Y,PcAA 


(A. 31) 


The  probability  of  exactly  two  classified  encounters  occurring  during  [0,f]  requires 
considering  all  the  permutations  where  two  classified  encounters  could  occur.  That 
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IS 


P  (C2 ,a)  —  P  (Cl,AAi  n  Ci,AA2)  +  P  (Cl,AAi  n  C\,AA3)  +  •  •  •  +  P  (Ci.AAl  H  Ci,AAn)  + 
P  (C\,AA2  n  Cl,AA3)  +  P  (C\,AA2  ^  Cl, am)  +  '  '  '  +  P  (Cl,AA2  n  Ci,AA„)  + 


+  P  (Ci,AA„_i  n  Ci,AA„)  + 

P  (C2,AA1)  +  P  {C-2,A A2)  +  '  '  '  +  P  (^2,A An)  (A. 32) 

Bayes’  rule  can  be  applied  to  the  joint  probabilities  P  (Ci.aa*  H  Ci.aaJ-  That  is, 

P  (Cl, A Ai  n  C\,AAk)  =  P  (Ci,AAi\Ci,AAk)  P  (C\,AAk)  (A. 33) 


So, 


P(C2,a) 


(N-  1) 
(N  —  1) 
(N  —  1) 
(N  —  1) 


PC2AA2  PclAAi 

Ab  Ab 
PCnAA  n  PC1AA1 

Ab  Ab 
pC3aa3  Pc2AA2 
Ab  Ab 
PCnAA  n  PC2AA2 
Ab  Ab 


+  (N-  1) 

+ 

+  (N~  1) 

+ 


PC3AA3  PclAAi 

Ab  Ab 


pC4aa4  pC2aa2 

Ab  Ab 


+  •••  + 


+  •••  + 


Pc  AAn  Pc  ,A4,_i 

(N  —  1)  Cnt _ -N  , _ —  + 


D  i—  1 


A 


B 


N\  1 


2  Al 


YJY,P^AiPckAA 


k  ~ 


B  i= 1  k= 1 


A 


B 


E  P»AA‘ 


1=1 


(A. 34) 


A-9 


The  last  step  involved  using  Fubini’s  Theorem.  In  general,  the  probability  of  exactly 
j  classified  encounters  occurring  during  [0,  t]  is 


P  ( C,.a) 


A, 


E  P«AA> 


i— 1 


(A. 35) 


Taking  the  limit  as  r*  — >  0,  i 


P  0 Cj,A ) 


1, . . . ,  n,  n  — >  oo  such  that  Y^i=i  A  =  t  yields 


/  N\ 

1  A 

G) 

X“  /  Pc{r)Q(r)dT 

L  ab  J o  J 

(A. 36) 


Equation  (A. 36)  can  be  used  with  Pc  (r)  =  1  —  Ptr  (t)  to  determine  the  probability 
of  no  target  attacks  occurring  up  through  time  t.  The  possible  mutually  exclusive 
events  are 

no  target  encounters 

1  misclassified  target  encounter  and  —  1  targets  not  encountered 

2  misclassified  target  encounters  and  N  —  2  targets  not  encountered 


N  misclassified  target  encounters 

Bayes’  rule  can  be  applied  to  calculate  the  joint  probabilities.  That  is, 

P  ( j  misclassified  encounters  D  N  —  j  objects  not  encountered)  = 

P  ( j  misclassified  encounters) N  —  j  objects  not  encountered)  x 

P  (N  —  j  objects  not  encountered)  (A. 37) 
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Therefore,  the  probability  of  no  target  attacks  occurring  up  through  time  t,  assuming 
N  targets,  is 


P(%,a)  = 


1  -  —  /  Q(r)dr 
ab  Jo 


N 


+ 


N 

£ 

j= i 


N\  r  l  r 

•  )  -7-  [!  -pTR  (t)]  Q  (r)  dr 

J  J  L  ab  Jo 


x 


1  fl 

N-j-] 

1  -  —  /  Q(t)  dr 
L  ab  Jo  J 

{ 

J 

N 

£ 

3= 0 


N\  r  l  r 

■  )  -7-  [1-Ptr  (t)]  Q  (r)  dr 
J  J  V  ab  Jo 


x 


1 


1  -  —  /  Q  (r)  dr 
Jo 


iV-j' 


1  -  -T-  [  Ptr  (t)  Q  (r)  dr 
Jo 


iV 


(A. 38) 


The  last  step  involved  using  the  binomial  theorem. 

Similarly,  the  probability  of  no  false  target  attacks  occurring  up  through  time 
t,  assuming  M  false  targets  is 


P{Fo,a) 


(A. 39) 


Equation  (A. 36)  assumes  the  same  classification  probability  is  used  for  each 
encounter.  If  j  objects  are  encountered,  suppose  l  are  classified  with  PCl  (r)  and  j  —  l 
are  classified  with  PC2  (r).  Thus, 


P  [l  classified  with  PCI  (r)  fl  j  —  l  classified  with  PC2  (r)|j  encounters]  = 


(D 

1  f* 

l 

1  P 

-7~  /  pci  (r)  Q  (t)  dr 

L ab  Jo  J 

-7~  /  PC2  (r)  Q  (r)  dr 

_ab  Jo 

(A. 40) 
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A. 3  Circular- Normal  Distribution 

There  are  N  normally-distributed  objects.  Normally-distributed  refers  to  a 
circular  normal  distribution  with  standard  deviation  a.  Thus,  the  search  area  As  is 
a  circular  disc.  Divide  the  disc  into  n  concentric  annuli  of  area  A1, . . . ,  An  such  that 

n 

Y,A*  =  A’  (A-41) 

i=i 


Each  Ai  is  calculated  using 


Ai  =  27rpiApi 


(A. 42) 


where  p*  is  the  radius  of  the  annulus  and  A p.L  is  the  thickness.  The  search  begins  at 
the  origin  of  the  disc  and  progresses  outward  such  that  pt  <  pt+ 1 ,  i  =  1, . . .  ,n  —  1. 
Let  the  mean  probability  of  classification  in  the  ith  annulus  be  PCi. 

No  encounters  during  [0,?’]  means  the  first  object  was  not  encountered,  the 
second  object  was  not  encountered,  ...,  the  Nth  object  was  not  encountered.  These 
events  are  independent,  thus  the  probability  of  no  encounters  (and  hence  no  classi¬ 
fication)  involves  multiplying  the  individual  probabilities. 


P{So,a) 


1 

27 rcr2 


P 

2a2 


pclpdd 


N 


(A. 43) 


Let  j  =  1.  In  the  ith  annulus,  either  encounter  the  first  object  or  the  second  object 
or  ...  or  the  Nth  object  is  encountered.  These  are  mutually  exclusive  events.  There¬ 
fore,  to  calculate  the  probability  of  one  encounter  in  the  ith  annulus,  the  individual 
probabilities  are  summed.  That  is, 


1  _  pi  1  _  p2i 

P{£i,aaJ  =  ^piApi^—^e  2P2  -1 - b  2'KpiApi—^ -e  ^ 


N 


AT  A  -  At  a 

=  N—e  A  pi 


(7 


(A. 44) 
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The  probability  of  exactly  one  classified  encounter  occurring  during  [0,  r]  requires 
considering  all  the  permutations  where  one  classified  encounter  could  occur.  That  is 

P  (Cl, a)  —  P  (CgAUi)  +  P  (CitAA2)  +  ■  ■  ■  +  P  (C\ ,AAn  )  (A. 45) 


So 


P(CiA)  =  N^e-£APlPcl+N^e-iAp2PC2  + 


G* 

n  2 


<7 


+  N^e  2an2  ApnPc, 


<y 


=  N\  Ai 


(A. 46) 


The  probability  of  exactly  two  classified  encounters  occurring  during  [0,r]  requires 
considering  all  the  permutations  where  two  classified  encounters  could  occur.  That 

is 


P(C2,a) 


p  (Ci,aax  n  Ci^aa2)  +  p  (CiAAl  n  CiAa3)  +  •  •  •  +  p  (C\AAx  n C\AAn )  + 

P  (Ci,aa2  n  CiAa3)  +  P  (CiAa2  hi  Ci,aa4)  +  •  •  •  +  P  (Ci,aa2  61  Ci,aa„)  + 


+  P  (Cl.AAn-!  n  C1AAn)  + 

P  (C-2,AA\  )  +  P  (' C2,AA2 )  +  •  •  •  +  P  (C-2,A An)  (A. 47) 

Bayes’  rule  can  be  applied  to  the  joint  probabilities  P  (CgAA;  H  Ci,aa*.)-  That  is, 

P  (Cl, A Ai  n  Ci,AAk)  =  P  (Ci,AAi\Ci,AAk)  P  (Ci,AAk)  (A. 48) 
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So, 


P(C2,a)  =  (N-l)^e-^Ap2PC2N^-e-^APlPCl  + 

u-  az 

(N  -  1)  ^e~-&  Ap3PC3N^e~£  APlPCl  +  •  •  •  + 
a-  crz 

(N  -  1)  ^e-S*  Ap^N^eS  APlPcl  + 
a-  a- 

(N  -  1)  ^e~A  Ap3PC3N^e~& Ap2PC2  + 
a-  aA 

(N  -  1)  ^e-^ApAPC4N^e~^Ap2PC2  +  •  •  •  + 
a-  az 

(N  -  1)  ^e-&ApnPCnN^e-iAp2PC2  + 

a-  a- 


(N  -  1)  ^  ApnPCnN^A-Ae  kpn^p  + 

a-  oz 

n  /  2  \  2 

2 ,  E 

i= 1  v 


N 


ST^  A  -At  a  p  Pk  -At  A  P 

2^2^e  ^Ap,PCi— e  ^APkPCk 


2=1  fc=l 


dfc  _-At 
O'2 


iv\  IvA  Pi  _1 £ 


Y.^e  ^^Pipa 


(A. 49) 


x  7  L*=i  J 

The  last  step  involved  using  Fubini’s  Theorem.  In  general,  the  probability  of  exactly 
j  classified  encounters  occurring  during  [0,  r]  is 


P  (Cj,A)  = 


N 


n  2 

i= 1 


(A. 50) 


Taking  the  limit  as  APl  — >  0,  i  —  1, . . . ,  n,  n  — >  oo  such  that  )T)”=1  A*  =  As  yields 


-P  (CyA(r))  — 


IV 


/•'■  p  _jA  1J 

/  p,c(p)4e  2<j2dd 

'0 


(A.51) 
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Equation  (A. 51)  can  be  used  with  Pc  (p)  =  1  —  Ptr  ( p )  to  determine  the  probability 
of  no  target  attacks  occurring  up  through  radius  r.  The  possible  mutually  exclusive 
events  are 

no  target  encounters 

1  misclassified  target  encounter  and  N  —  1  targets  not  encountered 

2  misclassified  target  encounters  and  N  —  2  targets  not  encountered 


N  misclassified  target  encounters 

Equation  (A. 37)  can  be  used  to  calculate  the  joint  probabilities.  The  probability  of 
no  target  attacks  occurring  up  through  radius  r,  assuming  N  targets,  is  then 


P  (^0,A(r))  — 


1  -  [  4e  ^ dp 
Jo  °T 


N 


+ 


N 

£ 


N' 


J  [do 


[1  -  PTr  (p)]  -4-e  2P  dp 


N 


3=0 


N 


<?i 


[1  Ptr  (p)]  —^e  2aTdp 


N 


1  - 


1  - 


r  P  -jpi 
/  AT-e  2aT  dp 

1 0 


r  P  —bl 

/  -Pre  2aT  dp 

1 0 


N-j 


N-j 


P 


1-  /  Ptr  (p)  ~Ye  2°t  dp 


On 


(A. 52) 


where  Or  is  the  standard  deviation  for  targets.  The  last  step  involved  using  the 
binomial  theorem. 

Similarly,  the  probability  of  no  false  target  attacks  occurring  up  through  radius 
r,  assuming  M  false  targets  is 


M 


p  {Po ,A(r))  =  <{  1  -  /  [1  -  Pftr  (p)]  — y— e  ^ dp 

.In  O  rrrr 


(A. 53) 


FT 


where  a  ft  is  the  standard  deviation  for  false  targets. 
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Equation  (A. 51)  assumes  the  same  classification  probability  is  used  for  each 
encounter.  If  j  objects  are  encountered,  suppose  l  are  classified  with  PCl  ( p )  and  j  —  l 
are  classified  with  PC2  ( p ).  Thus, 


P  [/  classified  with  PCl  (p)  fl  j  —  l  classified  with  PC2  (p)\j  encounters]  = 

p  P2  ll  r  rr  „  _2  lJ~l 


Pc 1  (p)  ->e  ^  dp 


L./0 


crz 


L./0 


pc2  (p)  4e  ^2dp 


(A. 54) 


A.  4  Discussion 

Now  that  expressions  for  the  probability  of  exactly  j  classified  encounters  oc¬ 
curring  up  through  a  specified  time  t  or  radius  r  exist,  the  12  elemental  probabilities 
for  each  of  the  scenarios  can  be  calculated.  This  is  done  in  Appendix  B. 
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Appendix  B.  Derivation  of  Elemental  Probabilities  for  Scenarios  1-7 

The  12  elemental  probabilities  needed  to  calculate  P  (m  =  m),  P  (n  —  n ),  P  (m  >  m), 
and  P  (n  >  h)  are  derived  for  Scenarios  1-7. 

B.l  Scenario  1 

This  scenario  is  the  specific  case  of  Scenario  3  when  N  =  1.  Hence,  the 
derivation  is  given  in  Section  B.3. 

B.2  Scenario  2 

For  this  scenario,  there  is  a  Poisson  field  of  false  targets  with  density  a  (t)  and 
a  Poisson  field  of  targets  with  density  /?(£).  The  air  vehicle  flies  along  a  straight 
path  with  area  coverage  rate  Q  (t).  In  a  Poisson  field,  there  is  an  infinite  number  of 
objects  available.  Unlike  uniform  and  normal  distributions,  events  in  a  Poisson  field 
are  independent.  Therefore,  there  is  no  need  to  account  for  every  object.  The  12 
elemental  probabilities  simplify  to 

-Pi  =  P  n  %.,<ia) 

P-2  =  P  (Pm, a) 

P3  =  P(Tn.A) 

Pi  =  P  (pm- 1,4  n  p \,dA) 

P5  =  P(PvO 

P6  =  P  (T0, A) 

P7  —  p  n  P\,dA  n  Po,Af) 

Ps  =  P  (P m,A  n  Po ,Af) 

P9  =  P  (fPn,A  H  Po,Af) 

P10  =  P  (Pm-1,4  n  P\,dA  n  Po,Af) 

Pll  =  P(P0,4S) 

Pl2  =  P  (%,AS) 
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To  calculate  joint  probabilities,  the  individual  probabilities  are  multiplied.  Using 
Eqs.  (A. 17)  and  (A. 18),  along  with  the  first  assumption  in  Section  A.l  multiplied  by 
the  corresponding  classification  probability  gives 


Pi 

=  e-»w. 

P2 

=  e~x{t) 

P3 

=  e-yW- 

Pa 

=  e~x{t)- 

P5 

=  e~x{t) 

Po 

=  e-y^ 

Pi 

=  e-y{P 

Ps 

=  e.-A'U 

Pg 

=  e-y^P 

Pi  0 

=  e-AU 

Pn 

=  e~x{P 

to 

=  e-y^P 

,  PrR(t)Q(t)j3(t)dt 

(n  —  1)! 

[*  (t)r 

m\ 

[y(tr 

n\ 

\x 

[1  —  Pftr  {t)j  Q  (t)  a  (t)  dt 


\y  (t)  l(n“1} 

,  PTR(t)Q(t)(5(t)dt 

[n  —  1)! 

[s  (*)P 

rn\ 

[y  M! 

n\ 

\x  fUl(m_1) 

_  1)1  —  Pftr  (i)]  Q  (t)  a  (t)  dt 


(B.l) 


where 

x  ( t )  =  vFT  ( t )  =  [  [1  -  Pftr  (t)]  a  (r)  Q  (r)  dr.  (B.2) 

Jo 


and 


y  (: t )  =  vT  (t) 


Ptr  (t)  (3  (r)  Q  (r)  dr 


(B.3) 
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B.3  Scenario  3 


For  this  scenario,  N  targets  are  uniformly  distributed  among  a  Poisson  field  of 
false  targets  with  density  a  (t).  The  air  vehicle  flies  along  a  straight  path  with  area 
coverage  rate  Q  (t). 

Scenario  3:  First  Elemental  Probability .  At  least  one  target  is  assumed  to  be 
encountered  in  dA.  Therefore,  at  most,  N  —  1  target  encounters  can  occur  in  A.  Let 
j  be  the  number  of  target  encounters  in  A.  For  n  —  1  target  attacks,  j  >  n  —  1. 
Applying  Bayes’  rule  to  the  equation  for  P1  gives 

Pi  =  P  n  T n  Nn n  Ti^a)  = 

P  n  T j-(n-i):A\.AfN-l-j,Af  n  %.,d.A)  P  (• J^N-l-j,Af  H  %.,d,A)  = 

P  (%i-l,A  n  T j-(n-l),A\J\fN-l-j,Af  n  Ti4a)  P  (A/jv-  1  —j,Af  \^l,dA  )P(ThdA)  (BA) 


For  a  given  j , 

P  (Tima)  =  N^PTR  (t)  Q  (i t )  dt  (B.5) 

ri  B 

using  Eq.  (A. 29), 
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using  Eq.  (A. 36),  and 


P  n  P j—(n—\),A \Af N—l—j,Af  n  Ti^a)  ~ 


J 

n  —  1 


-j-  [  Ptr  (t)  Q  (r)  dr 


iB  Jo 


n—  1 


X 


-r ~  f  Q  ir)  dT  ~  -T-  [  Ptr  (t)  Q  (t)  dr 
|_Ab  do  As  d0 


j-(n-i) 


(B.7) 


using  Eq.  (A. 40).  Define  one  new  state,  q  (t),  representing  the  probability  of  encoun¬ 
tering  a  specific  target. 

i  rf 

(B.8) 


9  It)  =  !  f  Q  M  <iT 
do 


The  state  y  (t)  is  defined  as 


V  (t)  =  J-  [  Ptr  (t)  (r)  dr 


(B.9) 


Jo 


Substituting  these  definitions  and  accounting  for  all  possible  values  of  j  gives 


N- 1  r 
j=n- 1 


^  /"  d  \  [,.OU»-l  LW  ..  UMJ-(n-l) 


d  d  \n  ~  1 


[» W]”"1  [9  (t)  -  y  (or [1  -  9  (t)]N"1_J  !•  x 


N—Ptr  (t)  Q  (t)  dt  (B.10) 

rl  B 


Scenario  3:  Second  Elemental  Probability .  Using  Eq.  (B.l), 


Po  =  e~x(t) 


x  (t)]r 


rni 


(B.ll) 


Scenario  3:  Third  Elemental  Probability.  Let  j  be  the  number  of  target  en¬ 
counters  in  A.  For  n  target  attacks,  j  >  n.  Applying  Bayes’  rule  gives 


P3  =  P  (%,A  n  T  j-n  n  A rN-j,Af)  =  P  fa, A  H  Tj  „  A'y  ,..!,)  P  {J^N-j,Af)  (B'12) 


B-4 


For  a  given  j, 


p  (rnA  n  t,  „  a”v  Q  [2/  (t)]n  [<?  (*)  -  y  {t)Y~n  (B.13) 

using  Eq.  (A. 40),  and 

P  (AT^a,)  =  P)  [1-,  («)]"-'  (B-14) 

using  Eq.  (A. 36)  where  q  (t)  is  defined  in  Eq.  (B.8)  and  y  (t)  is  defined  in  Eq.  (B.9). 
Accounting  for  all  possible  values  of  j  gives 

p3  =  £  { (’J)  Q  wrb  w -!/wrn[i- « (B.is) 

Scenario  3:  Fourth  Elemental  Probability .  Using  Eq.  (B.l), 

Pa  =  e~x (t)  [1  -  Pftr  (*)]  Q  (t)  a  (t)  dt.  (B.16) 

[m  —  1)! 

Scenario  3:  Fifth  Elemental  Probability .  Using  Eq.  (B.l), 

P5  =  e~x{t).  (B.l  7) 

Scenario  3:  Sixth  Elemental  Probability.  To  calculate  Pe,  Eq.  (B.15)  is  used 
with  n  —  0  giving 

p*  =  E  { (^)  [9  W  -  y  MF  [1  -  9  WA’}  •  (b.is) 

Applying  the  binomial  theorem  gives 

P,  =  [l-y(t)}N .  (B.19) 

B-5 


Scenario  3:  Seventh  Elemental  Probability.  Applying  Bayes’  rule  gives 


P7  —  P  (%1-1  ,A  n  P j-(ra-l),A  n  P AT-1  -j,Af  n  T\,dA)  — 

P  n  P j_(n_i)tA\T n -1- j, a f  n  P\aa)  p  (Tn-1  —j,Af  ^  7l  ,dA^) 
p  (%i-i ,a  n  T  j-(n-l),A\PN-l-j,Af  H  Pima)  P  (P  N-i-j,Af\Pi,dA)  P  (Pi  m)  •  (B.20) 


For  a  given  j, 


using  Eq.  (A. 29), 


P(Pi4a)  =  N—Ptr  (t)Q(t)  dt 
ab 


(B.21) 


P  (T N-l-j.A,\T1M)  =  { 


N-l-j 


N  -  1 
3 


Ab  Jt  PTR  ^  ®  ^  dT 

4 ~  [  Q  (T)  dr  ~  4~  f  Ptr  (r)  ^  (r)  dT 

ab  Jt  ab  Jt 


+ 4~  l  q  (r) dr  -  4~  l  Ptr  ^  dr 

AR  I n  A B  Jo 


*-B  Jo 

~4b  Jo  Q^(W+  4~  f  pTR(R)Q(R)dr 


B  Jo 
N-l 

j 


N-l-j 


[1  -q{t)-y(P)  +  y(t)] 


N-l-j 


(B.22) 


using  Eq.  (A. 36),  and 


p  (Tn-i,A  n  Pj-^AlP N-l-j, Af  n  Tim)  =  (n  3_  J  iv  (*)]n_1  [g  (*)  -  y  Wr(n”1} 

(B.23) 

using  Eq.  (A. 40)  where  q  ( t )  is  defined  in  Eq.  (B.8)  and  y  (t)  is  defined  in  Eq.  (B.9). 
Accounting  for  all  possible  values  of  j  yields 


N-l 


*=  E 

j=n- 1 


N-l\f  3 


j  J  An  -  I 


[1  -  q  (t)  -  y  (T)  +  y  (()] 


N-l-j 


N—PTR(t)Q(t)dt  (B.24) 

A  B 
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Scenario  3:  Eighth  Elemental  Probability.  Using  Eq.  (B.l), 


Po  =  e-X{T) 


X  (*)]r 


m\ 


(B.25) 


Scenario  3:  Ninth  Elemental  Probability .  Applying  Bayes’  rule  gives 


P9  —  P  {%i,a  n  T j_n  n  TN_jAf)  —  P  (%i,a  H  P ' j-n\P  N-j,Af')  P  (• P  N-j,Af )  •  (B.26) 


For  a  given  j, 


n  \  r  1  rT  ' N  1 

N-jjlAsJ  11  -Pr,0)}QO)dr 

u  L hi,  Q{T)iT~hSt  p™(r)<?(r)rfT 

+  -T-  /  Q  (r)  dr  -  [  PTR  (t)  Q  (r)  dr 

ab  Jo  ab  Jo 

"7-/  Q(r)dr+-7-  [  Ptr  (7)  Q  (r)  dr 
ab  Jo  ab  Jo 


N-j 


[l-q(t)-y(T)+y(t)]N-j  (B.27) 


using  Eq.  (A. 36),  and 


p  (rn,A  n  r,-  rY  ,f)  =  [y  (t)]n  [q (t)  -  y  (t)]j~n  (b.28) 

using  Eq.  (A. 40)  where  q  (t)  is  defined  in  Eq.  (B.8)  and  y  (t)  is  defined  in  Eq.  (B.9). 
Accounting  for  all  possible  values  of  j  yields 

p*  =  E  {  (~J)  (,))  1»  «]”  [?  (t)  -  y  MF”  [1  -  9  W  -  »  (T)  +  y  (t)]"-J}  (B.29) 
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Scenario  3:  Tenth  Elemental  Probability .  Using  Eq.  (B.l), 

Pw  =  e-*(T)  ^  '  [1  -  Pftr  (t)]  Q  ( t )  «  (t)  dt.  (B.30) 

(m  —  1) ! 

Scenario  3:  Eleventh  Elemental  Probability.  Using  Eq.  (B.l), 

Pu  =  e-x{T).  (B.31) 

Scenario  3:  Twelfth  Elemental  Probability .  To  calculate  Pi2,  Eq.  (B.29)  is  used 
with  n  —  0  giving 

=  E  {  ( ] )  1® W  -  y [i  -  i (<)  -y(T)  +  y (*)]""'}  •  (B.32) 

Applying  the  binomial  theorem  gives 

Pl2  =  [l-y(T)}N .  (B.33) 

B.f  Scenario  4 

For  this  scenario,  N  targets  and  M  false  targets  are  uniformly  distributed  in 
Ab-  The  air  vehicle  flies  along  a  straight  path  with  area  coverage  rate  Q  (t).  The 
elemental  probabilities  for  N  uniformly-distributed  targets  (Pi,  P3,  P6,  P7,  P9,  and 
P12)  were  derived  in  Section  B.3.  The  elemental  probabilities  for  M  uniformly- 
distributed  false  targets  (P2,  P4,  P5,  Ps,  P10,  and  Pu)  can  be  derived  in  the  same 
manner  replacing  j,  n,  N,  Ptr  (t),  and  y  (t)  with  i,  m,  M,  [1  —  Pftr  (£)],  and  x  (t) 
respectively,  where  x  (t)  is  defined  as 

x  (t)  =  — —  f  [1  -  PFTr  (r)]  Q  (r)  dr.  (B.34) 

ab  J  0 


B-8 


The  states  q  (t)  and  y  (t)  are  the  Scenario  3  states  defined  in  Eqs.  (B.8)  and  (B.9) 
respectively.  The  12  elemental  probabilities  become 


-t-PTr  it)  Q  (f)  dt  (B.35) 

Vi  B 

M  (  /  ft/r\  /  ■  \  ^ 

=  (i )  L  J  i*  «r  i«  w  - 1  «rm  i1  - « (*)]"■' }  (b.36) 

iV  ,  ,  -X  A 

=  E  L  J  [»  «]"  b «  -  9  (OF"  [1  -  9  (t)]N'j  (B.37) 

j=n  K  x  7  x  7  ^ 

p‘  =  E  { (M  “ *)  (m !_  J  [*  i*)]™-1  [« (*)  -  *  (i)r<m-,)  [i  - « W1M--}  ^ 

M 

——  [1  —  Pftr  (£)]  Q  (t)  dt  (B.38) 

Vi  B 

P5  =  [l-x(t)]M.  (B.39) 

P6  =  [l-y(t)]N .  (B.40) 

pi  =  E  { (iV; ') („ i  J  (9 «)]"■' tew -j/wr'”-1’ 

[1  -9(0  -!/(r)  +  !/(i)]'V_1_J}  N^bPtx  (*)<?(*)*  (B.41) 

Ps  =  5Z  1  (  ,•  )  (  *  )  ^  (t)]m  fa  (*)  -  x  WF™  t1  -  fi  (*)  —  x  (T)  +  x  {t)}M~l  \ 

l  \  /  \  /  J 

i=m  v  N  /  N  /  y 

(B.42) 

(V  ,  ,  -X  A 

a,  =  E  ,•  („  [»  wr  [9  w  -  v  wr”  (i  -  9  w  - » m + v  wi^  (b.43> 

J=n  k  X  7  X  7  ^ 


B-9 


=  E  t  { (M ,  ')  (m  1  j)  i*  mr 1  [9  m  -  *  mr  ("-,) 

[1  —  g  (S)  -x(T)  +  x(t)]M~1A  [1  -  Pftr  ( t)}Q(t)dt  (B.44) 

Pn  =  [1  -X(T)]M  .  (B.45) 

Pu  =  [l-y{T)]N.  (B.46) 

B.5  Scenario  5 

For  this  scenario,  N  targets  are  normally  distributed  among  a  Poisson  field  of 
false  targets  with  density  a  (p).  Each  target  has  a  circular  normal  distribution  with 
standard  deviation  ay.  The  search  area  As  is  a  circular  disc  of  radius  R.  The  air 
vehicle  searches  the  disc  using  concentric  annuli  of  radius  r  and  thickness  dr.  The 
search  begins  at  the  origin  of  the  disc  and  progresses  outward.  The  12  elemental 
probabilities  have  the  same  form  as  those  calculated  in  Section  B.3  for  Scenario  3, 
but  the  states  are  defined  as 


r  P  —4 

qy  (r)  =  /  —^e  2ay  dp 

J0  °y 

(B.47) 

pr 

x(r)=  [1  -  Pftr  (p)]  2npdp 

Jo 

(B.48) 

r  p  1 4 

y(r)=  /  PTr  {p)—2e  2aydp. 

Jo  ay 

(B.49) 

Also,  the  probabilities 

of  one  attack  in  dA  are 

2 

p  (Ti 4a)  =  NPtr  (r)  -^e  ^ dr 

ay 

(B.50) 

and 

P  (Pi ,dA)  =  [1  -  Pftr  (r)]  2tt ra  (r)  dr. 

(B.51) 

B-10 


With  these  substitutions, 


Pi  =  Y1  {  j  (n  J_  [y  (r)]n  1  (r)  -  y  (r)]J  (n  1}  t1  -  %  (r)]N  1  x 

_  r2 

NPTR{;r)^—e  dr  (B.52) 
ay 


p  =  e~x(r) 


x  (r)]r 
ml 


N 


j=n 


p  =  e~x(r) 


[y  (■ r)]n  [qy  (r)  -  y  ( r )]J  "  [1  -  qy  (r)]W  J 

\x 

— - — —  [1  -  PFTr  (r)]  27 rra  (r)  dr 

(■ m  —  lj! 

P5  =  e~x(r) 


Pe  =  [l-y  (r)[ 


N 


(B.53) 

(B.54) 

(B.55) 

(B.56) 

(B.57) 


u  =  E  { (N j *)  („  i  i)  i#  Mr1  fe  m  -  y  Mr'"-1’ 

^  y  _  r2 

[1  -qy{r)-y{R)  +  y{r))N-l-j\NPTR{r)^e  ^ dr  (B.58) 

ay 

p8  _  e~x(R)  Ig ^  (B.59) 

ml 

N  r  /n\  /  \  i 

Pg  =  \  (  •  )  (  3  )  [y  (r)f  idy  (r)  -  y  ( r)Y~n  [!  -  %  (' r )  -y(R)  +  y  ( r)]N~J  \ 

j=n  V\3  J  \nJ  J 

(B.60) 

N  fr 

Pi o  =  e~x{R)— - — —  [1  -  Pftr  (r)]  2nra  (r)  dr  (B.61) 

[yn  —  1)1 

Pu  =  e~x(R)  (B.62) 

P12  =  [l-y(R)]N  ■  (B.63) 


B-ll 


B.6  Scenario  6 


For  this  scenario,  N  targets  and  M  false  targets  are  normally  distributed  in 
Ab-  Each  false  target  has  a  circular  normal  distribution  with  standard  deviation  ax, 
and  each  target  has  a  circular  normal  distribution  with  standard  deviation  oy.  The 
search  area  As  is  a  circular  disc  of  radius  R.  The  air  vehicle  searches  the  disc  using 
concentric  annuli  of  radius  r  and  thickness  dr.  The  search  begins  at  the  origin  of  the 
disc  and  progresses  outward.  The  elemental  probabilities  for  N  normally-distributed 
targets  (Pi,  P3,  P6,  P7,  P9,  and  Pi2)  were  derived  in  Section  B.5.  The  elemental 
probabilities  for  M  normally-distributed  false  targets  (P2,  Pi,  P5,  P3,  P10,  arid  Pn) 
can  be  derived  in  the  same  manner  replacing  j,  n ,  N,  Ptr  (r),  qy  (r),  and  y  (?’)  with 
i,  m,  M,  [1  —  Pftr  (r)],  qx  (r),  and  x  (r)  respectively,  where 

qx  (r)  =  /  A^e  dp  (B.64) 

Jo  °x 

and 

pV  2 

x(r)=  /  [1  -  PFTr  (p)]  -^e  ^dp.  (B.65) 

Jo  °x 

The  states  qy  (r)  and  y  (r)  are  the  Scenario  5  states  defined  in  Eqs.  (B.47)  and  (B.49) 
respectively.  The  12  elemental  probabilities  become 

pi  =  j  X)  (n  t  [ V  (r)]n_1  Jh,  (r)  ~  V  (r)f (n_1)  [1  -  qy  (r)]^"1^  J  x 

F  —  r2 

NPtr  (r)  —e  dr  (B.66) 
ay 

P--  =  £  {  (f )  (!)  V  MP  fe  M  -  *  (r)]‘-m  [1  -  M]"-‘|  (B.67) 

X  =  £  {  (f )  (i)  [»  M]“  k,  (r)  -  V  MF"  [1  -  ?» 0-)]^}  (B.68) 
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p,  =  E  {  (M i  ')  (ra i)  l*  fe  M  -  1  Mf (m_1)  [1  -  9,  (r)]M-1-v}  x 

M  [1  —  Pftr  (r)l  —  e  dr  (B.69) 

°x 

P5  =  [l-x(r)]M  (B.70) 

P6  =  [l~y{r)]N  (B.71) 

p-,=  E  {('vj"1)(/1)wrr1[«,w-jwiH*_'1 

N  rp  _ 

[1  -  qy  (r)  -  2/  (P)  +  2/  (r)]^-1--7  j  iVPTK  (r)  —  e  ^ dr  (B.72) 

ay 

M  s  /  /  '  \  'j 

A  =  E  J  [*  MH  [Qx  (r)  -  X  (r)rm  [1  -  gx  (r)  -  *  (P)  +  z  ( r)]M~l 

(  \i  \m) 

i=m  /  \  /  s 

(B.73) 

N  ,  X^X  X  .X  X 

Pg  =  S  1  (  ■  )  (  J  )  ^  (r)]"  few  (r)  -  y  ir)Y~n  [!  -  (r)  -  2/  (P)  +  y  (r)]N~J  \ 

j=n  l  v  y  /  w  j 

(B.74) 

n„  =  E  {("rOL-i)  bMr'few-.wE"-1' 

i=m—  1  ^  ^  /  \  ' 

2 

[1 qx  (r)  -  x  (P)  +x  (r)]M_1_l)  M  [1  -  PFTR  (r)]  Y—e  dr  (B.75) 

Pu  =  [l-x(R)]M  (B.76) 

P12  =  [l-y(R)]N .  (B.77) 

B.7  Scenario  7 

This  scenario  is  the  specific  case  of  Scenario  5  when  N  —  1.  Hence,  the 
derivation  is  given  in  Section  B.5. 
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Appendix  C.  Nondimensionalization 

Consider  the  two  control  variables.  Probability  of  a  target  report,  Ptr,  has  no 
dimensions  and  is  between  0  and  1.  Area  coverage  rate,  Q,  has  dimensions  of  area 
per  time  and  is  between  Qmin  and  Qmax •  For  numerical  conditioning,  it  can  be 
beneficial  to  have  all  control  variables  with  no  dimensions  and  roughly  the  same 
order  of  magnitude. 

One  possible  way  to  nondimensionalize  area  coverage  rate  is  to  divide  it  by  the 
scaled  nominal  area  coverage  rate,  Qn.  Thus,  a  new  control  variable,  Q,  is  introduced 
such  that 

Q  =  A'  (C-1) 

Since  Q  <  Qn,  Q  will  be  between  0  and  1.  The  new  control  variable  must  be  inserted 
into  the  ROC  model.  Recall  the  ROC  parameter  c  was  defined  as 


(C.2) 


Thus, 


and 


(C.3) 


[1  —  Pftr  (£)] 


_ Q  CO  Ptr  (0 _ 

Q  (t)  Ptr  (t)  +  [1  —  Ptr  (£)] 


(C.4) 


The  problem  parameters  can  also  be  nondimensionalized.  If  the  false  target 
density,  a ,  is  assumed  constant,  one  can  define  the  nondimensional  parameter, 


(C.5) 


C-l 


which  represents  a  nominal  number  of  false  target  encounters.  Finally,  one  can  define 
the  nondimensional  parameter 

QnT 

K  =  — i - 

A-b 

which  represents  a  nominal  area  ratio. 

Consider  Scenario  1.  Table  C.l  lists  how  the  the  12  elemental  probabilities  and 
corresponding  state  definitions  would  be  impacted.  The  other  six  scenarios  could  be 
modified  in  a  similar  fashion.  The  final  result  is  an  optimal  control  problem  with 
nondimensional  problem  parameters  un  and  k,  nondimensional  control  variables  Ptr 
and  Q,  and  final  time  T. 


(C.6) 


Table  C.l  Scenario  1  Elemental  Probabilities  and  State  Definitions  with  Q 
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Appendix  D.  Further  Analysis  of  Scenario  1 

It  was  shown  in  Chapter  III  that  the  optimal  unconstrained  solution  was  increasing 
with  time  until  it  saturated  at  PTRmax  =  1  sometime  prior  to  the  final  time  T.  This 
strategy  of  starting  conservative  then  gradually  increasing  aggressiveness  until  the 
end,  where  an  all  out  effort  is  tried,  is  common  in  game  theory.  Without  a  constraint 
on  false  target  attacks,  the  system  has  nothing  to  lose  by  declaring  all  objects  as 
targets  near  the  end.  The  time  where  the  optimal  solution  starts  to  saturate,  tc.  is 
now  addressed. 

The  optimal  unconstrained  solution,  at  any  time,  was  shown  to  be 


Ptr  (t) 


c  -  a/Ax  (t)  acABex(P 
c  —  1 


At  tc, 

D*  U  \  _c~  a/  (A)  acABex<t°) 

*  TR  Ac)  —  1  ~  , 


Thus, 

Ax  (tc)  acABex{tc)  =  1. 


For  tc<t  <T, 


Ptr  (£)  —  !• 


The  state  and  costate  equations  when  tc  <  t  <  T  are 


(D.l) 


(D.2) 


(D.3) 


(D.4) 


x  (■ t )  =  aQ,  x  (' T )  =  x  (T) ,  tc  <  t  <  T  (D.5) 

A x(t)  =  ^  c.-AO.  Aa(T)  =  0,  tc  <  t  <  T  (D.6) 

A-b 

Integrating  the  state  equation  gives 

x  (t)  =  x  (T)  +  aQt  —  aQT ,  tc  <  t  <  T  (D.7) 


D-l 


Substituting  Eq.  (D.7)  into  the  costate  equation  gives 


A  At)  =  ~e 

A-b 


Q  „-[x(T)+aQt-aQT]  _  Q 


= —  e-m-aQT]e~aQt  t  <  t  <  T 

AB  ~ 


Integrating  Eq.  (D.8)  gives 


Ax  (t)  = 


1  r~x(T)  ,  1  —[x(T)+aQt—aQT]  f  <T 

aAB  aAB  ’  c ~  “ 


Substituting  Eq.  (D.7)  and  Eq.  (D.9)  into  Eq.  (D.3)  gives 


1  C~At)  _|_  _*_e-h(T)+oQtc-aQT]  ^  ac^Be[a:(T)+aQtc-aQT]  _  ^ 

aA  b  olA  b 


Solving  for  tc  gives 


1/1 

4  =  T  +  —  In  1 - 

aQ  \  c 


D-2 


(D.8) 

(D.9) 

(D.10) 

(D.ll) 
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